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1. Introduction 

In the standard SU (2) x U (1) gauge theory of Glashow, Salam and Weinberg predictions 
involving fermion masses and hadronic flavour changing weak transitions require a prior 
knowledge of the mass generation mechanism. The simplest method of giving mass to 
the fermions in the theory makes use of Yukawa interactions involving the doublet Higgs 
field. An as yet unconfirmed mechanism. These interactions give rise to the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix: Quarks of different flavour are mixed in the charged 
weak currents by means of an unitary matrix V . However, both the electromagnetic 
current and the weak neutral current remain flavour diagonal. Second order weak 
processes such as mixing and CP-violation are even less secure theoretically, since they 
can be affected by both beyond the Standard Model virtual contributions, as well as 
new physics direct contributions. Our present understanding of CP-violation is based on 
the three-family Kobayashi-Maskawa model of quarks, some of whose charged-current 
couplings have phases. Over the past decade, new data have enabled considerable 
reflnement of our knowledge of the parameters of this matrix V . Recent data based 
on the analysis of leptons with high center-of-mass momentum in B meson decays, 
indicate that the h u transition matrix element is nonzero. The complex phase of 
this matrix element is very important for the successful description of CP-violation 
within the framework of the CKM matrix. Results of experiments searching for the 
difference between CP violating decays of kaons to pairs of neutral and charged pions 
have been presented by FNAL and CERN, which support our understanding of CP 
violation through the CKM matrix. The top quark enters into several constraints on 
CKM parameters through loop diagrams, so that such an analysis necessarily implies 
a favored range of top quark masses. Over the past decade or so, many methods have 
been proposed for obtaining the three interior angles of the unitarity triangle of the 
matrix V, «, (3 and 7. Presently these CP phases are being measured in a variety of 
experiments at B-factories, KEK-B, SLAC-B, HERA-B and will be measured at LHC-B 
and B-TcV. As always, the hope is that these measurements will reveal the presence of 
physics beyond the Standard Model. 

The small visible branching ratio of B decays to CP eigenstates, C(10~^), requires 
a large number of B mesons to be produced in order to study CP violation in enough 
detail. The two complementary methods are e+e~ colliders tuned to the T(4S) resonance 
or high energy hadron machines where the hb cross-section is large. 

Current e+e~ colliders have achieved peak luminosities of ~ 3 x lO^^cm^^s"^ 
producing hb pairs at the rate of ~ 1 Hz, although raw data rates are considerably 
higher. The B mesons are produced in a coherent state and it is necessary to measure 
the time separation of both decay vertices to measure CP asymmetry in B^B^ mixing. 
To this end both, both PEP-H and KEK use asymmetric beam energies to boost the 
distance between the decay vertices. Yet with typical separations of only ~ 250 //m, the 
detectors used to resolve the vertices must still be as close as possible to the beam fine 
to achieve suitable vertex accuracies of ~ 50-100 //m. Symmetric energy machines can 
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use the decays of charged B mesons for CP investigations. In all cases, the known initial 
beam energies, even in the presence of initial state radiation, is an important constraint, 
improving reconstructed B mass resolutions by an order of magnitude. The B^ system 
can be investigated by moving to the T(5S') resonance, although with a large reduction 
in cross-section. 

The hadron colliders have the advantage of much higher production rates, C(IO^) 
Hz at the Tevatron and 10 times greater again at the LHC as well as producing both 
B and B5 mesons without then need for altering beam energies. As the centre of mass 
energy increases, the ratio of the hb to inelastic cross-section increases. The challenge is 
however to cope with the very high rates of background events and the large numbers 
of tracks in all events. These very high interaction rates require the use of sophisticated 
triggers operating at high rates. The beam crossing rate at the LHC for example is 
40MHz while B decays that are interesting occur at a rate of a few Hz. 

The accurate reconstruction of decay vertices and the ability to cleanly identify 
hadrons are major design requirements of both current and future detectors. Good 
vertex resolution requires the use of high precision vertex detectors very close to the 
beam pipe while the high particle fluxes place stringent requirements on the radiation 
hardness of the devices used. The large number of B hadrons produced means that, 
particularly at hadron machines, it becomes possible to make use of B decays with 
small cross sections where good differentiation between charged pions and kaons is 
required. High precision measurements in these final states will be the forte of the 
future experiments. This has prompted the development of high precision Chcrcnkov 
counters. The use of final states containing neutral particles requires the use of finely 
segmented calorimeters and will be particularly difficult at hadron coUiders given the 
high backgrounds that will be encountered. 

A complete review on the present status of CP measurements and their future 
prospects is beyond the scope of the present report. In this report we summarize the 
discussions in the working group on CP violation at the UK phenomenology workshop 
2000 in Durham. In the following we give a short outline of the various topics which 
were addressed. In Section 2 we describe the measurements of sin2j3 by the BaBar and 
the Belle experiment and discuss some prospects for the future. We also summarize 
CP violation measurements from B decays observed at CDF during the last Tevatron 
run (Run-I). Prospects for measuring CP violation in the period between March 2001 
and March 2003 (Run-II) at the upgraded CDF-II detector are also summarized. It is 
expected that the current generation of experiments will make the first observation of 
CP violation in the B system. We also concentrate on the anticipated performance of 
the next generation experiments, namely LHCb at CERN and BTeV at Fermilab, which 
will measure CP violating observables with extremely high precision, thereby thoroughly 
testing the Standard Model description of CP violation and searching for new physics 
beyond. In section 3 we focus on the infiuence of new physics on CKM phenomenology 
and on CP-violating observables. We also review various methods for determining the 
weak phase angle 7. Section 4 is devoted to the CP violation parameter e'/e within 
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and beyond the SM. In Section 5 discussions on specific -B^— decays and on the width 
difference in the B^-Bd system are included. In Section 6 we give an overview of the B- 
Physics trigger strategy for Run 2 at the Tevatron, with an emphasis on CP- Violation. 
Also the LHCb trigger strategy is reviewed. In addition a method to separate B events 
from continuum background in BaBar is presented. In Section 7 the phenomenological 
impact of the QCD-improved factorization approach is discussed while Section 8 deals 
with statistical issues relevant to heavy flavour physics including confidence level and 
the new technique due to Feldman an Cousins. 

2. What ought we to be measuring? 

2.1. The BABAR measurement of sin2/3 and its future prospects 

James Weatherall, Univ. Manchester (representing the BABAR Collaboration) 

2.1.1. Introduction The BABAR experiment consists of an asymmetric electron- 
positron collider operating at the T(4S) resonance. More details on the detector can be 
found in The aim is to overconstrain the unitarity triangle by measuring its sides 
and angles. The analysis reviewed here measures sin2/3 by studying time-dependent CP 
violating asymmetries in B^ J/ipK^ and B° ip'Kg decays. 

2.1.2. Overview of the sin2/3 analysis There are five main parts to measuring the 
CP violating asymmetry: 

• Selection of signal CP events 

• Measurement of the distance Az between the two B° decay vertices along the T(4S) 
boost axis 

• Determination of the flavour of the tag-side B 

• Measurement of dilution factors for the different tagging categories 

• Extraction of sm2(3 via an unbinned maximum likelihood fit 

Event Selection The sample used for the analysis is 9.8 fb~^ of data recorded 
between January and July 2000 of which 0.8 fb^^ was recorded 40 MeV below the 
T(4S) resonance. Particle identification uses mainly the Csl calorimeter for electrons, 
the Instrumented Flux Return for muons and the DIRC for kaons. Extra information 
is provided by dE/dx measured in the tracking system. The selection for the CP events 
proceeds as follows. Pairs of electrons or muons coming from a common vertex are 
combined to form J/ip and ip' candidates. The ip' is also reconstructed from its decay 
into J/ip TT+TT^. The Ks candidates are made from either a pair of charged tracks or 
a pair of tt^ candidates. In addition there are various event shape and topological cuts 
designed to reduce continuum and BB background. Full details of the selection can be 
found in . The final event sample is shown in figure |T[ 
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There are two other B decay samples. One consists of fully reconstructed 
semileptonic (5° D*^l^vi) and hadronic (5° D^*^~ -k^ , D^*^~ p'^ , D^*^' a^) decays 

(*)0 

as well as a control sample of D tt"'' events. The selection of this sample is 

described in and 0]. The other is a charmonium control sample containing fully 
reconstructed neutral or charged B candidates in two-body decay modes with a J/V' in 
the final state (e.g. B+ Jjip K+,B° J/i) K*^{K*^ fsT+vr-)). 
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Figure 1. CP signal event distributions for J / t/j K s (tt^ 7^ ) (left), J/iPKs{tt'^tt'^) 
(middle) and ip' Ks{Tr^TT^) (right). 



Measuring Az The time-dependent decay rate for the Bqp is given by 

/±(r, Anid, Vsm2(3, t) = ^Te"^!*! [1 ± Vsm2(3 x sinAm^t] (1) 

where the + or - sign indicates whether the Btag was tagged as a B^ oi B^ respectively. 
The dilution factor T> is given by P = 1 — 2w, where w is the mistag fraction (the 
probability that the Btag is identified incorrectly). To account for finite detector 
resolution, the time distribution must be convoluted with a resolution function: 

2 

n{Az; a)=y] — ^ exp (-(A^ - 5,)V2a,') , (2) 

which is just the sum of two Gaussians where the 6i and ai are the normalizations, 
biases and widths of the distributions. In practice two scale factors Si and S2 are 
introduced such that = Si x a^^t where (Taj is an event-by-event calculated error on 
At. They take account of underestimating the uncertainty on At due to effects such as 
hard scattering and possible underestimation of the amount of material traversed by the 
particles. The resolution function parameters are obtained from a maximum likelihood 
fit to the hadronic B^ sample and are shown in table |1|. The parameter represents the 
width of a third Gaussian component, included to accommodate a small (~1%) fraction 
of events which have very large values of Az, mostly caused by vertex reconstruction 
problems. This Gaussian is unbiased with a fixed width of 8 ps. Further details can be 
found in [§. 
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Parameter 


Value 


6i (ps) 


-0.20 ±0.06 


from fit 




1.33 ±0.14 


from fit 


(%) 


1.6 ±0.6 


from fit 


/i (%) 


75 


fixed 


^2 (ps) 





fixed 


S2 


2.1 


fixed 



Table 1. Resolution function parameters. Those, labeled 'from fit' are measured from 
data and those marked 'fixed' are determined from Monte Carlo. 

B flavour tagging Each event with a CP candidate is assigned a S° or S° tag if it 
satisfies the criteria for one of the several tagging categories. The figure of merit for 
each tagging category is the effective tagging efficiency Qi — ei{l — 2wiY where is 
the fraction of events assigned to category i and Wi is the probability of mis-tagging an 
event in this category. The statistical error on sin2/9 is proportional to l/\/Q where 
Q — 'YliiQi- There are five tagging categories: Electron, Muon, Kaon, NTl and NT2. 

The first three require the presence of a fast lepton and/or one or more charged 
kaons in the event and depend on the correlation between the charge of a primary lepton 
or kaon and the flavour of the h quark. If an event is not assigned to either the Electron 
or Muon categories, it is assigned to the Kaon category if the sum of the charges of all 
the identified kaons in the event is different from zero. If both lepton and kaon tags are 
available but inconsistent the event is rejected from both categories. 

NTl and NT2 arc categories from a neural network algorithm, this approach 
being motivated by the potential fiavour-tagging power carried by non-identified leptons 
and kaons, correlations between leptons and kaons and more generally the momentum 
spectrum of charged particles from B meson decays. The output of the neural network 
tagger xnt can be mapped onto the interval [-1,1] with xnt < representing a 5° tag 
and xnt > a, B tag. Events with \xnt\ > 0.5 are classified in the NTl category and 
events with 0.2 < \xnt\ < 0.5 in the NT2 category. Events with \xnt\ < 0.2 are excluded 
from the final analysis sample. 

Measurement of tagging performance The effective tagging efficiencies and 
mistag fractions for all the categories are measured from data using a maximum 
likelihood fit to the time distributions of the 5" hadronic event sample. The procedure 
uses events which have one B fully reconstructed in a flavour cigcnstate mode. The 
tagging algorithms are then applied to the rest of the event, which represents the 
potential B^ag- Events are classified as mixed or unmixed depending on whether the 
Btag is tagged with the same or opposite fiavour as the Bqp- One can express the 
time-integrated fraction of mixed events x as a function of the B'^B mixing probability, 
X — Xd+{^~ '^Xcd'w where Xd — i^d/ (1 + with Xd — Am^/r. Thus an experimental 
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Tagging category 


e(%) 


w (%) 


Q{%) 


Lepton 


11.2 ±0.5 


9.6 ±1.7 ±1.3 


7.3 ±0.3 


Kaon 


36.7 ±0.9 


19.7± 1.3± 1.1 


13.5 ±0.3 


NTl 


11.7±0.5 


16.7 ±2.2 ±2.0 


5.2 ±0.2 


NT2 


16.6 ±0.6 


33.1 ±2.1 ±2.1 


1.9 ±0.1 


all 


76.7 ±0.5 




27.9 ±0.5 



Table 2. Tagging performance as measured from data. 



value of the mistag fraction w can be deduced from the data. 

A more accurate estimate of w comes from a time-dependent analysis of the fraction 
of mixed events. The mixing probability is smallest at low At so that this region is 
governed by the mistag fraction. Figure shows the fraction of mixed events versus At. 
The resultant tagging performances are shown in table |^. 



1 .00 I 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 




0.00 I ' ' ' 1 ' ' ' 1 ' ' ' 1 ' ' ' 1 

2 4 6 8 

At (ps) 

Figure 2. The fraction of mixed events as a function of |Ai| for data events in the 
hadronic sample for neutral B mesons (full squares) and charged B mesons (open 
circles). The dot-dashed line at tcut = 2.5 ps indicates the bin boundary for the 
time-integrated single-bin method. 

Extracting sin2/3 A blind analysis technique was adopted for the extraction of sin2/3 
to eliminate possible experimenter bias. The technique hides both the result of the 
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Source of uncertainty 


Uncertainty on sin2/3 


uncertainty on 


0.002 


uncertainty on Am^ 


0.015 


uncertainty on A2; resolution for CP sample 


0.019 


uncertainty on time-resolution bias for CP sample 


0.047 


uncertainty on measurement of mistag fractions 


0.053 


different mistag fractions for CP and non-CP samples 


0.050 


different mistag fractions for and 


0.005 


background in CP sample 


0.015 


total systematic error 


0.091 



Table 3. Summary of systematic uncertainties. The different contributions are added 
in quadrature. 



likelihood fit and the visual CP asymmetry in the At distribution. This method allows 
systematic studies to be performed while keeping the numerical value of sin2/3 hidden. 

Possible systematic effects due to uncertainty in the input parameters to the 
fit, incomplete knowledge of the time resolution function, uncertainties in the mistag 
fractions and possible limitations in the analysis procedure were all studied. Details can 
be found in [Q. The systematic errors are summarized in table ^. 

Results and checks The maximum likelihood fit for sin2/?, using the full tagged 
sample of 120 5° J/i^K^ and ip'K^^ events yields: 

sin2/5 = 0.12 ± 0.37 (stat) ± 0.09 (syst) (preliminary). (3) 

The log likelihood is shown as a function of sm2(3 in figure |^. The raw asymmetry as a 
function of At is shown in figure § 

The probability of obtaining a statistical uncertainty of 0.37 is estimated by 
generating a large number of toy Monte Carlo experiments with the same number of 
tagged CP events as in the data sample. The errors are distributed around 0.32 with a 
standard deviation of 0.03, meaning that the probability of obtaining a larger statistical 
error that the one observed is 5%. From a large number of full Monte Carlo simulated 
experiments, we estimate that the probability of finding a lower value of the likelihood 
than the one observed is 20%. 

Several cross-checks are performed to validate the main analysis. The charmonium 
and fully-reconstructed hadronic control samples are composed of events that should 
exhibit no time-dependent asymmetry. These events are fitted in the same way as the 
signal CP events to extract an "apparent CP asymmetry". The results are shown in 
table |. 
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Figure 3. Variation of the log likelihood as a function of sin2/3. The two horizontal 
dashed lines indicate changes in the log-likelihood corresponding to one and two 
statistical standard deviations. 



Sample 


Apparent CP asymmetry 


Hadronic charged B decays 


0.03 ±0.07 


Hadronic neutral B decays 


-0.01 ±0.08 


J/i}K+ 


0.13 ±0.14 


J/^IjK*\K*° K+TT-) 


0.49 ±0.26 



Table 4. Summary of systematic uncertainties. The different contributions are added 
in quadrature. 



Constraints on the unitarity triangle The Unitarity Triangle in the {'p,r]) plane 
is shown in figure |^. The two solutions corresponding to the measured central value are 
shown as straight lines. The cross-hatched regions correspond to one and two times the 
one-standard-deviation experimental uncertainty. The ellipses represent regions allowed 
by all other measurements that constrain the triangle. They are shown for a variety of 
choices of theoretical parameters. More details can be found in p . 

2.1.3. Future prospects The preceding pages describe only a preliminary 
measurement of sin2/? by the BABAR experiment. More data will allow extra channels 
to be included in the final fit as well as providing more events for the currently used 
decay modes. The new channels will bring extra experimental and theoretical challenges 
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Figure 4. The raw B° - B° asymmetry {Nbo - 7V^)/(A^ijo + N^o). The time- 
dependent asymmetry is represented by a solid curve for the central value of sin2/3, 
and by two dotted curves for the values at plus and minus one statistical standard 
deviation from the central value. The curves are not centered at (0,0) because the CP 
sample contains an unequal number of B° and B events (70 versus 50 B ). The 
between the binned asymmetry and the result of the maximum likelihood fit is 9.2 
for 7 degrees of freedom. 

with them. Such present and future issues are discussed in the sections that follow. 

Available modes The B decay modes that have been used to measure sin2/9 up to 
now are clean in that they are vector-scalar, b ccs transitions which have no significant 
pollution from penguin diagrams. The next step is to add vector-vector modes such as 

— s> J/ipK*. These modes require an angular analysis of the vector meson decay 
products, due to the different partial waves and therefore admixture of CP odd and CP 
even that is present in the final state. Such an angular analysis has already yielded 
preliminary results for the J/ipK* modes. Once one has measured the polarizations 
in these modes, they are as clean, theoretically, as the vector-scalar modes. Another 
obvious addition is B° — > J/ipKi decays where the challenge here is to understand the 
background well enough to make the channel feasible. Work is ongoing in this area. 

A different kind of difficulty is presented by channels with a significant degree of 
penguin contamination, such as 6 — ccd scalar-scalar modes (e.g. — > D^D~). Here 
the fit must take into account the fact that the true value of sin2/5 is shifted by an 
amount proportional to the ratio of tree to penguin contributions. This ratio is model 
dependent and subject to large theoretical uncertainties. 
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Figure 5. Present constraints on tlie position of the apex of the Unitarity Triangle 
with tlie BABAR result indicated by the cross-hatched regions. 

Finally, modes such as D*D* and B — ^ J/ipp^ which are vector- vector, 

h — > ccd transitions face the theoretical challenges of the penguin contaminated modes 
described above, as well as requiring an angular analysis to solve the vector-vector CP 
admixture problem. 

Experimental considerations There are also experimental analysis issues which 
need to be resolved or studied in greater depth in the future. The tagging algorithms that 
BABAR uses should be developed and extended to include extra tagging categories such 
as the using the soft pion from D* decays and incorporating leptons at an intermediate 
momentum (i.e. from a cascade). It would be useful to take account of correlations 
within an event, such as when two different tagging categories report an answer. This 
can give more information about the event if the correlations are well understood. There 
is also an open question when it comes to measuring the tagging performance from the 
hadronic or semileptonic B decay samples. One then needs to be absolutely sure that 
using exactly the same numbers for the CP signal event sample is a valid thing to do. 

The measurement of Az is another crucial part of the analysis and it is important 
that the errors and biases to this distribution are understood. The distribution tends 
to be biased by the decays of particles which fly significantly from the original B decay 
vertex, such as D^s. These can be rejected by looking explicitly for cascade decays. 
The parameterization of the resolution function incorporates detector effects such as 
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misalignments and electronics readout effects. All contributions to the width should be 
studied in order to fully understand the error on A2;. 

Backgrounds to the various CP modes can also be a problem. The channels vary in 
terms of how much background they experience and this background can be particularly 
dangerous if it has a significant structure in Az. For charmonium channels, much of the 
background comes from events containing a real J/ip. In that case, one needs to study 
exactly which modes contribute and what their shape is in the final distributions (if they 
cannot be removed otherwise). Non- resonant backgrounds to vector- vector modes such 
as the J/ipK^Ti^ contribution to J /il)K*^{K*^ — > K^ti^) are in principle dangerous since 
they can have CP violating properties but no angular structure. However, the branching 
ratios for these non-resonant modes are typically poorly known and consistent with zero 
making it difficult to simulate them in the correct proportions. 



Study of statistical error It seems anomalous that both BABAR and Belle record 
higher statistical errors than one would expect. The fitting procedure is, and continues 
to be a vigorously studied part of the analysis as we need to be certain that the likelihood 
function is of exactly the correct form for the final fit. 



2.1.4- Conclusions A preliminary measurement of sin2/3 by BABAR has been 
presented. The errors on the final result make it difficult to express its significance in 
terms of constraints on the Unitarity Triangle. However, results based on a much larger 
data sample (~20 fb~^) will soon be available. Combined with a better understanding of 
systematic effects, this should make the next measurement of sin2/3 even more interesting 
than the current one. It is also expected that other CP modes will soon be available 
for analysis including — > J /ipK*^{K*^ Kstt^) and B^ J/ipKi. The larger data 
sample with additional CP modes should yield a value of sin2/3 for which the statistical 
and systematic errors are about one-half of their current values. 
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2.2. The first results from Belle 

Yoshihito Iwasaki, KEK-IPNS (representing the Belle Collaboration) 



2.2.1. Introduction Observation of CP violation in the B meson system is one of 
the most exciting physics targets dX a, B factory experiment. In the standard model. 
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CP violation is a natural consequence of the complex phase of quark mixing in the 
weak interaction as described by the Kobayashi-Maskawa matrix|I[] . This phase can be 
detected in physics processes where amplitudes with different KM phase interfere. 

In the decay of the neutral B meson to a CP eigenstate, at least two amplitudes 
fcp) and A{B° ^ B° ^ fcp)) exist due to mixing. These amplitudes 

interfere. The time dependent CPasymmetry in the decay rate of B^ and B^, acp{t), 
can be written asp| 

... nB'jt) ^ fcp) - n'B'it) ^ fcp) ... 

acpyt) = =0 (4) 

T{B'^{t)^fcp) + T{B\t)^fcp) 

= r/cp sin 201 ■ sin( Ami), (5) 

where t is the proper decay time of the B, rjcp is the CP eigenvalue of the final state 
fcp, Am is the mass difference between the two B^ mass eigenstates, and 0i (also 
known as (3) is one of the three angles of the CKM unitarity triangle formed by the d 
and b quark, 

01 = 7r-arg(— — — ). (6) 

In a -B factory experiment, a BB pair is created from T45 decay in a coherent 
quantum state. At the decay of one B, the other B oscillates starting with the opposite 
flavour of the B. Experimentally we measure At = tcp — ttag instead of t, where tcp is 
the decay time of the neutral B meson decaying to a CP eigenstate(i?c'p), and ttag is the 
decay time of the other B{Btag)- The flavour of Btag specifies the flavour of Bcp at the 
start of the B^-B^ mixing. To extract sin20i , we measure the proper time difference 
distributions instead of acp{t) : 

^ = 7r^exp(--^^)(l - r7cpsin20i sin(AmrfAt)) (7) 
at It go r go 

2.2.2. KEKB accelerator and Belle detector KEKB is an asymmetric energy 
e'^e~ collider that produces a boosted T45 in the laboratory frame. The beam energies 
of e"*" and e~ are 3.5 and 8.0 GeV, respectively. The boost factor of the T45 is 0.425. 
The beam size at the interaction point is 2/im and 100/xm in the vertical and horizontal 
direction, respectively. 

The Belle detector is located at the interaction point of KEKB in the Tsukuba 
experimental hall. Construction was started in 1994 and completed in 1998. Belle is a 
general-purpose detector with a Silicon Vertex Detector (SVD), a Central Drift Chamber 
(CDC), an Aerogel Cherenkov Counter (ACC), a Time Of Flight scintillation counter 
(TOP), an Electro-magnetic Calorimeter (ECL), a solenoid magnet, and a Ki^ muon 
catcher (KLM). Because the beam energies are asymmetric, the detector shape is also 
asymmetric in the beam direction in order to cover a large solid angle in the T4S' rest 
frame. 

Charged tracks are reconstructed from CDC and SVD hits as the particle spirals 
in the solenoidal 1.5 Tesla magnetic field. The transverse momentum resolution is 
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^ptlVt = O.OOlQpt ©0.0034. The impact parameter resolution in the plane perpendicular 
to the beam axis is = 21 sin^^/^e where p is the momentum in GeV and (3 is the 
velocity. The impact parameter resolution along the beam direction is cr^ = 39® ^^ sin^/^g 
/im. 

Photons are reconstructed from the energy deposited in the ECL with a resolution 
of aE/E = 0.013 © 0.0007/^ © 0.008/E^/^ where E is the measured energy in GeV. 
Kaons are identified using probabilities calculated from dE/dx measured by CDC, TOP, 
and hits in ACC. The time resolution of the TOP is 95 ps. The refractive index used 
by the ACC is chosen to provide good vr/K separation for 1.5 < p < 3.5 GeV. The 
efficiency is ~ 80% and the fake rate is ~ 10% for momentum up to 3.5 GeV. Electron 
identification is done using dE/dx, ACC hits, and energy deposited in the ECL. The 
efficiency is above 90% for p > 1.0 GeV. The fake rate is below 0.5%. Muon identification 
is done with hits in the KLM. The efficiency is above 90% and the fake rate is below 
2%. 



2.2.3. Event selection of CP eigenstates The CP eigenstates we search for are 
charmonium plus as listed in table p.2.3| . All decay modes include h — > ccs transitions 
except for J /ip{l^l~) + 7r°, which is a 6 ^ ccd transition. 

Candidate J/ip mesons are formed from pairs of oppositely charged tracks where 
at least one track is positively identified as a lepton (/i or e) and the other is consistent 
with a lepton. Por the J/ip ^ e+e" channel we also include 7's within 0.05 radians 
of the electron direction to recover events with initial state radiation from the e. We 
require the invariant mass of a lepton pair to be —0.006 GeV < M^+^- — Mj/^ < 0.036 
GeV and -0.150 GeV < M^+^- - Mjj^ < 0.036 GeV for and e+e" pairs, 

respectively. Por tp' —>■ we require the invariant mass of the lepton pair to be 
-60 MeV < M^+^- - M^, < 36 MeV and -150 MeV < M^+e- - M^, < 36 MeV for 
fi~^fi~ and e~^e~ pairs, respectively. Por ip' J/ip-K^Ti^ we combine a J/ip candidate 
and a tt+tt^, where the invariant mass of the pion pair is required to be greater than 400 
MeV. We require the mass difference of iIj'{1^1^7i^7i^) and J/ip to be between 0.58 GeV 
and 0.60 GeV. Por Xci J/i'li we combine a 7 and J/ip where the 7 is not consistent 
with forming 7r° when combined with any other 7. 

Por Ks — > 7r+7r~ we select a pair of oppositely charged tracks where the closest 
distance of two tracks in the z coordinate is consistent with zero. The invariant mass 
of the candidate is required to be within ±3(T (±16 MeV) of the Ks mass peak. Por 
Ks —>■ 7r°7r° we use pairs of vr^'s where 7r° is reconstructed from two 7's. We require the 
invariant mass of Ks and 7r° to be within the range 0.3 GeV to 1.0 GeV and 118 MeV 
to 150 MeV, respectively. 

Por J /ipTT^ we select a high momentum 7r° reconstructed by two 7's where the energy 
of each 7 is required to be greater than 100 MeV. The invariant mass requirement is 
identical to the vr" requirements in Ks reconstruction. 

Por selection of B candidates in all modes except for B —>■ J/ipKi, we use the energy 
difference AE = £"5- £^cm/2 and the beam-constrained mass Mbeam = a/ {Ecm/'^Y — Pb- 
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Decay mode CP AE (MeV) Nbg Ntagged 



5° - 


J/^Pil+l-) + Ksin+TT-) 


-1 


±40 


70 


3.4 


40 






-1 


±100 


4 


0.3 


4 




^ij'{m-) + Ks{7i+7i-) 


-1 


±40 


5 


0.2 


2 




f{J/ij{l+l-)n+n-) + Ks{n+n-) 


-1 


±40 


8 


0.6 


3 




-^Xci{J/Hl^l-h) + Ks{'K+n-) 


-1 


+30 
-40 


5 


0.75 


3 




J/ijil+l-) + Kl 


+ 1 




102 


56 


42 




J/^{l+l-)+Tl^ 


+ 1 


+50 
-100 


10 


1 


4 



Table 5. The decay modes of B going to CP eigenstate. I represents e or /i. 



In figure the scatter plot of /S.E versus Mheam is shown for B J /ipKsiT^^T^')- We 
define the signal region to be \A.E\ < 40 MeV and \Mi,eam — {Mbeam) I < 10 MeV where 
{Mbeam) is the mean of the observed Mbeam- The signal region in AE is varied depending 
on the decay mode (see table |2.2.3D . However, the signal region in Mbeam is the same 
as that for B J/ipKs because the error on Mbeam is dominated by the beam energy 
spread. 

For B J/ipKi we use tighter cuts for J/ip reconstruction by requiring positive 
identification for both leptons and the momentum of the J/iIj{Pj^^) in the CMS frame 
to be 1.42 < Pj/^ < 2.00 GeV. We reconstruct the momentum of the from the 
momentum of the J/ip with the assumption of a two body decay B J/ipKi- We also 
require associated KLM hits in the direction of the momentum. To select signal 
events we require the momentum of the B in the CMS, P^, to be in the range 200 MeV 
^ Pb ^ 450 MeV. A true candidate should peak around 340 MeV corresponding to the 
initial momentum of the B from the T45. In figure |^, the P^ distribution for J/ipKi 
candidates is shown with the expectation obtained by a full MC simulation study. In 
the signal region, we have 102 J/ ipK^ candidates where we expect 8 background events 
from B J/iIjKl7c^ (a mixture of CP + 1 and CP — 1) and 48 background events from 
other sources. 



2.2.4- Flavour tagging To determine the flavour {B^ or B^) of the Bqp candidates, 
we partially reconstruct the other Btag in the event. The flavour of the Bcp should 
be opposite to that of the Btag at the decay time of the Btag{At = 0). We apply four 
methods sequentially. 

(i) the charge of a high momentum lepton (P^* > 1.1 GeV) : This tags the flavour via its 
primary leptonic decay. We assign Bcp = B^{B^) if the charge is positive(negative). 

(ii) sum of the charge for positively identified kaons : This relies on s quark pop- 
up in cascade decays. We assign Bcp = B^{B^) if the sum of the charge is 
positive (negative) . 
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Figure 6. (a) The scatter plot of AiJ versus Mbeam for — > J/iIjKs candidates, (b) 
The projection in Ai? with the cut \Mheam — Mgo] < 0.01 GeV. (c) The projection in 
Mbeam with the cut |A£;| < 0.04. 



(iii) the charge of a medium momentum lepton : This is similar to method (1) but for 
the case where there is large missing momentum in the CMS(P*jj^^,). Wc require 
the lepton momentum to be in the range 0.6 < Pj* < 1.1 GcV and the missing 
momentum in the CMS should satisfy -P;*+-P^iss — ^-0 GeV. The flavour assignment 
is identical to method (1). 

(iv) the charge of slow pion : This tags the flavour of D* coming from Btag- We require 
the momentum of the slow tt to be less than 200 MeV. We assign Bcp = B^[B^) if 
the pion charge is negative (positive). 

The flavour tagging efficiency, e, and the wrong tag fraction, uj, are measured from 
data with self-tagging decay modes. We exclusively reconstruct B^ — > D^*'>~l^v, and 
apply the flavour tagging methods for the rest of the event. Because of B^-B^ mixing. 
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Figure 7. The distribution for B J /'ipKi, candidates. The upper solid line is the 
fit of a sum of the signal and background expected from MC simulation studies. The 
lower solid line is the total background. The dotted line is the background component 
from sources other than B — > J/tpKLir^. 



the probability to find the opposite or same flavour for the exclusively reconstructed B 
and the result of the flavour tagging is 

Poppositei^t) (X 1 + (1 - cos(AmrfAt) (8) 
Psamei^t) ocl - (1 - 2cj) cos(AmdAt). (9) 

Then, u can be extracted from the amplitude of the B^-B^ mixing 

A = - ^^-"^^ = {l-2uj) cos( Am^At) . (10) 

Popposite ~l~ Psame 

The vertex position of the D'^*^^l^u is determined by requiring the / and D^*^ 
to form a common vertex. The determination of the vertex position of the Btag is 
described in the next section. At is calculated from the difference of the two vertices in 
the z direction. We perform an unbinned maximum likelihood flt to the amplitude of 
the B^-B^ mixing to obtain u with allowed to be free. In table p.2.4| , e and u are 



summarized. We obtain Am^ = 0.488 ± 0.026 ps ^, which is consistent with the world 



averageP]. In table 2.2.3, the number of tagged events for CP eigenstates are listed. 



We flnd 98 tagged events in total. 

2.2.5. Proper-time difference The proper-time difference is estimated from the 
difference of z coordinate of the vertices of Bcp and Btag in good approximation, 

At = 4f (11) 

IPC 
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Method 


e 


UJ 


High momentum lepton 


0.014 ±0.021 


0.071 ±0.045 


Kaon 


0.279 ± 0.042 


0.199 ±0.070 


Medium momentum lepton 


0.029 ±0.015 


0.29 ±0.15 


Slow pion 


0.070 ±0.035 


0.34 ±0.15 



Table 6. The flavour tagging efficiency(e) and the wrong tag fraction (w) 



where 7/5(=0.425) is the boost factor of T45. The Bcp vertex is determined from the 
two lepton tracks in the J/ip decay. The Btag vertex is determined by the tracks used 
for the flavour tagging after poorly measured tracks are removed. The expected vertex 
resolutions are ~ 40 /im and ~ 85 /im for the Bcp and Btag vertices, respectively. 

The resolution of At, i?(At), is parametrized by two Gaussian distributions, where 
the first Gaussian is for the intrinsic vertex resolution and the effect of the secondary 
charmed mesons in the Btag side, and the second Gaussian accounts for the tail due to 
poorly measured tracks: 

R{^t) = — 7^exp( — )± 7^exp( ). (12) 



aV2n ' ^'^ ^ 2a^ ' ■ (TtaiiV2^ " " ' '^(^Li 
where ftaii = 1 — /- The means {fi and fitaii) and widths (cr and crtaii) of the two Gaussians 
are calculated event-by-event from the errors on the two vertices. The fraction of the first 
Gaussian, /, is 0.96 ± 0.04, determined from full MC simulation studies and B D*lv 
data. 

2.2.6. Extraction of sin20i The probability density function with a CP eigenvalue 

of Tjf is 

S{At,r]f,q) = — exp(-At/rBo)x{l-g(l-2cj)r7/sin20i sin(AmrfAt)}(13) 

Tbo 

where q = 1(— 1) if Btag = B^{B^). The wrong tag fraction, u, is a function of the 
tagging method as listed in table p. 2. 4 We fix tbo and Am^ to the world averages]^. 



1.548 ± 0.032 and 0.472 ± 0.017 ps ^, respectively. The probability density function 
for backgrounds is B{At) = -^exp{—At / r^kg) where the lifetime of the backgrounds, 
0.73 ± 0.12 ps, is obtained from the side bands of the signal in AE and Mheam- 
To extract sin20i , we define the likelihood of an event: 

/+00 /"+00 
RiAt - s)S{s, r]f, q)ds ± (1 - fs) / R{At - s)B{s)ds (14) 
-00 J —00 

where fa is the signal fraction. The extraction of sin20i is done by minimizing the 

log-likelihood, — ^ ■ In function of sin 20i . The results are summarized in table 



To verify our analysis, we analyze control data sample which should not have 
any CP asymmetry. The control data samples are B^ — >■ J/ipK*^{K*^ K~^7r~), 
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10 15 




Figure 8. (a) The fitted proper-time difference distribution for all candidates : 
Y.^j=-i,+i{dNld{r]fM)\BO + dN/d{-r]fM)\Bo). (b) The log-likelihoods for CP + 1, 
CP - 1,'and both CP±l. 



Decay mode sin20i 

All CP - 1 modes +0.81;! 

All CP + 1 modes -0.61^ 

All CP ± 1 modes +0.45;! 



Table 7. Rcsuhs of the CP fit 



— > J/ipK^ , B~ D^TT^ , and D*^l^v. All results obtained from the control 

samples are consistent with zero asymmetry. 

In table |2.2.6| , the systematic errors on sin 20i are listed. The largest error is due to 
uj, which is obtained from data. The total systematic error on sin 20i is the quadratic 
sum of all sources, -1-0.07 for positive side and —0.09 for negative side. 



2.2. 7. Conclusion We collected 6.2 fb ^ of data on the T45 by the end of the summer 
2000. Using this data sample, we made a preliminary measurement of sin20i using 
B ^ J/^Ks{iT+7r-), B J/ijKs{7r%'>), B ^ ^' Ks{7r+7r-), B ^ xci^5(vr+vr-), 
B J^vrO, and B J/^Kl. We found sin20i = +0.45+°:^^(stat.)+°:°^(syst.). Our 
result is consistent with the standard model prediction. We expect to improve the 
statistical errors as more data become available. 
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Source 


+a 


—a 


VViUiig, i-L city LlUiil I 




u.uuu 


Resolution function 


+0.026 


-0.025 


Background shape 


+0.029 


-0.042 


Background fraction 


+0.029 


-0.032 




+0.005 


-0.006 


IP profile 


+0.004 


-0.000 


Total 


+0.07 


-0.09 



Table 8. Summary of the systematic errors on sva2(j)i 
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2.3. CP Violation in i?-Decays at CDF: Results and Prospects. 

Farrukh Azfar, Oxford University 

2.3.1. Introduction Charge-Parity (CP) violation in particle decays is necessary 
to explain the matter-anti-matter asymmetry observed in the universe today. The 
amount of Standard Model (SM) CP violation is too small to account for the observed 
asymmetry. A detailed study of CP violation provides us an excellent opportunity to 
search for physics beyond the Standard Model. 

The CKM matrix CP violation in the Standard model has its origins in the complex 
couplings of the Cabibbo-Kobayashi-Maskawa quark mixing matrix. The interference 
of various decay amplitudes are expected to give rise to large CP violating effects in the 
B system. 

The CKM matrix can be expressed in terms of 3 real parameters (p, A, A) and one 
imaginary [rj) parameter in a representation known as the Wolfenstein parameterization. 
This matrix is unitary and several orthogonality relations between its rows and columns 
can be derived. We represent one of these KidKTfc + ^cdVcb + ^tdYtb = 0, as a triangle in 
the complex plane. 

At CDF we intend to measure one side and two angles of the unitarity triangle, by 
measuring the time dependent CP asymmetry in the modes —>■ 7T~^7r~ , 5° —>■ K~^K~ , 
B^ — > J/ipKg and the B^ — B^ mixing parameter Xg as illustrated in Fig. |^. A 
measurement of the time dependent CP asymmetry utilizing B'^ —>■ J/ ipKg decays with 
100 pb~^ of data collected during Run-I has already been done at CDF 

The time dependent CP asymmetry will also be measured in modes where the SM 
asymmetry is expected to be small e.g. B^ J/'4^<P and 5° DfDj. A measurement 
of the width difference between the weak eigenstates of the 5° is complimentary to 
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Figure 9. CDF plans for measuring the unitarity triangle (red). 



such measurements and would utilize the same two modes. Comparisons of the width 
difference to the mixing parameter are also sensitive to non-SM Physics. 

Mixing and CP violation in the neutral B meson system \n p — p collisions 
B and B mesons are created as strong interaction eigenstates which then "mix" into 
each other due to second order weak interactions represented by the box-diagram. The 
heavy and light weak interaction eigenstates are very nearly CP eigenstates each with 
its distinct mass and width, the quantities used to describe this two state system are: 
Mh, Ml, Fh, Pi, AM = Mh- Ml, AP = P^, - P^, P = th,l = and 

_ AMd^s 

'^'■'^ r n n ■ 

a ^ _ 

If B and B can decay to a CP eigenstate | / >, CP violation can occur if there 
is more than one amplitude contributing to the decay. If the complex CKM phases in 
both amplitudes are different then they will interfere causing an asymmetry in the rates 
B ^\ f > vs. 5— s-|/>. The interference is caused by second order weak-interactions 
that are represented by the box or penguin diagrams. 

Depending on whether we know (tag) the initial and final flavours of the B and if 
the decay is to a CP eigenstate or a flavour specific final state, we will get a particular 
time evolution, these are summarized below: 



TL,H 



1. The flavour of the B is unknown and the decay is to a CP eigenstate: P{t) 
Examples are B^ D^D~ , B^ tt+tt" (untagged). 

2. Flavour at birth is known, but not at decay, and the final state is a CP eigenstate: 

{ — ^) 

P{t) ^ ^^^{1 ± A^'p'sinAMt ± A^^cosAMt) Examples: Ba J /i^K^ 
A^^Jf = sin 2/5, where A^% = 0, and 5° ^ n+n- A^Jp" ^ and Agj^' =^ 0, (tagged). 

3. Flavour at birth and at decay is known (mixing): P(t) ^ - — ^(1 ± cos(AMt)) 
Examples 5° ^ Dfl^v{v) 5° ^ Dfn^, 5° ^ D^l^v{v). 

Experimentally the path length before decay of the B meson, ct is measured. 
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2.3.2. The CDF detector, run-I and the run-II upgrade. In this section, we 
briefly describe components of the CDF detector relevant to 5-physics. A partially 
instrumented detector with several upgraded components, was tested withp— p collisions 
at a centre-of-mass energy of 2 TeV in October 2000. We expect the upgrade described 
to be completed by March 2001 0. The detectors present in Run-I are described along 
with their Run-II successors. 

Tracking in the central region Tracking in the Central Region is provided by wire 
drift chambers. In Run-I the Central Tracking Chamber (CTC) was used, this had ~ 
6000 axial and stereo sense wires with a transverse momentum resolution of = 0.3 
%, covering the region — 1 < < 1. For Run-II the Central Outer Tracker (COT), with 
the same coverage and Pt resolution will be used, this is already installed, and has ~ 
30,000 axial and stereo sense wires. The increased number of wires provide better 
somewhat improving particle identification. 

Silicon microstrip detectors Silicon Microstrip Detectors for precise vertex 
determination are needed for B lifetime measurements. The Run-I Silicon Vertex 
detector (SVX), had four axial layers and a coverage of —1 < 77 < 1. The proper decay 
length resolution was cTct'- 35 /im = 0.12 ps. This detector provided measurements in 
the r — plane only. The Run-II Silicon Vertex Detector is known as the SVX-II, and 
has seven stereo and axial layers providing coverage in —2 < r] < 2 with an additional 
layer (Layer 00) at a distance of 1.4 cm from beam-pipe. Axial and Stereo strips provide 
the capability to do 3 dimensional tracking using only the SVX-II. The SVX-II acr is 
expected to be 15 /xm = 0.045 ps. 

Particle identification A Time Of Flight (TOF) detector built from scintillators 
has been mounted on the outer surface of the COT, this provides us with some particle 
identification capability at low momenta. Differentiating kaons from pions is crucial 
for tagging the flavour of Bs. In addition to the TOF, the coverage of the Muon 
Chambers has been extended from rj =1.0 in Run-I to rj = 1.4. High momentum muons 
are easily identifled, and allow us to select events with b — > fiu^X and b cc with 

CC ^ J /lp ^ 

2.3.3. Triggers and data selection At the Tevatron with p — p collisions at a 
centre-of-mass energy of l.STeV the b production cross-section is crj,^ =100yub, however 
the inelastic cross section is three orders of magnitude higher. Although cr^f; is much 
higher than for e"^ — e~ at the T(4S') and Z resonances (1 and 7 nb), the high inelastic 
cross section requires specialized triggers for the selection of bb events. 

CDF run-I triggers The selection of bb events during Run-I at CDF relied on high 
transverse momentum (Pt) leptons. The inclusive lepton i = (/x, e) trigger selected 
events with Pt{B) >8 GeV, the quark level decays being b — > iucX or b —>■ cX c — > iuX. 



CP Violation, UK Phenomenology Workshop 2000 



23 




Figure 10. A side view of the CDF detector. 



Dilepton triggers selected (/ze) and (/x/i) events with Pt{B) ^10 GeV were also in use. 
These were crucial in selecting b —>■ J/ipX with J/ip fifi and b —>■ fi^X, b e^X. 
The sin 2/5 analysis requires the reconstruction of the mode J/ipKg, for which 

the data sample was selected using the J/tp —>■ ^~ trigger. The long-lifetime of B 
mesons was not utilized in any trigger during Run-I. 

In Run-I the single lepton trigger data samples were used to measure B^ mixing. 
The presence of a neutrino in the final state introduces uncertainties in the B decay 
vertex determination. Despite this the CDF measurement of the B^ mixing parameter 
Xd = 0.77 ± 0.040 compares well with world average Xd= 0.739 ± 0.023 0. However 
this uncertainty will degrade the B^ mixing parameter Xg measurement significantly, 
since the oscillation period is much smaller. Clearly it is crucial to be able to trigger on 
fully reconstructible modes. 

CDF run-II triggers The Run-II CDF triggers used to select B decays will include 
all Run-I triggers, in addition a new high impact parameter track trigger using SVX-II 
hit information will be used in order to select displaced hadronic tracks from B decays. 
Therefore in Run-II CDF it will be possible to select B decays with fully reconstructible 
purely hadronic decay states, such as B^ DfDj, B^ tt+tt^ and B^ Dfn^. 
These decays allow us to measure the B^ width difference and mixing. 

2.3.4- The CP asymmetry, sin 2/3 in B^ J/ipK^: run-I To measure the CP 
asymmetry in B^ — > J/ipKg, the flavour of the B^ at production has to be tagged, and 
the decay fully reconstructed. The tagging of the flavour is neither fully efficient nor 
correct each time, we define the dilution variable D = where Nji and Nw are 

the number of correct and incorrect tags. If the tagging efficiency is given by e then a 
CP asymmetry measured using tagged events will have the statistical power of eD^N 
fully efficiently and correctly tagged events. 
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Flavour tagging: opposite and same side Opposite side tagging algorithms use 
decay products from the h quark on the opposite side of the reconstructed decay of 
interest. The decay products used are the jet associated with the hadronization of the 
h into a 5 or the sign of the lepton in case the B decays leptonically. In case of non- 
leptonic decays a weighted sum of charges of all tracks is used to tag the B flavour. The 
lepton tagging is known as the Soft Lepton Tag, has high dilution but low efficiency, 
the figure of merit, eD^ is 0.91 ± 0.1 %. The Jet Charge technique is more efficient but 
has worse dilution, with ^0"^= 0.78 ± 0.12 %. 

There are two disadvantages of the opposite side tagging techniques: the opposite 
B is within detector acceptance only 40 % of the time, and if the opposite B is neutral 
then it may have mixed. These factors degrade efficiency and dilution. 

The same side tagging (SST) algorithm uses the hadronization products on the 
same side as the reconstructed decay to determine the flavour of the B. This has a 
flgure of merit eD^= 1.8 ± 0.4%, as measured in Run-I. 



Results from run-I The candidate events for reconstructing the decay 5° J / ipKg 
with J/ip ^ and Kg tt+tt^ were selected using the Run-I J/ip ^ trigger. 

Roughly 400 such events were reconstructed with some 200 within the acceptance of 
the SVX. The tagging algorithms described were tested and tuned on a sample of 
B^ J/ipK^ decays, and the dilutions determined. The flavour of the B is already 
known in this decay mode, so it can be used for tuning the algorithms. 

All 400 events were used in an unbinned likelihood multi-parameter flt of the mass 
and decay length distributions, and tag information of the reconstructed B^ — > J/ipKg 
decays, each tagging algorithm was used for each event with an appropriate weighting 
accounting for agreement or disagreement between the taggers. 

The reconstructed mass spectrum and time dependent asymmetry for B'^ — > J/ ipKg 
are shown in Fig. |Tl|. The best flt value of the parameter sin 2/5 is 0.79 ± 0.39 (stat) ± 
0.16(sys) where the systematic error is dominated by uncertainties in the dilution. 



2.3.5. Prospects for run-II The plan for Run-ll is to constrain the unitarity 
triangle by measuring one side using B^ mixing, two angles 7 using CP asymmetries in 
B^ vr+vr" and B^ K^K~ and /3 (sin 2(3), as in Run-I, by measuring the asymmetry 
in B^ J/ipKg, this is shown in Fig. ^. 

Preparing the ground for projections We expect at least 2fb^^ of data in Run-II 
which is a factor of 20 higher than Run-I, this flgure has been used as the basis for 
all projections in this article, it is important to note that this is a very conservative 
estimate, the laboratory director's recent stated goal being 15 fb^^. 

In addition to the increase in luminosity, the increased coverage of the SVX-II and 
Muon chambers will increase our statistics for all Run-I modes by a factor of 50-60. 

All Run-I analyses such as the CP asymmetry in B^ — >■ J/ipKg will be repeated 
with much higher statistics and several new modes, selected using the hadronic displaced 
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track trigger will be used to search for CP violation and mixing. 

sin 2/9 in run-II We expect at least 10,000 events in the channel —>■ J/ifjKg, 
therefore the statistical error crsm2f^ will be ~ 0.067. We expect a similar decrease in the 
systematic error, since it is dominated by errors in measuring dilution, due to the larger 
B^ — > J/ipK^ calibration sample. We have not included the contribution of other 
decays that the hadronic displaced track trigger will select: such as B^ D*^D*~, 
B^ —>■ D^D~ and B'^ —>■ (pK^, where the time dependent CP asymmetry is also 
proportional to sin 2/3. 

7 using B^ — > TT+TT^ and B^ K^K~ In the absence of penguin-graphs the 
asymmetry in the decay B^ tt+tt^ will measure sin 2(7 + (3) and similarly the 
asymmetry in the decay 5° K^K^ should only measure sin 27. However penguins 
are present in both decays, the tt+tt^ mode is tree dominated the K^K~ mode is 
penguin dominated. The CP asymmetry for B^ —>■ tt+tt" is A^p^^^. cos AM^t + 
^cpV^- sin AMrft and for the mode 5° K+R- it is cos AM,t + 

^CPK+K- sin AMgt. The four Acp terms can be expressed in terms of 4 parameters, d 
the ratio of the penguin to tree amplitudes (assuming SU(3) symmetry), 9 the ratio of 
the phases of the penguin and tree amplitudes, and the CKM angles 7 and (3. These can 
be solved for by using the 4 asymmetries and the measurement of sin 2/3 0] . We expect 
the hadronic displaced track trigger to select 5000 B^ tt+tt^ and 10000 B^ K^K~ , 
and a simulation based on these estimates was performed. An estimated SU(3) breaking 
of 20 % was put into the simulation and the resulting inaccuracy is part of the combined 
systematic and statistical error of 10 degrees in measuring 7. All estimates of background 
were based on Run-I data. 
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Determination of Xs the mixing parameter The ratio of the mixing parameters 
1^ of the and i?° neutral mesons is proportional to |^^- Since Vts » Vtd, the -B° 
oscillates much faster than the B^, and we expect Xs ~ 25 {xd ~=0.75) 

The measurement of Xg oscillation is challenging. In Run-I we measured the 
oscillation parameter of the B^ meson, this means measuring an oscillation of period 
2.12 ps with a ct resolution of 0.1 ps. To measure the Xg we will try to measure a period 
of 0.07 ps with the SVX-II resolution of 0.045ps. 

We expect 20,000 5° Dfir^ and 5° Dfn^n^Ti^ with Df 07r± which the 
hadronic trigger will select. We have estimated our accuracy for measuring Xg with a 
signal to noise of 2:1 and 1:2 based on Run-I estimates. We expect to measure a Xg of 
up to 63 with a significance of 5cr. This accuracy in measuring Xg corresponds to an 
accuracy of 7 % in measuring | \ . In addition to measuring this side of the unitarity 
triangle the value of Xg can be used as a constraint in CP asymmetry measurements in 
CP violating decays. 

2.3.6. Non-SM surprises from the B'^ width difference and CP violation in 5° 
decays In the standard model the B^ width difference AP can be written in a CKM 
independent form in terms of the top, bottom, W and charm masses and the B^ mixing 
parameter Xg [Q. The fractional width splitting between the light and heavy 5^ ^ is 
expected to be large ~ 0.15. If we measure both Xg and (7r^)_B0 we can test the SM 
prediction. In particular, we are interested in cases where Xg is too large to measure 
but ^ is measurable or conversely, if ^ is too small to measure and Xg is measurable. 
In the presence of new Physics a non-SM CP violating phase (pBSM can appear, and 
will reduce the width difference: AT Measured = AT sm cos (pBSM- Complimentary to 
this effect the CP asymmetry in modes such as B^ J/'4"P be proportional to 
sm2(t)BSM- 

In Run-I approximately 60 5° J/ 4^4' decays were reconstructed from a data 
sample selected using the J/ip ^ trigger and a single lifetime was measured 0. 

The final state J / ip(f) is however to a mixture of CP even and odd final states, it therefore 
contains two lifetimes, for the heavy (~ CP odd) and light (~ CP even) 5° states. 
Using angular variables to disentangle the CP content of the final state we can fit for 
two lifetimes utilizing a likelihood function normalized over invariant mass, lifetime and 
the angular variable [0. We expect a signal of 4000 B^ J/4'4> Run-II. A simulation 
based on 4000 signal events and Run-I observations of signal and background shapes 
has been used to predict that we can measure a (^)50 of 0.15 with a precision of 0.05. 
This does not include any contribution from the purely CP even mode 5° D+D^ 
which will be selected using the hadronic displaced track trigger, the lifetime measured 
in this mode can be compared with the results of a single-lifetime fit to semi-leptonic 
5° decays and (■^)so can be extracted. 

The SM prediction for CP asymmetry in the B^ —>■ DfDj and B^ —>■ J / ip(f) modes 
is of order 3%, we do not expect to be able to see this at CDF, however if there are 
non-SM sources of CP violation, we may be able to observe this, by doing a conventional 
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tagging analysis in 5° DfD^ and a tagging analysis of the disentangled CP states 

of 5° ^ J/tlj(j). 

2.3.7. Conclusion During Run-II at CDF we have been able to measure of 
sin 2/? =0.79 ± 0.39 (stat) ± 0.16(sys) and the mixing parameter Xd = 0.77 ± 
0.040(stat) ±0.039 (syst), and have fully reconstructed the largest sample of 5° decays. 

Thus the viability of a rich CP violation program for Run-II, during which we expect 
a factor of 50 increase in data, has been established. The introduction of the displaced 
track trigger will allow us to select fully reconstructible hadronic modes of all B mesons. 
We expect to be able to constrain one side and two angles of the unitarity triangle. We 
expect to measure | | to ~ 7% accuracy, using mixing. We emphasize that 
for the next few years CDF is unique in its ability to analyse B^ decays. We expect 
to use — >■ 7r+7r~, B^ —>■ K^K~ to measure 7 with a precision ~ 10", and using 
B^ — > J/ipKg to measure sin 2/5 to a precision of 0.067. 
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2.4. B physics potential of future experiments at hadron machines 

V Gibson, Cambridge 

2.4.1. Introduction CP violation remains one of the enigmas of particle physics 
today. Experimentalists have just started a full programme of research to study CP 
violation in B decays. The experiments can be divided into two phases. The main goal 
of the current phase of experiments will be to observe CP violation in the B system 
in the exciting exploratory phase, with the potential to establish a breakdown in the 
Standard Model [|I|]. This review concentrates on the next generation of experiments 
which become operational around 2006. The experiments considered are the dedicated 
B physics experiments, LHCb at the Large Hadron Collider (LHC) and the proposed 
BTeV experiment at the Tevatron, as well as the general purpose experiments 
ATLAS and CMS at the LHC. The main aims of these experiments will be 
to make precision measurements of CP violating observables using many different decay 
channels and species of B hadrons. They will thoroughly test the internal consistency 
of the Standard Model description of CP violation and have the sensitivity to search for 
the necessary new physics beyond. 
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CP violation in the Standard Model and beyond CP violation arises naturally 
in the Standard Model through the presence of a single phase in the unitary Cabibbo- 
Kobayashi-Maskawa (CKM) quark- mixing matrix [Q. The unitarity of the CKM 
matrix is clearly exposed when using an explicit parameterization. A very popular 
parameterization is the perturbative form suggested by Wolfenstein [|^, which can be 



expanded to order (A = sin 6, 



( 



V, 



CKM 



c ~ 0.22 where 6c is the Cabibbo angle): 
A A\^{p - iri) \ 



X i-\x^ 

- p - it]) -AX^ 



AX^ 
1 



+ 



o\ 








-iA^X^T] 
y AX^{p + ir])/2 AX^{l/2 - p - ir]) 

The parameter rj represents the CP violating phase in the Standard Model and appears 
in three of the matrix elements. The unitarity of the CKM matrix implies that there 
are six orthogonality conditions, which can be represented geometrically as triangles in 
the complex plane. Two such unitarity triangles, shown in Figure are expected to 
have angles, which are all non-trivial. 

The angles of the unitarity triangles are all related to the single CP-violating phase 
in the CKM matrix and are designated by a, /?, 7 and ^7; 



(3 = tan 



7 = tan^ 



-1^ 
P 



S'j = r]X^ 



l-p(l-AV2). 

and a = TT — /5 — 7. The angles f3, ^7 and 7 are commonly referred to as the Bj mixing 
phase, the 3° mixing phase and the weak decay phase respectively. By 2006, it is expected 
that a measurement of sin 2/3 will have been made with a precision of ~ 0.02 There 
will be no good or direct measurement of 7 and there will be no sensitivity to ^7. 




Figure 12. The unitarity triangles cliaracterized by the relations VudVub + ^cdV*fj + 
ytdVfl = and KdVf^ + VusVfl + VubV^l = 0. The phase convention is chosen such 
that VcdV*i^ is real and all sides are normalized to |VcdVj^j|. 
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It is expected that the single Standard Model CKM phase is insufficient to explain the 
observed baryon asymmetry in the universe |^ and that new physics must intervene. A 
large source of CP violation at the electroweak scale could be provided for in extensions 
to the Standard Model |Ty] . A summary of physics beyond the Standard Model and its 
effects on the CKM parameters can be found in references ITT] , |12[. In order to search 
for new physics it is essential to measure and calculate as many processes as possible 
and to compare the resulting CKM parameters with each other. 



The next generation experiments The LHC and Tevatron colliders will provide an 
intense source of the full spectrum of B hadrons (B^, B^J, B^, B^ and b-baryons). The 
running parameters of the colliders and experiments are given in Table |^. The expected 
cross-section for the production of bb pairs at the LHC is approximately five times that 
expected at the Tevatron, with a relatively smaller inelastic cross-section. 





Tevatron 


LHC 


Energy /collision mode 


2.0 TeV pp 


14.0 TeV pp 


bb cross-section 


~ 100/ib 


~ SOO^ib 


Inelastic cross-section 


~ 50mb 


~ 80mb 


Ratio bb/inelastic 


0.2% 


0.6% 




Bunch spacing 


132 ns 


25 ns 




Bunch length 


^ 30cm 


~ 5cm 




BTeV 


LHCb 


ATLAS/ CMS 


Detector configuration 


Two- arm forward 


Single-arm forward 


Central detector 


Running luminosity 


1.5 X 10^2 cm-^s-i 


1.5 X lO^^cm-^g-i 


1 X 10^'^ cm~^s^^ 


bb events per 10^ sec 


2 X 10^^ X accep 


1 X 10^^ X accep 


5 X 10-*^^ X accep 


(Interactions / crossing) 


~ 2.0 


~ 0.5 (30% single int.) 


~ 2.3 


Mass resolution B|| — > J/^/^Kg 


9.3 MeV/c2 


7 MeV/c^ 


18/16 MeV/c^ 


Proper time res. Bg ^t^^ 


43 fs 


43 fs 


50/65 fs 



Table 9. Comparison of the LHC and Tevatron collider running parameters and 
experiments. 



The forward detector geometries of the LHCb and BTeV experiments exploit the 
expected forward peaked and strongly correlated production of b and b hadrons. LHCb 
is instrumented on one arm with a dipole spectrometer and is designed to run at 
a low LHC luminosity of 1.5 x 10^^cm~^s^^. It employs a precision ministrip silicon 
detector and has an efficient multi-level trigger, which includes a vertex, trigger at the 
second level. The experiment employs two RICH detectors for particle identification 
and has hadronic and electromagnetic calorimetry. The key design features of the 
BTeV detector include a forward two-arm spectrometer, a precision silicon pixel vertex 
detector, a vertex trigger at the ffist level, a RICH detector and a lead-tungstate 
calorimeter for neutral particle reconstruction [§. The geometrical acceptance, the use 
of a vertex trigger at the ffist level and of multi-bunch interactions mostly compensate 
for the smaller bb cross-section at the Tevatron. ATLAS and CMS are central detectors 
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designed for general-purpose use at the LHC ^. During the first three years of 
running, the LHC will operate at a low luminosity, 1 x 10^^ cm~^s~^, thereby enabling 
ATLAS and CMS to pursue a B physics programme. ATLAS and CMS have tracking 
capabilities in the central region and employ precision silicon pixel and microstrip 
detectors. Specialist B triggers are achieved by reducing the lepton pt thresholds to 
a minimum. 

The requirements of the dedicated experiments are governed by the need to measure 
time-dependent CP asymmetries for B and B hadrons decaying to the same final state. 
The experimental acceptance and trigger efficiencies mainly cancel. However, precise 
measurements require good decay time and mass resolutions, efficient triggers for low 
and high multiplicity B final states, particle identification for vr/K separation and photon 
detection for neutral final states. An example of the B|^ J/ipY^^ mass resolution is 
given in Table The flavour of the B hadron at production needs to be determined 
either using the signal B or the other B in the event and good control of systematic 



uncertainties is crucial [1^. The experimental need for good flavour tagging and B 
proper time resolution is demonstrated through the measurement of B^ oscillations. 
The bench-mark decay mode for this measurement is B^ — > D^tt^ in which the flavour 
of the Bg is given by the charge of the Dg. The proper time resolutions obtained by 
the experiments is summarised in Table |^. The event yields and expected measurable 
values for the B^ mixing parameter are given in Table |T0|. 

2.4-2. Direct measurements of CKM angles and search for new physics 

The internal consistency of the Standard Model description of quark mixing in weak 
interactions can be thoroughly tested by measuring CP violating observables in the 
decays of B mesons. If the Standard Model is correct then all such measurements will 
be describable with a single set of CKM parameters. New physics outside the Standard 
Model could lead to additional phases in the CKM matrix and an inconsistency between 
the measurements. This review discusses the current status of studies performed by the 
new generation experiments to extract the CKM angles a, /?, 7 and ^7. The event yields 
and sensitivities presented are mainly taken from references and ||12| where further 



details can be found. Since the experiments are at different stages in their preparation 
for physics, the results quoted should only be considered as a current snapshot. 

B^ mixing phase The decay B^ J/ip'K^ is a transition into a CP eigenstate and 
is dominated by only one CKM amplitude. Hence, the time-dependent CP asymmetry 
is governed by the B° mixing phase, (3. The decay is also experimentally clean and 
can be reconstructed with relatively low backgrounds. Examples of the B^ — > J / ip'K^ 
reconstructed mass distribution and CP asymmetry are shown in Figure 0. The 
expected event yields and sensitivities for the next generation experiments are given 
in Table Ultimately, a combined statistical precision of 0.005 should be achievable 
at the LHC. 
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Figure 13. The ATLAS B[J -> J/f/iKg, J/?/; e+e" reconstructed mass distribution 
after 3 years of data and the LHCb CP asymmetry with one year's data. 



The Bg mixing phase The decays 3° — ^ J / ip(p and 3° — J/ z/^t^^'^ are 3° counterparts 
to the decay mode B[| ^ J/'j/'Kg. Once again these decays are dominated by only 
one CKM amplitude and are sensitive to the 3° mixing phase, ^7. Experimentally, 
the decay mode 3° — > J/ipcf) is very clean, but requires a full angular analysis to 
disentangle the mixture of CP eigenstates in the final state. The reconstruction of 
the decay 3° J/ipr]^'^ benefits from good electromagnetic calorimetry and, because 
it is a transition into a CP eigenstate, the ^7 can be extracted directly from the decay 
asymmetry. The expected sensitivities of the experiments for the extraction of ^7 
are summarised in Table |10| and depend strongly on the proper time resolution and 
3g mixing parameter. It is interesting to note that the expected sensitivity for the 
extraction of from 3^ J / ipr]^'^ decays in one year of running (3TeV) is comparable 
to 5 years of running for the full angular analysis of 3° J/'^4> (LHCb). It has been 
suggested that some of this difference can be recuperated using a transversity angle 
analysis |13| . A particularly interesting feature of these decays is that they exhibit very 
small CP violating effects within the Standard Model, = X^rj ^ 0.03, and hence are 
very sensitive to new physics (see []T2[ and references therein). 



The angle a 

3j p^Ti^: A full three-body analysis of the decay 3^ tt+tt^tt'^ in the 
p resonance region, taking into account interference effects between vector mesons 
of different charges, has been proposed to extract all parameters that describe both 
the tree and penguin contributions to 3° — >• p^tt"^, including the CKM angle a |n 



Although the method is theoretically clean, it is experimentally difficult due to the need 
to reconstruct vr^'s and the presence of a large combinatorial background. The invariant 
mass for 3^ p^n^ candidates reconstructed in 3TeV and the tt+tt^tt*^ Dalitz plot from 
LHCb are shown in Figure H. The event yields, given in Table are expected to be 
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Figure 14. The tt+tt 7r° invariant mass distribution for background and signal 
events for p°7r° decays from BTeV and the LHCb DaUtz plot. 



Bj — ^ TT+TT^: In principle, the decay mode tt^ti^ allows the CKM angle a 

to be probed. However, penguin contributions to the transition amplitude introduce a 
direct CP violation contribution, Aq^, into the decay asymmetry, 

A{t) = cos(AMt) + sin(AMt). 

Experimentally, the analysis of B^ — ^ tt+tt^ decays is complicated by the existence 
of backgrounds with similar topologies and with unknown CP asymmetries, such as 
BO ^ TT+K-, BO ^ K+K-, BO ^ tt+R-, Ab ^ pvT" and Ab pK". This is 
illustrated in Figure |l5| (left) which shows the reconstructed two pion invariant mass 
distribution obtained in ATLAS. In order to reject the non tt+tt^ background, LHCb 
exploits the powerful RICH particle identification detectors, shown in Figure |15| (centre) . 
The expected vr+vr" event yields and the experimental sensitivities for the CP 

observables are given in Table |10|. The observables can be written to leading order in 
the ratio of penguin to tree amplitudes, |P/T|, 



= - sin 2a -2 



sin a cos 2a cos(5) , A Jp 



sin a sin (5), 



where 6 is the CP conserving strong phase, 6 = arg{PT*). Unfortunately, a theoretically 
reliable prediction for |P/T| which would allow the extraction of a, albeit with a four-fold 
discrete ambiguity, is very challenging. Figure ^ (right) shows the expected sensitivity 
to a for different values of the uncertainty on |P/T| and integrated luminosities. It can 
be seen that for values of a around 90°, the sensitivity to a is already limited after one 
year of running, if the uncertainty on |P/T| is not better than 10%. 



The weak decay phase 7 
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Figure 15. The ATLAS tt+tt" invariant mass distribution (left), the LHCb tt+tt" 
invariant mass distribution after the application of RICH information (centre) and 
the expected LHC combined sensitivity to a as a function of a (right), for 5 = 30°, 
|P/T| = 0.2 after one year (dashed lines) and 5 years (solid lines). 



Bj — > TT+TT and Bg K"^K : A strategy has been proposed to combine the 
CP observables from the decay B^ tt+vt^ with those from B^ — K^K~ to provide 



a simultaneous determination of the angles 13 and 7 |15[. The decays B^ — tt^^tt and 



Bg — i> K+K~ are related to each other by interchanging all down and strange quarks 
(U-spin symmetry). Assuming the B^ mixing phase ^7 is known, the four observables, 
depend on four unknowns: two hadronic parameters, the B° mixing phase (3 and the 
weak phase 7. If (3 is also fixed from external measurements then the weak phase, 7, 
can be extracted in a theoretically clean way. This approach is unaffected by penguin 
topologies and final state interaction effects and is only limited theoretically by U- 
spin breaking effects. Moreover, since penguin processes play an important role in the 
decays of B^ — > 7r"*"7r~ and B^ — >■ K+K~, this strategy is very promising with regard to 
the search for new physics. Experimentally, the experiments expect large statistics for 
both B° — i> TT+vr^ and B^ K+K" decays. The expected event yields and sensitivity 
for 7 are given in Table |10|, where the range of sensitivity quoted reflects the dependence 
of the result on the B^ mixing frequency. 



B± ^ D^K^ and B^ ^ D°K*°: The decays B± ^ D°K± and B° ^ D°K*° 
are pure tree decays and can be used to give a theoretically clean determination of the 
weak decay phase 7. The approach is through the measurement of six exclusive decay 
rates, either B+ DiK+, B+ ^ DOK+ and B+ ^ D°K+ + c.c; or B^ ^ DOK*°, 
B° ^ D°K*° and B^ DiK*° + c.c. Here, Di = (D° + D0)/V2 is the CP even state 
and leads to relationships between the decay amplitudes that can be used to extract 
27 with a two-fold ambiguity. Experimentally, the separation between the decay modes 
B+ D°K+ and B+ ^ WK^ is extremely difficult as the D° and decay to the 
same hadronic final state. Also the decays of neutral D mesons into CP eigenstates, 
such as Di TT+TT^, K+K", require excellent particle identification. A method has been 
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suggested to make use of the large interference effects by measuring at least two different 
final states of the D° and |]T6[ . BTeV have investigated this approach and have studied 



the two decays modes; D° K+vr" and D° — K"'"K~. The annual event yield for this 
type of decay assuming branching ratios of 10^^—10^® is given in Table [1^. The expected 
sensitivity for 7, quoted in Table |10], depends strongly on the strong phase and the value 
of 7 itself. LHCb, benefiting from the hadron trigger and RICH particle identification, 
have investigated the possibility to determine 7 in the family of B° D'^K**' decays 
and have demonstrated that it will be possible to reconstruct samples of such events. 
However, the visible branching ratios are very small, 10^^ — 10^'', leading to low event 
yields. 



B° D(*)^7r'F: The decays Bj D(*)=^7r'F are transitions into non-CP 
eigenstates and receive contributions only from tree diagrams which lead to interference 
effects between the B° — B° mixing and decay processes. Measurements of the time- 
dependent asymmetries for the final states D"??"*" and D"'"7r~ lead to a measurement 
of 2(3 + 7. LHCb have investigated the potential of measuring 7 using B^ — > D'^*^^7r'^ 
decays where the D^*^^ decays strongly. Two methods have been studied: a conventional 
exclusive reconstruction with K^n" and an inclusive approach. For the inclusive 

approach, the momentum of the B is reconstructed using the momenta of the pion 
coming directly from the B, the momenta of the pion from the D^*-*^ decay and the 
direction of the B. The expected error on 2/5 + 7 depends strongly on the value of 
2/? + 7 and is > 11° for one year's data. 

Bg —>■ D^K^: The decays Bg —>■ D^K^ are the Bg counterparts to the 
BJI — > D^*'''''7r^ decay modes and likewise receive only tree diagram contributions, thereby 
probing the CKM angle combination —2^7 + 7 in a theoretically clean manner. The 
interference effects in B^ D^K^ are much larger as they are not doubly-Cabibbo 
suppressed as in the case of B^J D(*)^7r^. Experimentally, the event selection is 
challenging as B^ —>■ D^vr^ events, which come with a 20 times larger branching ratio, 
need to be rejected. The LHCb and BTeV RICH detectors are therefore crucial to the 
analysis of this decay. The precision with which the CKM phase combination can be 
measured after one year's operation, given in Table |TU|, depends strongly on the B^ decay 
width difference and mixing frequency. 



B ttK: Due to the dominant role of QCD penguins, B ttK decays 
potentially offer a determination of the weak decay phase 7 which is sensitive to new 
physics |T^. Experimentally, the strategy involving K+vr" and K°7r+ final states provides 
the cleanest channel. However, the measurements require the knowledge of the trigger 
and reconstruction efficiencies for the different final states and hence will incur an 
additional source of systematic error in contrast to most CP asymmetry measurements. 
A value for 7, with a four-fold ambiguity, can only be extracted once the ratio of tree 
to penguin decay amphtudes is determined from theory. This is limited by rescattering 
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effects and colour suppressed electroweak penguins. A preliminary study by LHCb 
shows that a precision ~ 2° and ~ 7° for two of the solutions can be obtained if the 
ratio of tree to penguin amplitudes is 0.18 ± 10%. 



Bd(s) J/'j/'Kg, 'Qd{s) D^^D^^: A strategy to extract the weak decay phase 
from B^(s) — i> J/'j/'Kg and 'Qd{s) Dj^D^^ decays has recently been proposed |[T8[ . 
The method makes use of the U-spin symmetry of the decays and is sensitive to new 
physics due to the presence of penguins. For example, the decay 3° — > J/ipK.^ is used 
to extract three observables from the time-dependent CP asymmetry. The ratio of 
the Bg — > J/^/^Kg to B° — > J/ipK^ untagged time-integrated decay rates provides a 
fourth observable. These are then used to extract the four unknowns : two hadronic 
parameters, 7 and ^7. Experimentally, although the J/ip'K.^ final state is clean, the 
isolation of — * J/^/'Kg events is challenging due to a low event yield, a large 
combinatorial background and the close B° mass peak. A study by CMS using an 
event selection tailored to B^ J/^/'Kg decays indicates that a measurement of the 
CP asymmetry in B^ J/^/'Kg is feasible at the LHC and that 7 can be measured 
with a precision of ~ 9° in 3 years operation. The strategy for extracting 7 from 
B(i(s) D(s)D(s) decays is analogous to the strategy for 'Qd(s) ^ Jl'4^'^s decays. The 
time-dependent asymmetries are measured using B° D+D^ decays and the overall 
Bj] D"'"D~ normalization is fixed through the CP-averaged B^ D+D^T decay rate. 
The method benefits experimentally as it is neither necessary to resolve the rapid Bg 
oscillations, nor does flavour tagging reduce the already suppressed yield in B^ events. 
However, as the final state consists of six hadrons, a hadron trigger and RICH particle 
identification are vital for this analysis. A very preliminary study by LHCb suggests 
that a precision on 7 of a few degrees should be achievable. 



Rare B decays Flavour-changing neutral current decays involving 6 — s or 6 — > d 
transitions occur only at a loop-level in the Standard Model with very small branching 
ratios and therefore provide an excellent probe of indirect effects of new physics. The 
study of rare B decays at the next generation experiments is at a preliminary stage and 
only a few remarks will be made here. A first assessment of the physics potential of 
these experiments shows that it will be possible to: 

• observe Bg — /i^/i" and measure its branching ratio ((9(10^^) in the SM), 

• perform a high sensitivity search for B^J fJ'~^fJ'~ (branching ratio O{10~^^)), 

• measure the branching ratio and decay characteristics of Bj^ K*'^7, 

• measure the branching ratios of B[} —>■ K*°/i"'"/i~, B^J — pji^ jj,' and B^ (pji^ jj,'' 
and study their decay kinematics; and 

• measure the forward-backward asymmetry in Bj^ K*'^/x+/i^ decays, allowing the 
distinction between the Standard Model and a large class of SUSY models. 

The expected number of signal and background events at the LHC is given in Table 0. 
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Measurement 


Channel 




Event Yield 






Sensitivity 










BTeV 


LHCb 


ATLAS 


CMS 


BTeV 


LHCb 


ATLAS 


CMS 


Mixing Phase 


sin 2/3 


Br 




80.5k 


88k 


165k 


433k 


0.025 


0.021 


0.017 


0.015 


B^ Mixing Phase 


sin 2(^7 


B°- 


J/v</'(t) 




370k 


300k 


600k 




0.03 


0.05 


0.03 


sin 2^7 


B°- 




9.2k 








0.033 








CKM Angle a 


m 


Br 




23.7k 


12.3k 


2.3k 


0.9k 


0.024 








A mix 
^CP 


Br 














0.07 


0.21 


0.14 


Adir 
^CP 


Br 














0.09 


0.16 


0.11 


sin 2a 


Br 




10.8k 


3.3k 








2.5° 






Weak Decay Phase 


7 




-^D"K± 


2.1k 








2-10° 








7 


Br 






0.4k 








~ 10° 






7 


Br 


K+K- 




4.6k 


0.54k 


0.96k 






4-8° 




2/3 + 7 


(Br 

B^- 


TT+TT") 

D(*)±7r=F 




340k 








> 11° 






-2^7 + 7 


Br 




13.1k 


6k 






6-15° 


2-12° 






Bg mixing 


.Ts(5fT) 


Br 


Df ttT 


103k 


86k 


3.5k 


4..5k 


75 


75 


46 


42 


R<ir(> B Dccavri (f) 


S/B 


Br 






33/10 


27/93 


21/3 










S/B 


BS- 






22.4k/1.4k 


2k/0.3k 













Table 10. A comparison of the expected physics performance for one year of running 
for all experiments proposing to measure CP violating observables in the B system. 
The numbers quoted for channels marked with a f are for 5 years of data for LHCb 
and 3 years of data for ATLAS and CMS. 



2.4-3. Concluding remarks The experimental search for the origin of CP violation 
has just entered an exciting new era. If the Standard Model description of CP violation 
is correct, then it is expected that CP violation in B decays will be discovered around 
the year 2001. However, in order to thoroughly test the internal consistency of the 
Standard Model and search for the necessary physics beyond, CP violating observables 
in the B system must be ultimately measured with very high statistics and in many 
different decay channels. The next generation of experiments will become operational 
around 2006 and will be the ultimate source of CP violation studies in the B system for 
many years to come. 
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3. CKM phenomenology and new physics 
and what can we learn from sin2f37 

3.1. Determination of the CP violating weak phase 7 

I 

C.S. Kim, Yonsei Univ. 



3.1.1. Introduction The source for CP violation in the Standard Model (SM) with 
three generations is a phase in the Cabibbo-Kobayashi-Maskawa (CKM) matrix. One of 
the goals of B factories is to test the SM through measurements of the unitarity triangle 
of the CKM matrix. An important way of verifying the CKM paradigm is to measure 
the three angles ||l[], 

= Arg[-(V,dV:J/(V,dV:b)], 
/? =Arg[-(V,dV:j/(VtdV:b)], 

and 7 = Arg[-(V,dV:b)/(VedV:j], (15) 

of the unitarity triangle independently of many experimental observables and to check 
whether the sum of these three angles is equalf]] to 180°, as it should be in the paradigm. 
It is well known that among the three angles, 7 would be the most difficult to determine 
in experiment. There have been a lot of works to propose methods measuring 7 using B 
decays, but at present there is no gold-plated way to determine this angle. In particular, 
a class of methods using B DK decays have been proposed [0, ^, || . 

§ The work was supported in part by Seo-Am (SBS) Foundation, in part by BK21 Program, SRC 
Program and Grant No. 2000-1-11100-003-1 of the KOSEF, and in part by the KRF Grants, Project 
No. 2000-015-DP0077. 

II The sum of those three angles, defined as the intersections of three lines, would be always equal to 
180", even though the three lines may not be closed to make a triangle. 
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3.1.2. Methods to extract 7 In Ref. 0, Gronau, London, and Wyler (GLW) 
suggested a method for extracting 7 from measurements of the branching ratios of 
decays D^K^, D^K^ and B^ DcpK^, where Dcp is a CP eigenstate. 

However, the GLW method suffers from serious experimental difficuhies, mainly because 
the process B~ —>■ D^K~ (and its CP conjugate process B~^ —* D^K^) is difficult to 
measure in experiment. That is, the rate for the CKM- and color-suppressed process 
B^ D^K~ is suppressed by about two orders of magnitudes relative to that for the 
CKM- and color-allowed process B^ D^K~, and it causes experimental difficulties 
in identifying through K^tx~ since doubly Cabibbo-suppressed K^n^ 

following B~ D^K~ strongly interferes with —>■ K^tt' following the rare process 
B- D^K-. 

To overcome these difficulties in the GLW method, a few variant methods have 
been proposed. The Atwood-Dunietz-Soni (ADS) method uses the processes B~ — > 
D^{D^)K^ fK~ with two neutral D decaying into final states / that are not CP 
eigenstates, such as / = K^ir'., Knir, etc. In this method, large CP asymmetries 
are possible since magnitudes of the two interfering amplitudes are comparable; i.e., 
the process B^ D^K^ fK^ is CKM- and color-suppressed in B decay, while the 
process B^ D^K^ fK^ is doubly Cabibbo-suppressed in D decay. The extraction 
of 7 can be allowed without measurement of the branching ratio for B^ D^K~ . Note 
that the decay amplitudes of B~ — > D^K~ and D*^ f contain the CKM factors V*gVcb 
and V*^Vus, respectively, while the amplitudes of B^ D^K~ and ^ f contain the 
CKM factors l/^Kb = \V*,Vub\e-''^ and 1/*^^. = iKTd^c.!, respectively. We define the 
following quantities: {i = 1, 2) 



a 


= A{B- - 


D'^K') = \A{B- 




b 


= A{B- - 


D^K-) = \A{B- 




Ci 


= A(D° - 


. U) = \A{D' ^ 




1 


= A(D° - 


. U) = \AiD' - /.) 


\e 


di 


= A{B- - 







where A denotes the relevant decay amplitude and (5's are the relevant strong rescattering 
phases. Similarly, we also define a, b, Ci, c[ and di as the CP-conjugate decay amplitudes 
corresponding to and di, respectively, such as di = A{B~^ [fi]K^), etc. 

Here [fi] in di denotes that fi originates from a or decay. Note that |x| = |x| 
with but in general \di\ 7^ \di\, as shown below. Then, the amplitude di 

can be written as 

di = |aQ|e(^"+^=») + Mle-'^e^'"^'^'^, (17) 
which leads to 

Idjp = |acjp + l&c^P + 2\abcic[\ cos(7 + Aj), 

|J.|2 = |aci|2 + \bc'f + 2\abcic'i\ cos(7 - A^), (18) 

where Aj = Sa — Sb + —^c'.- We see that \di\ 7^ \di\, unless Aj = nn {n = 0, 1, ...). Now 
the four equations for i = 1,2 in Eq. ([T8| ) contain the four unknowns 7, Ai, A2, 
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assuming that the quantities |a|, |cj|, |c^|, \di\ and \di\ are known, but \b\ is unknown. 
By solving the equations one can determine 7, as well as the other unknowns such as 
\h\ = \A{B- ^ D'>K-)\. 

In Ref . m , Gronau suggested a method to determine 7 using only the color-allowed 
processes, B~ D^K~ and B~ — > DcpK~ , and their CP-conjugate processes. 

In Ref. two groups, Gronau and Rosner (GR), Jang and Ko (JK), proposed 
a method to extract 7 by exploiting Cabibbo-allowed decays B D^*^K and using 
the isospin relations. In GR/JK method the decay modes B DK with the 
quark process h ucs contain the CKM factor |KifeV^s|e~''^ and their amplitudes can 
be written as 

A{B^^D^K°) =Aie'^'e-'\ (19) 

where Ai and 5i denote the amplitude and the strong re-scattering phase for the isospin i 
state. In this method, three triangles are drawn to extract 27, using the isospin relation 

A{B- D^K~) + Ai.B- D-R^) = A{B^ D^K^) (20) 

and the following relations 

A{B- DiR-) = A{B° DiK°) + -^A{B° D+K^), 

v2 

A{B+ DiK^) = A{B^ DiK^) + -^A{B^ D-K+), (21) 

v2 

where Di is a CP eigenstate of D meson, defined by Di = :^(-D° + D^). 

Recently a new method has been presented by Kim and Oh (KO) 0, which is 
similar to the ADS method, but uses B — > Dn decays instead of -B DK decays used 
in the ADS method. In fact, CLEO Collaboration have observed [0] that the branching 
ratio for B^ D^it~ is much larger than that for B^ —>■ D^K^ , 

B(B- D^K-) , , 

^ ^0^-) = 0-055 ± 0.014 ± 0.005 . (22) 

This new KO (Kim and Oh) method considers the decay processes B^ D^n^ /^r^, 
B^ D^Ti^ fiT^ and their CP-conjugate processes, where and decay into 
common final states / = K~^tt~, K'^p", Ktitt, and so forth. The mode B^ D^n^ 
is much suppressed relative to the mode B^ D^n^ , and this fact causes serious 
experimental difficulties in using B" D^n^ decays for the GLW-type method. 
However, in KO method one needs not to perform the difficult task of measuring the 
branching ratio for B~ —>■ D^vr", similar to the case of the ADS method. Note that 
the decay amplitudes of B^ — > D°7r~ and ^ f contain the CKM factors V*^Vcb 
and V*^Vus, respectively, while the amplitudes of B^ D^7t~ and ^ f contain the 
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CKM factors V^^V^b = {V^^V^ble-'^ and V^^V,, = |V;*^K.|, respectively: (i = 1,2) 
a = A{B- D\-) = \A{B- D\-)\e'^\ 
b = A{B- D\-) = \A{B- L)°7r-)|e-^V^', 

d, = A{B-^D%^f,)n-). (23) 
Then, the amphtude di can be written as 

d, = A{B~ ^ D°7r-)A(D° ^ /,) + AiB- ^ D\-)A{D'' ^ /,) 



= aci + be'- 

= |aQ|e('^"+'^^-) + Ibc'.le-'^e'^'"^'^'^. (24) 
Thus, and are given by 

= lacil"^ + |6c-|^ + 2\abcic[\ cos(7 + Aj), 

= |aQ|^ + l&c^l^ + 2|a6cjCj| cos(7 — Aj). (25) 

The expressions in Eq. (^) represent four equations for i = 1,2. Now let us assume 
that the quantities |a|, |cj|, \c[\, \di\ and \di\ are measured by experiment, but |6| is 
unknown. Then there are the four unknowns 7, Ai, A2 in the above four equations. 
By solving the equations one can determine 7, as well as the other unknowns such as 
\b\ = \A{B- ^ D''7r-)\. 

3.1.3. Experimental considerations Now we study the experimental feasibility of 
the ADS and KO method, by solving Eqs. (4,11) analytically. 



cosii ± A,) = r\-n ' " ■ (26) 



\d^\ 


2 _ 


aci\ 


2 _ 






2 


2 


\acibc[\ 





To make a rough numerical estimate of possible statistical error on 7, we use the following 
experimental result: 

BR{B- D\-) > BR{ D^R-) > BR{ D^^R-) > BR{ D\- 
oc C(IOO) : 0(10) : C(l) : O(O.l). 

Therefore, if we assume the precision of 1 % level experimental determination for the 
product of branching ratios, BR{B^ D^n^) x BR{D^ —>■ R'n^), then we roughly 
get 

|acip(7r) = 100 ± 1, \bc'-\'^{7i) = 0.1 ± 0.03 
\aCi\'^{R) = 10 ± 0.3, \bc'f{R) = 1 ± 0.1. 

Then, we can estimate the statistical error as 

A[cos(7± A)(5^ ^ D(^ /)7r±)] -0.3, 

A[cos(7 ± A){B^ f)K^)] ~ 0.15. 

Even though the ADS method can give approximately twice better precision statistically 
for determination of 7 than KO method, however, KO method can have other 
advantages: 
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• The value of \di\^ oc BR{B~ [/j]7r^) is an order of magnitude bigger than 
\di\^ oc BR{B^ —>■ [fi]K^). Therefore, if B-factories can produce only handful of 
such events, KO method can be only the possible option. 

• Systematic errors could be much smaller for KO method due to final state particle 
identification. 

We can extend the GLW method to B^. decay |Q from the relations, 

V2AiBt ^ DtD',) = A{Bt ^ DtD') + A{Bt - DtD^), 

V2A{B- ^ D^Dl) = A{B- ^ + A{B: ^ (27) 

and we can obtain 7 from A{B'^ DfD'^) = e^^^ A{B^ DjD^). Here the advantage 
is that the amplitude with the rather small CKM element V^b is not color suppressed, 
while the larger element Vcb comes with a color suppression factor. Therefore, the two 
amplitudes are similar in size, as 

A{B+ DjD^) : A{B+ DfD^) oc : X^N^, 

because of [Vufc/Kbl ~ •^/2. 



3.1.4- 7 and new physics Finally, I would to make a short note on new physics 
effects on determination of weak phase 7. There can be two independent approach to 
find out new physics beyond the Standard Model, if it exists. 

• We can assume the unitarity of CKM matrix. In this case, new physics effects can 
only come out from new virtual particles or through new interactions in penguin 
or box diagrams in B meson decays. If this is the case, all the methods, which I 
described, will give the exactly same 7. 

• We can generalize CKM matrix to the non-unitary matrix. In this case, new physics 
effects can appear even in tree diagram decays. And the values of 7 extracted from 
each method can be different. Therefore, I will describe in more detail for this 
second case. 

In fact, in models beyond the SM, the CKM matrix may not be unitary; for instance, in a 
model with an extra down quark singlet (or more than one), or an extra up quark singlet, 
the CKM matrix is no longer unitary |Q . If the unitarity constraint of the CKM matrix 
is removed, the generalized CKM matrix possesses 13 independent parameters (after 
absorbing 5 phases to quark fields) - it consists of 9 real parameters and 4 independent 
phase angles. The generalized CKM matrix can be parametrized as |1^ 




(28) 

Then, the GLW method would measure the angle {'y — 6), where 

7 = -5i3 and = 622, (29) 

instead of 7. The ADS method would still measure 7, but GR/JK method would 
measure (7 — 6) rather than 7, and KO method would measure (7 + 0) 
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3.2. The unitarity triangle and new physics 

Paolo Gambino, CERN 



3.2.1. Introduction The determination of Cabibbo-Kobayashi-Maskawa matrix has 
enormously improved in the last few years |l|]. All data agree remarkably with the 
Standard Model (SM) within experimental and theoretical errors. Only the recent direct 
measurement of sin 2/? from at BaBar and Belle [0, seem to show a mild 

discrepancy with SM expectations and to hint to new physics. In the case these results 
are confirmed, the main theoretical interest will be in discriminating among different 
new physics scenarios. It is therefore important to have a clear picture of the ways the 
experimental determination of the CKM matrix would be affected by physics beyond 
the SM. 

Let us briefly recall the standard way of determining the CKM matrix in the 
Wolfenstein parameterization. A and A are determined from semileptonic K and B 
decays sensitive to the elements \Vus\ and \Vch\ respectively; as these decays are tree 
level processes, this determination is to an excellent approximation independent of new 
physics. Q and fj are determined by constructing from various decays the unitarity 
triangle. 

The standard construction of this triangle involves the ratio |Vuh/K6| extracted 
from inclusive and exclusive tree level B decays and flavour changing neutral current 
processes such as B^ — B^ mixing (the mass difference {AM)^) and indirect CP violation 
in Kl decays (the parameter e), both sensitive to the CKM element Vtd- There is also 
a constraint coming from the lower bound on the mass difference (AM)^ describing 
B^ — B^ mixing. In the case of B^^ — B^^ mixings {AM)d,s is given by 0] 

(AM),,, = ^VBmB,jB,,FlM^wFu\Vti,,s)\' (30) 
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Here Fu is a function of and resulting from box diagrams with top quark 
exchanges, 5^ is a non-perturbative parameter, Fb is the B meson decay constant 
and r]B the short distance QCD factor common to (AM)^ and (AM)^. Similarly, the 
experimental value for e combined with the theoretical calculation of box diagrams 
describing — mixing gives the constraint for (^, ff) in the form of the following 
hyperbola [Q: 

[(1 - -Q)A^'n^Fu + Pc{s)\ A^Bk = 0.226 . (31) 

Here Bk is a non-perturbative parameter analogous to Bb^^, f]2 is a short distance QCD 
correction, F^t is the same function present in (|30|) , and Pds) = 0.31 ±0.05 summarizes 
charm-charm and charm-top contributions. 

Combining these and other observables which will become available in the future 
in a global fit, one determines the range of values of {g,fi) consistent with the data |jl|]. 
If the SM is correct all these measurements will result in a unique value of {g,fi). 

This procedure of testing the SM can be applied to its extensions as well. The 
determination of the unitarity triangle is affected by new physics in the following ways: 

(i) New contributions to Fu and to similar short distance functions entering rare decays 
that in the SM depend only on and M\y. These new contributions depend on 
masses and couplings of the new particles. 

(ii) New contributions that are not proportional to the same combination of CKM 
elements as the SM top contribution (disruption of the GIM cancellations). This 
means, e.g., new contributions to Pc in eq. ( pTD or new contributions to {AM)d,s 
proportional to | V,*rf(s)^cfeP- 

(iii) New phases beyond the CKM one. For instance, the CP asymmetry in i? — ip^s 
will no longer measure f3 but (5 + O^p where 9]^p is a new phase. 

(iv) New local operators contributing to the relevant amplitudes beyond those present 
in the SM, e.g., with different chirality. This would introduce additional non- 
perturbative factors Bi and new box and penguin functions. 

It is clear from (|30D and (0) that any modification of the function Fu will change 
the values of the extracted {g,f])- A recent analysis of this type in the MSSM can 
be found in p. The presence of new physics and of new phases will be signaled by 
inconsistencies in the {g,f]) plane when observables are calculated in the SM. 

3.2.2. The universal unitarity triangle While in principle a global fit of all 
experimental data can be used to test the SM and its extensions it is desirable to 
develop strategies which allow to make these tests in a transparent manner. In order to 
sort out which kind of new physics is responsible for the possible inconsistencies between 
different observables, it is useful to introduce a broad class of models very similar 
to the SM, that do not have any new operators beyond those present in the SM and in 
which all fiavour changing transitions are governed by the CKM matrix and there is no 
phase other than the CKM phase. Furthermore, we ask that in these models the only 
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sizable new contributions be proportional to the same CKM parameters as the SM top 
contributions. 

This class of models represent the slightest modification of the SM, or the most 
drastic simplification among new physics scenarios. It includes the Two Higgs Doublet 
models 1 and II and the MSSM with minimal fiavour and CP violation (MFV) — see 
1^. The kind of new physics it characterizes affects the determination of the unitarity 
triangle only through (i) above, namely through the function Fu. Therefore, the CKM 
parameters extracted from a set of data independent of the loop functions like Fu are 
universal in this class of models. Correspondingly, there exists a universal unitarity 
triangle (UUT) [0. In fact, there are as many UUT as are the sets of observables 
meeting the above requirement. 

For example, from (^) one finds that the ratio 



\Va\ _ , /(AMX, ^ , _ Fb.\IBb. 



depends only on the measurable quantities {AM)d,s, t^b^s ^^"^ non-perturbative 
parameter ^. To good accuracy this fixes one of the two unconstrained sides of the 
unitarity triangle, Rt = k/A, independently of new parameters characteristic for a given 
model. There are other quantities which may allow a clean measurement of Rt within 
our class of extensions of the SM, like the ratios of branching ratios for B Xdvv over 
B —>■ Xgiyiy Notice that the hadronic uncertainties in ([3^ ) are only due to SU{3) 



breaking effects, which should be eventually known well from lattice simulations. 

Having measured Rt from (|32|) , in order to complete the determination of the UUT 
one can use sin 2/? extracted either from the CP asymmetry in — > 'ipKs 0, ^ or 
from K — > TTuu decays ||^. Both extractions of sin 2/? are to an excellent accuracy 
independent of the new physics parameters and hadronic uncertainties have been found 
to be negligibly small B . An alternative to sin 2f5 is represented by a measurement of 



7 from tree level decays (see [IC] and refs. therein). Another possibility is to use the 
measurement of \/^~+^ by means of |Kib|/|Kfe| but this strategy suffers from much 
larger hadronic uncertainties. 

In our class of models all these different methods determine the values of f] and 
g independently of the parameters of new physics. Using them, one can calculate e, 
e'/e, {AM)d, (AM)^ and the BRs for rare decays As these quantities depend on 



the parameters characteristic for a given model, the results for the SM, the MSSM and 
other models of this class will generally differ from each other. Comparison with data 
will then single out models within our class or exclude them all, in the case new physics 
goes beyond the simple paradigm we have considered. 

3.2.3. Lower bound on sin 2/3 An interesting application of the strategy described 



so far is the derivation of a lower bound on sin 2/? in the models of our class [12 . Direct 



measurements of sin 2/3 from B^ lead to sin 2/? = 0.42 ±0.24 ^, which is only 
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marginally consistent with the SM fit |I| . Although it is certainly too soon to conclude 
that a signal of new physics has been detected, it is intriguing to speculate what kind of 
new physics could be responsible for the low value of sin 2(3 [|1^]. Two obvious possibility 
would be (a) a new CP violating phase and/or (6) a modification of the function Fu 
entering e and AMd^s- In our class of models only the second possibility is open and we 
can parametrize the new physics contribution by Ftt = F^^^^{1 + /), where / represents 
the variation relative to the SM. From and e given in (|30|j3l|) we can determine 



sin 2/3 indirectly and in model dependent way |T2 

1.26 



sin 2(3 



0.226 



-vPcie) 



(33) 



A^Bk 

where Rq depends only on AM^ and relative hadronic parameters, ri2 is a short distance 
QCD correction factor which depends on new physics only at NLO. A mild dependence 
on new physics comes mostly in an indirect way through f/. The unitarity of the 
CKM matrix alone places a constraint on /, —0.4 < / < 5.5. However, in a minimal 
supergravity scenario with universal soft terms (which belongs to our class of models) 
< / < 0.75 [i, while in a generic MSSM model with MFV one finds < / < 1.13 
Using this information together with the present ranges for the hadronic parameters, it 
is possible to derive lower bounds on sin 2/5 for the various scenarios [[12 
supergravity sin 2/5 > 0.53, in the MSSM with MFV sin 2/3 > 0.40 and 
model of our class — sin 2(3 > 0.34. Assuming a reduction by a factor 2 of all theoretical 
errors, the latter absolute bound becomes much more stringent, sin 2/? > 048. 



: in minimal 
for a generic 



3.2.4- Conclusions In summary, the UUT provides a transparent strategy to 
distinguish between models belonging to the class we have considered and to search for 
physics beyond the SM. It usefully decouples the determination of the CKM matrix from 
that of new physics parameters, in a way essentially free from hadronic uncertainties. A 
very simple and clear application of this strategy leads to the derivation of a stringent 
lower bound on sin 2(3 in a broad class of models. If the preliminary result of BaBar, 
sin 2/5 = 0.12 ±0.38, were confirmed with smaller errors, the whole class of models with 
MFV would be ruled out. 
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3.3. CP-violation, the CKM matrix and new physics 

D. Wyler, University of Zurich 



3.3.1. Introduction Observation of novel phenomena often paves the way to new 
physics. For instance, (3 decays, parity and flavor violation required the existence of 
a new force, the weak interactions. At present, it is often thought that CP-violation 
could signal new physics beyond the standard model. Although the latter can indeed 
account for the observed effects Jj] (even e'/e may be described by the standard model) 
its predictions are not well tested (compared to physics at LEP) and therefore a 
comprehensive study of CP- violation experiments is important. As sketched in figure 1, 
CP- violation manifests itself in many areas; only a comparison between them can 
determine the correct description. In the standard model, all CP- violation resides 




Figure 16. CP-Violation 



in the CKM matrix which describes the couplings of the W-bosons to the quarks of 
different charges. Thus all appreciable CP- violation occurs within flavor physics. Thus, 
one obvious strategy to search for new forces and particles would be to look for non- 
zero CP- violating effects in electric dipole moments or asymmetries in nuclear reactions. 
Unfortunately, the effects of new physics are judged to be quite small (apart from the 

^ a notable exception is the baryon asymmetry in the universe 
+ I do not discuss the so-called 9 term 
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dipole moments). Therefore more chance is given to the flavor sector instead, that is the 
physics of Kaons and mostly B-mesons. For a recent extensive review of CP- violation, 
seeref. {^). 

The unitarity of the CKM matrix 



(v. 



V 



\ 



ud 

Vtd 



Vus 
Vts 



Vub 



\ 



J 



(34) 



implies among others the triangle relation 

KdKb + KdK; + Vuv:, = (35) 

which relates observable products of matrix elements and gives stringent tests of the 
validity of the standard model. Using the Wolfenstein parameterization, and scaling as 
usual the bottom side to one, we can write for the other sides of the scaled triangle 



Rh 



ud ub 



g + ir] 



Rt 



V,,V:, = l-{-Q + in). (36) 



Here, following ref. 

^ = ^(1 - 



], the quantities 

_) f] = 7]{l-—) 

are introduced to take into account even higher powers of A. 



(37) 
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Figure 17. Unitarity triangle in the complex (fp, 77) plane. 

An elaborate analysis of superallowed (3 decay, semileptonic Kaon and Z?-meson 
decays and decays of B mesons into charmed and charmless flnal states yields 0] 

= 0.9736 ±0.001 = 1.010 ±0.16 

= 0.2205 ±0.0018 = 0.224 ±0.016 (38) 
Vub = 0.04 ± 0.002 Vcb = 0.0036 ± 0.006 
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These are (apart from corrections) all tree-level processes and therefore thought to 
be governed by the standard model They are however not sufficient to check unitarity 
(unless very precise data from t decays would be available, or if the sum of the squares 
would be significantly away from 1). 

Further input comes from loop-induced observables. They can be calculated within 
perturbation theory and input from hadronic physics. While the former are rather 
reliable and usually give results accurate to 10 percent or so, the latter are generally 
difficult to estimate. One usually considers the Kaon-mixing quantity ex-, the mass 
difference of the B and the B mesons (and also of the B^ and B^ mesons). This analysis 
has resulted in the range of values for the three angles a, (3 and 7 of the unitary triangle 
and its sides. The hadronic uncertainties are summarized in and are reflected by 

= 0.39 ±0.07 = 0.98 + 0.04- 0.22 (39) 

and by §, |, ^ 

(sin 2/?)sM = 0.75 ±0.20. (40) 

The new results this summer concern the angle (3. It was found that the coefficient a of 
sin{AMBj) in the asymmetry for B J/'^Ks is 

a = 0.79±0.4(CDF)[8] (41) 
a = 0.45 ±0.4(5e//e)|| (42) 
a = 0. 12 ± OA^BaBar) |T0| (43) 

In the standard model, one has a = sin{2(3)] comparing eqs. (^) and (10) 
we see a surprising inconsistency. Of course, this is a preliminary result, and may 
disappear as experiments collect more statistics. However, it makes it mandatory to 
investigate CP- violation in a (standard) model independent way. Unless CP- violation 
within the standard model is grossly wrong, this program essentially amounts to making 
many measurements and extracting discrepancies between quantities thought to be the 
same in the standard model. Many authors have discussed this situation; see e.g. 



11, 12, 13, 14, 15 



3.3.2. A more general framework New physics may affect every process. Because 
the standard model describes the most important weak decays, we will assume that it 
accounts for semileptonic and tree-level quark decays, at least to the required accuracy. 
This assumption can be tested, by investigating the consistency of different semileptonic 
decays, bounds from LEP etc. As an example consider the strengths of the effective 
Hamiltonians 

Heff = GF{cLl^.hL){-SLl^CL) (44) 
Heff = GF{uLl^.hL){sLYuL). (45) 

In the standard model, they are difficult to estimate. One usually considers the Kaon- 
mixing quantity ex-, the mass difference of the B and the B mesons (and also of the 
* of course, the small b ^ u transition could be due to new physics 
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Bs and mesons). On the other hand, a new neutral intermediate boson, say Z\ 
may exist, coupled to the current {sLl^hi). and {cl'~^^cl). If it also couples to quark 
and lepton pairs, such as (ul'J^iUl) and {cl'J^iCl), it would contribute to the above 
interaction, to Bg mixing, to B^ l^l^ etc. If the couplings are the same for all these 
pairs, the effective strength would be the same for the two terms in eqs. ( PD and (12). 
Therefore a new Z'-mediated interaction would induce a deviation from the standard 
model result that the couplings of the two interactions have a relative strength of A^. 
Thus detailed studies could in principle also test the first assumption. But of course, 
there are various experimental and theoretical difficulties to overcome before one will 
obtain accurate enough results. 



^ — « >- 



iy 





interfere 

B -> DK B^-> D D s 

Figure 18. two quark diagrams wliose interference gives 7 

From fig. 3 we see that the determination of the angle 7 from tree level processes 
involves the interference of amplitudes proportional to Vub and Vub respectively. This 
is achieved in processes where the two diagrams of fig. 3 contribute. A well known 



example are the decays B — > DK [|16|, more recently the advantage of Be DDg 
was stressed [|^. The idea is the same as in the previous papers on i? — > DK : One 
needs to measure the six amplitudes shown in Fig. 4. One usually considers the Kaon- 
mixing quantity e^, the mass difference of the B and the B mesons (and also of the 
Bg and Bs mesons), the sides of the triangles in Fig. 3 are of similar length and an 
extraction of 7 seems possible with the 10^° or so i?c-niesons expected at LHC. This 
method does not suffer from hadronic uncertainties. 

The experimental difficulties associated with these decays have lead to other 
possibilities. The decays B — > Kit are sensitive to the interference of the tree level 
diagram (with Vub) and the penguin diagram. This also yields the angle 7 if the penguin 



graph has no extra phase. This decay has been discussed by many people [T^ 



A third possibility that was discussed are the decays 5° D'^n^ [^. The usual 
mixing-decay formalism yields for the time dependent asymmetries the coefficients 

a ~ Im{e-'^'^'^^'^-^^^)const (46) 
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A(Bt ^ DjD^) - A{B- ^ D;D°) 



Figure 19. The extraction of 7 from Df{D", D'^,D^} decays. 



a ^ Im{e-'^'^*^'^^^^)/ const. (47) 

where const is an unknown hadronic number. It cancels in the product which then 
yields the combination 

20™x + 7 (48) 

The BB mixing angle 0mjx can be determined as usual from the decay B —>■ J/'^Kg. 

The other angles of the triangle cannot be determined by a tree level analysis 
independently. But we see, that the tree level analysis allows difficult to estimate. One 
usually considers the Kaon- mixing quantity ex, the mass difference of the B and the B 
mesons (and also of the Bg and Bg mesons), to determine the unitary triangle of the 
standard model. It yields, in principle, also the unknown side Rt and the angle (3. Any 
further independent measurement of these quantities checks the standard model with 
high accuracy, but it requires loop effects. 



3.3.3. New physics: phenomenology Among the CP- violating observables, the 
mixing-decay asymmetry is the cleanest theoretically ||21|. It is therefore reasonable 
to start an investigation of new physics with this quantity. Denoting the coefficient of 
sin(Amt) by a, one has in general 

aM^F = Im{{^)M-{^)F) (49) 
q a q 

where (|) are the mixing parameters and a, a the amplitudes for M ^ F and M ^ F, 
respectively. 

Setting for the S-meson mixing element M12 

Mi2 = rV*'^""e2^^|Mf2^| (50) 

to account for a possible new phase and magnitude of the mixing, the asymmetry 
coefficient is given in the table below: 
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The phase takes into account a possible new phase in the decay. The entries left out 
receive possibly sizeable contributions from penguin diagrams and cannot be brought to 
the simple form. This result tells us that comparing the different asymmetries, we can 
check the consistency of the standard model and determine the phases of new physics. 

New physics will also influence other CP- violating observables, such as the direct 
asymmetries of, say charged B- meson decays. In cases such as -B ^ J^tt, where 
the asymmetry is small in the standard model new physics may give rise to sizeable 
asymmetries. Of course, one need to continue the search for these, but because of the 
difficulty of calculating direct asymmetries, only quantitative statements are possible. 



3.3.4- New physics: analysis If new physics is associated with a scale A much above 
the weak scale (~ M]^), the total Lagrangian density may be written in the form []22| 

C = C''' + J2^^^^'^ (51) 
where the Oi are operators of dimension six induced by new physics and their 
coefficients di are of order(l/A^). This 'effective' Lagrangian is not renormalizable, and 
therefore one usually uses the new operators only at tree level (see a discussion by). 
The CP- violation induced by effective operators OfP can in most cases only be seen 
when they are in loops, because the imaginary part (discontinuity) of the corresponding 
Feynman graph is responsible for CP- asymmetry. || At low energies we then have an 
effective Hamiltonian 

n = J2c.Or + j:d^Or (52) 
The amplitudes for a process I ^ F and the CP conjugated one I F then are 

A{I ^F) = J2 c^Rj + ^hf^' + ^K^J- + 'h)'''' (53) 
where R and I are the dispersive and absorptive parts of the matrix elements. For the 
charge-conjugated process we have similarly 

A{I ^F) = Y, c*{Rj + + + '^i)'''' (54) 

When we calculate the CP-violating asymmetry a ~ {\A{I P)p — \ A{I P)P), we 
obtain in leading order in QCD and in NP 

a ~ Im{cd*){R'^n^'' - P^^J^*0- (55) 
ji an exception is the electric dipole moment 
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R^^ is a (finite) tree level amplitude, however also the loop I^^ is finite. Therefore 
the problems associated with a the non-renormalizable theory ^ diOf^ disappear and 
exact predictions are indeed possible for the the CP-violating asymmetry. Therefore, an 
analysis of the effects of new operators is possible also at for CP-violating asymmetries, 
and not just at tree level! 



3.3.5. New physics: models Virtually any model beyond the standard one carries 
new sources for flavour and CP-violations. It is therefore more economical to look at 
them in increasing complexity. 

The simplest one are the minimal flavour violating ones (MFV) where all sources 
of flavour violation reside in the CKM matrix. This results in many cases in a simple 
modiflcation of the coefficients in the usual loop expressions. However, there still is a 
unitary triangle, but its sizes and angles may change. It was analyzed by Ali and London 
0; recently Buras and Buras |jT5|] found a clever lower bound on sin{2(3). The idea is 
simple. For both e and the P-meson mass difference, the standard model contribution 
consists mostly oi — W — t — t box diagram; its value might be denoted by Fu- The 
MFV modify this to 

Ftt = 5oK) (1 + /) . (56) 
Then we can write for ek 

eK--n [(1 - Q)A%Fu + Pc{e)\ A^Bk (57) 
while the P-meson mass difference yields the relation 

Po 1 



where 



Rt = 1.26 



Pn 



With 



one gets 



sin 2/5 



sin 2(3 




1.26 



0.226 



R'oV2 VA^Bk 



(58) 



(59) 



(60) 



(61) 



Since unitarity implies t] < Rf, , there exists a lower bound on sin 2(3. A careful numerical 
analysis implies 

f] > 0.34. (62) 

The lower bound in fact corresponds to a Fu which is three times larger than the 
standard model value. 

Supersymmetry is a attractive candidate for new physics. In general, there are many 
new CP- violating phases. Since they can directly affect observables such a the electric 
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dipole moment, it is natural to take them to be small (approximate CP- violation, []1|). 
In this situation, also CP- violating effects in the P-system are small. This implies a 
small angle j3. This is in contrast to the standard model, where the ffavour structure 
suppresses CP-violation. 

The problem with this scheme is that it is hard to get ek right and that e'/ e tends 
to be to small. 

Similarly, models with left-right symmetry tend to have small CP- violating phases, 
thus the effects tend to be small also. 



3.3.6. CP-violation in P-mesons In the standard model, CP- violation is small 
in the P-System. This is partly due to the rather large tree-level decay rates and 
small coupling of the third generation. Therefore one would expect new physics CP- 
violation mostly in the mixing (see [0] for a more detailed discussion). Recent studies 
of time-dependent decay rates of P° K^ti~ by the CLEO collaboration and 
measurements of the combination of P° — * K^K~ and P° — > K^t:^ rates by the 
FOCUS collaboration gave first information on the mixing. 
As usual, one define the mixing quantities 

m2-mi r2-ri 

CP-violation in the mixing is defined by the angle (j). The experiments find that the 
quantity ycos4> is significantly larger than the expectation in the standard model. The 
errors being large, this result is not yet significant, but it shows the potential of P-meson 
physics. 



3.3.7. i^-physics Finally let me mention i^'-physics. Of course, efforts continue in 
calculating e'/e. However, the rare decays n^uu and —>■ tt^uj? provide a 

theoretically clean way to measure (in the standard model) \Vtd\ and ImVtd Clearly, 
this can be used as a test of the unitary triangle, however the measurement of the neutral 
decays is not easy and probably many years away. 



3.3.8. Conclusions The new results on sin 2/3 are surprising; they may indicate 
a failure of the standard model. Several parameters have to be stretched beyond 
their reasonable values to account for them. One can modify the standard model to 
accommodate the small value of sin 2/3, but it is not clear that these modifications are 
consistent. 

Nevertheless, the result brings back the (old) view, that a (standard) model 
independent and broad analysis of CP- violation is required in order to fully understand 
this phenomenon and the need for new interaction. this implies in particular 
measurements of many decay channels. 

I have sketched strategies to determine the source of CP- violation for the case 
that the standard model accounts for tree level processes and given a phenomenological 
framework to calculate the effects of new operators. Needless to say that all of this will 
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take many years of hard work on both the experimental and the theoretical side and 
that also less perfect measurements have to be pursued. 
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4. What can we learn from e'/e? 

CP violation in the kaon system 

V. Martin, University of Edinburgh 

CP violation was first observed in the neutral kaon system by Christenson, Cronin, 
Fitch and Turlay in 1964 when they observed the long-lived neutral kaon decay into 
a TT+TT" state. To date, the neutral kaon system is the only place where CP violation 
has been conclusively observed. 

In the neutral kaon system CP violation is classified into two types: indirect 
and direct, parameterised by e and e' respectively. Indirect CP violation may be 
observed through the asymmetry in the mixing of the neutral kaons. Direct CP 
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Figure 20. Recent theoretical estimates of TZe{e'/e) in chronological order 
[H H' 0! 1' ^- The red/dark-grey error bars correspond to a Gaussian treatment of 
the input parameters. The yellow/light-grey error bars correspond to a flat scanning 
of all the input parameters. The cross lines corresponds to the most favoured value. 
Where there are two different results from one group, these correspond to different 
renormalisation schemes. 



violation may be observed through the asymmetry in the decays of the neutral kaons. 
As direct CP violation is a much smaller effect than indirect CP violation, the ratio 
TZe{e'/e) ^ e'/e is usually considered as the measurement of direct CP violation. The 
last generation of Tle{e'/e) experiments announced final results in 1993. The NA31 
experiment at CERN had a result which was compatible with direct CP violation: 
ne{e'/e) = (23.0 dz 6.5) x 10"^ ||, whereas the E731 experiment at FNAL found a result 
compatible with no direct CP violation: 7^e(e7e) = (7.4 ± 5.9) x 10"^ 0. 

This disagreement prompted a new generation of direct CP violation experiments 
with the aim of measuring TZe{e' /e) at the 10~^ level. 

Theoretical estimates for the value of TZe{e'/e) are based on QCD calculations. 
The short distance contributions have been reliably calculated using perturbative- 
QCD 0]. Various methods have been used to estimate the long distance effects, such 
as chiral perturbation theory and lattice QCD. Figure ^ shows some recent theoretical 
predictions of 7le{e'/e). 

Experimentally, the parameter TZe{e'/e) is obtained by measuring the double ratio 
of decay rates (R) of the long-lived and short-lived neutral kaons Kl and Ks into two 
pion final states: 

_ r(KL ^ 7r°7r°)/r(Ks ^ vr^vr") .^^/^ (f^A^ 

^ = r(K,^vr+vr-)/r(Ks^vr+vr-) = ' " ^^^^^ ^^'^ 

Currently there are three dedicated experiments analysing or taking data with 
the aim of measuring 'Re{e'/e): KTEV at FNAL, NA48 at CERN and KLOE at the 
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DA$NE $ factory at Frascati. Each of these experiments was designed specifically to 
measure TZe{e' /e). In particular, they were designed to record high statistics and to 
minimise systematic errors in the measurement of TZe{e'/e). The following discussion 
here will focus on the measurements from NA48 and KTeV, which were presented at 
this workshop. 

In the measurement of the double ratio R, many systematic effects cancel, such as 
the number of Ks and Kl particles produced, and the efficiency of detecting the charged 
and neutral final states. 

The NA48 and KTeV experiments use a similar experimental set-up, aiming to 
exploit many of the cancellations in the measurement of R. In order to be able to 
detect all four TZe{e' /e) decay modes simultaneously, both experiments employ a double 
beam technique: two neutral kaon beams are produced at different distances from the 
detector. Decays from the beam with the longer path length are predominantly Kl 
decays whereas decays from beam with the shorter path length are predominantly Kg. 

KTeV uses two identical, parallel, neutral kaon beams. A plastic scintillator 
regenerator is placed in one of the beams to regenerate Kg particles. The position of the 
decay products in the detectors can be used to determine whether the decay came from 
the vacuum or regenerator beam. The two NA48 beams are created by coUiding the 
SPS protons on identical beryllium targets. The so-called Kl and Kg targets are placed 
respectively 240 m and 122 m upstream of the detector. The targets are arranged such 
that the two beams converge, meeting at the centre of the electromagnetic calorimeter. 
The times of the protons in transport to the Kg target arc recorded by the proton 
tagging detector. The proton times can be compared offline to the time of the decay 
products in the detectors to determine which of the beams the decay came from. 

To detect the decay products of the neutral kaons both the NA48 and KTeV 
experiments have: a magnetic spectrometer to detect charged tracks and reconstruct 
K — > TT'^TT" decays; a high-resolution electromagnetic calorimeter to detect photons from 
7r° — >• 27 decays and a muon detector to reject background from Kl —>■ Ti'^fi^u^ decays. 

NA48 uses a quasi-homogeneous Liquid Krypton calorimeter which contains 9 m^ 
of liquid krypton kept at 121 K. KTeV has a Cesium-Iodide crystal calorimeter, made 
of 3,100 blocks of pure Cesium-Iodide. The quality of these two calorimeters is reflected 
in the resolution in measuring reconstructed masses from the K — > tt^tt" — 47 decays. 
NA48 achieves a resolution on the reconstructed 77 mass of 1.1 MeV/c^, and KTeV 
achieves a resolution on the reconstructed kaon mass of 1.7 — 1.9 MeV/c^. 

The analysis of the data from the TZe{e'/e) experiments is, in principle, 
straightforward. The number of events in each of the four decay modes in counted 
and a initial value of R is obtained. Corrections must be made to this value of R 
for backgrounds, trigger efficiencies, detector energy scale, beam scattering, accidental 
effects and detector acceptance. In addition both experiments make corrections for 
effects due to mis-identifying which of the beams the decay came from. Most of these 
corrections are small due to the cancellations present in R. 

The largest potential correction to R is due to the differential acceptance of the 
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experiment 


Kl ^ ttOttO 


Ks ^ TX^Tl^ 


Kl — ^ TT + TT 


Ks TT+TT 


NA48 
KTeV 


1.80 X 10^ 
0.86 X 10^ 


1.14 X 10^ 
1.42 X 10^ 


7.46 X 10^ 
2.61 X 10^ 


4.87 X 10^ 
4.51 X 10^ 



Table 11. Numbers of decays for each TZe{e' /e) decay mode collected by NA48 in 
1998 and by KTeV for their preliminary result. The numbers have been corrected 
for background decays. The KTeV numbers actually refer to the numbers of decays 
counted from the vacuum beam (Kl) and the regenerator beam (Ks). 



four decay modes. To reduce this correction NA48 uses the technique of weighting 
the Kl events so that the decay distribution of weighted Kl events is the same as the 
(unweighted) Ks events. This reduces the correction to R due to the acceptance from 
around 13% to around 3 x 10~^, but leads to a reduction of statistical precision of 
35 — 40%. KTeV instead relies on a detailed Monte Carlo simulation of the experiment 
to calculate the acceptance. They calculate a correction to R of around 4%, which they 
measure with a precision of 13.4 x 10^^ [p!0| . 

KTeV collected data for the 7le{e'/e) measurement during 1996, 1997 and 
1999. A preliminary result based on some of the data taken during 1996 
and 1997 is described in [1^. The result of that analysis is 7le{e'/e) = 
(28.0 ± 3.0(stat) d= 2.6(syst) ± 1.0(MC stat)) x 10"^ 

NA48 has collected data for the TZe{e'/e) measurement in the years 1997-99 and 
will take more data in 2001. Preliminary results based on the data taken in 1997 
and 1998 have been announced ||Tl|, |12|. The combined result from these data sets 
is: TZe{e' /e) = (14.0 ± 4.3) x 10~^, where the statistical and systematic errors have 
been added in quadrature, and the small correlation between systematic effects has 
been taken into account. The error on this result is currently dominated by systematic 
effects, but as the systematic error includes some statistical component, the error will 
improve significantly when more data is analysed. 

Table [ll| shows the number decays in each TZe{e'/e) mode analyzed by NA48 (in 
1998) and KTeV for these preliminary results. 

Figure 21 shows the preliminary NA48 and KTeV results along with the results 
from E731 and NA31. A weighted average of these four results gives 7?.e(e'/e)giobai = 
(19.2 ± 2.5) X 10""^ with a xV^df = 10.4/3. The yellow/light-grey band on figure |T 
shows the ila allowed region of this average. As the P^r degree of freedom on the 
average is larger than 1, the approach proposed by the PDG can be used to rescale the 
error by {x^/ndfy/^ = 1.8, resulting in 7^e(eVe) global = (19.5 ± 4.6) x 10"^ 

The theoretical calculations shown in figure |2^ generally predict smaller values 
for lZe{e'/e) than the experimental results, leading to some speculation that the 
measurement of TZe{e' /e) may be a sign of new physics. However due to the relatively 
large uncertainties reported by both theoretical calculations and experiments it is 
obviously too early to draw any conclusions on this issue. 

Concluding on CP violation in kaon physics: the experimental picture is now much 
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NA31 E731 KTEV NA48 



Figure 21. Recent measurements of TZe{e' /e) |||, ||, |T^, |T2j. Where both a statistical 
and systematic error is given, these have been added in quadrature. The yehow/hght- 
grey band shows the ztlcr allowed region of the average of these results. 



clearer than it was in 1993. Direct CP violation in the neutral kaon system is now firmly 
established. However, due to the large variation in the measurements, it is too early to 
settle on an exact value for TZe{e' /e). New, more precise experimental results will be 
produced by NA48, KTeV and KLOE in the future, allowing TZe{e' /t) to be known with 
an uncertainty of around 1 x 10^"^. 
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4-2. Supersymmetry predictions for e' /e 

Shaahan Khalil, University of Sussex 



4-2.1. Introduction The most recent results of e'/e , which measures the size of the 
direct CP violation in — ^ vrvr, reported by KTeV and NA48 lead to a world 
average of Re e'/e =(21.4±4.0) x 10~^ 0. This result is higher than the Standard Model 
(SM) predictions 0|, opening the way to the interpretation that it may be a signal of 
new physics beyond the SM. The SM predictions for e'/e suffer from large theoretical 
uncertainties such that one can not draw a definite conclusion if this observed high 
value of e'/e can be accommodated in the SM. In any case, one may wonder if the 
supersymmetry (SUSY) can be responsible for enhancing e'/e . 

In the minimal supersymmetric extension of the SM (MSSM) there is no way of 
generating a sizable SUSY contribution to e'/e even if one assume that the SUSY CP 
violating phases are large and the electric dipole moments (EDM) of the electron and 
neutron are less than the experimental bounds due to the cancellation between the 
different contributions. This is mainly due to the assumption of universal boundary 
conditions of the soft-breaking terms [|, 0, ||, P]. It has been shown that, without new 
fiavor structure beyond the usual Yukawa couplings, general SUSY models with phases 
of the soft terms of order 0(1) (but with a vanishing CKM phase Sqkm = 0) can not give 
a sizeable contribution to the CP violating processes 0, ||, |^, This means that the 
presence of non-universal soft breaking terms besides large SUSY phases is crucial to 
enhance these CP violation effects. In agreement with this, it has been explicitly shown 
that contributions to ex are small within the dilaton-dominated SUSY breaking of the 
weakly coupled heterotic string model [1^, where A-terms as well as gaugino masses 
are universal. On the other hand, it is well-known that the strict universality in the 
soft breaking sector is a strong assumption not realized in many supergravity and string 
inspired models |TT|, |T2[. All these arguments indicate not only that the presence of 
non-universal soft terms can solve the problem of too large contributions to EDMs but 
also that it allows for large SUSY contributions in CP violation experiments. Hence, in 
this work we will follow this avenue and analyze the effects of non-universal soft terms 



in CP violation in the i^'-system [13, Q 



4-2.2. CP violation in minimal supergravity model It is well known that in 
SUSY models there are new possibilities for CP violation. In particular, the soft SUSY 
breaking terms contain several parameters that may be complex, as can also be the 
/i-parameter. In the minimal supergravity model there are only two new CP-violating 
phases. This can be seen as follows. The parameters M, A and B and fi can be complex. 
But of these four phases only two are physical. First, by an R-rotation with R-charge 
Qr = 1 for lepton and quark superfields and Qr = for the vector and the Higgs 
superfields, the gaugino mass parameter M can be made real. Second, Bfi can be made 
real by a change of phase of the Higgs superfield. This ensures that the Higgs vacuum 
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expectation values are real. The remaining phases cannot be defined away and violate 
CP. One is in A = Aqc'"^^ and the other m. B — B^e^'^^ . The ji parameter then has a 
fixed phase ^l — //oe~*'^^ . In any phase convention 

0A = arg(AM*), 0^ = arg(SM*). (65) 

These phases can cause at one loop level an electric dipole moment (EDM) for the 
quarks and leptons, and therefore also for the neutron. It has been known for a long 
time that in the SUSY models the contributions to the neutron electric dipole moment 
are larger than the experimental limit 6.3 x 10^^^ e cm unless either the new 'SUSY 
phases' are tuned to be of order 10^'^, or the SUSY masses are of order a TeV. Such 
small phases can not generate sizable CP violation. Also they constitute a fine tuning. 
This is known as " SUSY CP problem" . However, in the last few it has been suggested 
that a natural cancellation mechanism exists whereby the electric dipole moment of the 
neutron may be made small without such fine-tuning. In this case large SUSY phases 
are expected and still satisfy experimental bounds on the values of EDM of the electron 
and neutron. 

In this section we will study the effect of these phases in CP violation observables 
as e and e' /e . We assume that 6ckm = to maximize this effect. The value of the 
indirect CP violation in the Kaon decays, e, is defined as e = e'^^lmMu / V^Attik , where 
ArriK = 2Re{K^\Heff\K'^) = 3.52 x 10"^^ GeV. The amplitude Mu = {K'^\Heff\K^) . 
The relevant supersymmetric contributions to — are the gluino and the chargino 
contributions, {i.e., the transition proceeds through box diagrams exchanging gluino- 
squarks and chargino-squarks) . It is usually expected that the gluino is the dominant 
contribution. However, as we will show, it is impossible in the case of degenerate A- 
terms that the gluino gives any significant contribution to e when the CKM matrix is 
taken to be real even with large phase of A. The amplitude of the gluino contribution 
is given in terms of the mass insertion 5ab defined by 5ab = Aab/'^'^ where m is an 
average sfermion mass and A is off-diagonal terms in the sfermion mass matrices. The 
mass insertion to accomplish the transition from diL to dji^ {i,j are flavor indices) is 
given by 



d ).. 
LL)i3 



(A' 
(AL)., 



(A 



d V. 
RLJtJ 



87r2 
1 

8^ 
1 

8^ 



{37h' + \X\ 



v"^ siv? P cos P 
sin^ P cos P 



ln(^)X, 



Ml 



w 



where X 



iKRh = 0, 

Ad — fj. tan l3. 



(66) 

(67) 

(68) 
(69) 

It is clear that Ajj in general are complex due to the 



w 



complexity of the CKM matrix, the trilinear coupling A and fx parameter. Here we 
assume the vanishing of Sckm to analyze the effect of the SUSY phases. We notice that 
(A^^)i2 is proportional to |Xp i.e., it is real and does not contribute to e whatever the 
phase of A is. Moreover, the values of the (A^^)i2 and (A^^)i2 are proportional to rris 
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and rrict^ hence they are quite smalL Indeed in this case we find the gluino contribution 
to £ is of order 10~®. 

For the chargino contribution the amphtude is given by [1^ 

1 

-( 

-3 



(2vr)^ 

+ + ^^3(.)i^3(.))]. (70) 

The complete expression for these function can be found in Ref. [|^. For low and 



moderate values of tan j3 the value of C3 is much smaller than Ci since it is suppressed 
by the ratio of rUs to Mw- However, by neglecting the flavor mixing in the squark mass 
matrix Ci turned out to be exactly real 0. The imaginary part of Ci is associated to 
the size of the intergenerational sfermion mixings, thus it is maximal for large tan/?. In 
low tan/5 case, that we consider, the imaginary part of Ci is very small, and the gluino 
contribution is still the dominant contribution e. In particular, we found that e is of 
order 10~^, which is less than the experimental value 2.26 x 10~^. 

Now we consider the effect of these two phases {(pA and (p^) on the direct CP 
violation parameter e' /e . Similar to the case of indirect CP violation parameter e , 
in the gluino contribution the L-L transitions are almost real and the L-R transitions 
are suppressed by two up Yukawa couplings and a down quark mass. Moreover, the 
analysis of chargino contribution is also the same as in the indirect CP violation. Even 
the experimental bounds on the branching ratio of 6 ^ 37 decay impose sever constraint 
on the LR transition. Therefore we do not flnd any signiflcant SUSY CP violation effect 
in e' /e too. 



4-2.3. Non— universal soft terms and SUSY CP violation In the previous 
section, we have shown that CP violation effects are always very small in SUSY 
models with universal soft SUSY breaking terms. Recently, it has been shown that 
the non-universality of A-terms is very effective to generate large CP violation effects 
0, H, 0, |18|. In fact, the presence of non-degenerate A-terms is essential for 

enhancing the gluino contributions to e'/e through large imaginary parts of the L- 
R mass insertions, Im((5x,R)i2 and Im(5/^2.)i2- These SUSY contributions can, indeed, 
account for a sizeable part of the recently measured experimental value of e'/e [jl], @. In 
the following, we will present an explicit realization of such mechanism in the framework 
of a type I superstring inspired SUSY model. Within this model, it is possible to obtain 
non-universal soft breaking terms, i.e. scalar masses, gaugino masses and trilinear 
couplings. 

Type I string models contain nine-branes and three types of five-branes (5a, 
a = 1,2,3). If we assume that the gauge group SU (3) x U (1) is on one of the branes (9- 
brane) and the gauge group SU{2) is on another brane (5i-brane). Chiral matter fields 
correspond to open strings spanning between branes. Thus, they have non-vanishing 
quantum numbers only for the gauge groups corresponding to the branes between which 
the open string spans. For example, the chiral field corresponding to the open string 
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between the SU (3) and SU (2) branes can have non-trivial representations under both 
S'f/(3) and SU{2), while the chiral field corresponding to the open string, which starts 
and ends on the S'[/(3)-brane, should be an S'f/(2)-singlet. There is only one type 
of the open string which spans between the 9 and 5i-branes, which we denote it as 
C^^^. However, there are three types of open strings which could start and end on 
the 9-brane, that is, the Cf sectors (i=l,2,3), which corresponding to the i-th complex 
compact dimension among the three complex dimensions. If we assign the three families 
to the different Cf sectors we obtain non-universality in the right-handed sector. In 
particular, we assign the Cf sector to the third family and C| and C| to the first and 
second families respectively. Under these assumption the soft SUSY breaking terms are 
obtained, following the formulae in Ref . |[Tl|, [T^ . The gaugino masses are obtained 



M3 = Ml = V3m3/2sin^ e-^°^ 

M2 = VSms/scos^ Gie^*"!. 
While the A-terms are obtained as 

Acf = -\/3m3/2sin^ e-'"^ = -M3, 
for the coupling including Cf, i.e. the third family, 

= -y3m3/2(sin^ e'^^^ + cos^ (Gie 
for the coupling including C|, i.e. the second family and 

Ac9 = -V3m3/2{sme e'^^^ + cosO (Gie^^^i - 636 



-lai 



626-*"^)), 



-l«3 



)), 



(71) 
(72) 

(73) 

(74) 

(75) 



for the coupling including C|, i.e. the first family. Here 7/13/2 is the gravitino mass, 0:5 
and ttj are the CP phases of the F-terms of the dilaton field S and the three moduli 
fields Tj, and 9 and 6j are goldstino angles, and we have the constraint, X] = 1- 
Thus, if quark fields correspond to different Cf sectors, we have non-universal A-terms. 
Then we obtain the following A-matrix for both of the up and down sectors. 



A 



A 



A 
A 



(76) 

\ ^cl ) 

Note that the non-universality appears only for the right-handed sector. The trilinear 
SUSY breaking matrix, iY^)ij = (Y)ij{A)ij, itself is obtained 



(77) 



in matrix notation. In addition, soft scalar masses for quark doublets and the Higgs 
fields are obtained, 
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;i_rcos'^ (l-Qi)). 
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The soft scalar masses for quark singlets are obtained as 

1 -3cos2 e^). 
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(79) 
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if it corresponds to the Cf sector. 

In models with non-degenerate A-terms we have to fix the Yukawa matrices to 
completely specify the model. In fact, with universal A-terms the textures of the Yukawa 
matrices at GUT scale affect the physics at EW scale only through the quark masses and 
usual CKM matrix, since the extra parameters contained in the Yukawa matrices can be 
eliminated by unitary fields transformations. This is no longer true with non-degenerate 
y4-terms. Here, we choose our Yukawa texture to be 

y" = — ^diag (m„, mc, mt) , = — ^K"^ ■ diag (m^, m^, m^) ■ (80) 
V cos jj vsmfj 

where K is the CKM matrix. In this case one find that the mass insertion 5f^^ can be 
written as [14] 

(^LRlij =^^i ('^aAU + c-gnfg- /ie^'^nan/?) 

+ K*, CA (A^« - A*^,) + K,, K*, CA (A^. - A^.)) (81) 

where m? is an average squark mass and rrij the quark mass. This expression shows 
the main effects of the non-universal A-terms. In the first place, we can see that the 
diagonal elements are still very similar to the universal A-terms situation. Apart of 
the usual scaling with the quark mass, these fiavor-diagonal mass insertions receive 
dominant contributions from the corresponding A^-g terms plus an approximately equal 
contribution from gluino to all three generations and an identical fi term contribution. 
Hence, given that the gluino RG effects are dominant, also the phases of these terms 
tend to align with the gluino phase, as in the minimal supergravity. Therefore, EDM 
bounds constrain mainly the relative phase between /i and gluino (or chargino) and give 
a relatively weaker constraint to the relative phase between Aqs (the first generation A- 
term) and the relevant gaugino. Effects of different A(j9 in these elements are suppressed 
by squared CKM mixing angles. However, flavor-off-diagonal elements are completely 
new in this model. They do not receive signiflcant contributions from gluino nor from 
/i and so their phases are still determined by the A^s phases and, in principle, they do 
not directly contribute to EDMs . 

In flgure ^ we show the allowed values for as, 0.2 and 03 assuming ai = = 0. 
We have imposed the EDM, e and b s'-f bounds with the usual bounds on SUSY 
masses. We can see that, similarly to the minimal supergravity, ip^ is constrained to be 
very close to the gluino and chargino phases (in the plot as — 0, vr), but 02 and as are 
completely unconstrained. 

Finally, in flgure ^ we show the values of /m(5^'^^)2i versus the gluino mass in 
the same regions of parameter space and with the same constraints as in flgure ^ As 
we can see due to the effect of the off-diagonal phases a large percentage of points are 
above or close to 1 x 10~^, hence, sizeable supersymmetric contribution to e'/e can be 
expected in the presence of non-universal A-terms. 
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Figure 22. Allowed values for a2-as (open blue circles) and a^-as (red stars) 



4- 2.4- Conclusions Non-universal Supersymmetry soft breaking terms are a natural 
consequence in many supergravity or string inspired SUSY models. Moreover, the non- 
universality of the A-terms has a significant effect in the CP violation. We have shown 
that in these models a sizeable supersymmetric contribution to CP observables e' /e and 
e can be easily obtained. 
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5.1.1. Disentangling the CP eigenstates The decay rate of the mode Bg J/iprj 
as a function of time is approximately given as 

^'^(Bs^J/'^v) ^|^(o)|2 [e-r.t + 2e-f*sin(Am,t)57] , (82) 

where T = (F^ + r//)/2. The non-oscillatory part of (|82D gives the value of Ti whereas 
the oscillatory part gives 6-y. 

The decay mode Bg J/ip{i~^i~)(l){K^K~) is another candidate for the 

determination of ^7. It has the advantages of having a larger branching ratio and 
a higher efficiency of detection (since all the four final state particles are charged). 
However, the final state J/ip4' is not a pure CP eigenstate, but an admixture of CP 
odd and CP even components. Here, we argue that in spite of this, the mode J/ip(f) is 
a better candidate than J/ipr] for determining ^7. 

The most general amplitude for this decay can be written in terms of the 
polarizations ej/^, e<^ of the two vector mesons as Q: 

AiB, Jm = Ao (9^) '*J%'t-^ ^m<-'^^ ^}/^x4-p , (83) 

where is the energy of and p is the unit vector in the direction of in the J j ip rest 
frame. The superscripts L and T represent the longitudinal and transverse components 
respectively. Since the direct CP violation in this mode is negligible, the amplitudes Aq 
and A\\ are CP even whereas A^ is CP odd. 

The final states with different CP parities can be separated through their different 
angular distributions. The decay Bs —>■ J/tp{i~^i~)(f){K'^K~), however, has four final 
state particles, which implies that the angular distributions will, in general, be in terms 
of three physical angles [jl], 0]. Such a three angle analysis has been performed 0], where 
preliminary studies indicate that that the values of ^7 as low as 0.03 may be accessible, 
depending on the value of Xg = Amg/f . 

However, more the number of angular terms, more formidable the task of 
disentangling their coefficients, which are the physical quantities of interest. The method 
of angular moments [0, ^ may be of some help, but with as many as 6 angular terms in 
the complete three angle distribution, the errors in the determination of the coefficients 
are large. Moreover, it involves a fit to as many as 8 independent parameters, some of 
which seem to have strong correlations. Getting rid of these problems would increase 
the sensitivity on the value of ^7, and might bring its SM value within the domain of 
measurability. 

The angular analysis can be simplified a lot if the angles are chosen such that the 
angular distribution that separates the CP odd and even terms can be written in terms 
of a single angle. Here we show that such a angle may be defined for this decay. It 
follows on the lines of the arguments given in The transversity angle distribution 
reduces the problem to disentangling only two angular terms (as opposed to six earlier) 
and only five independent parameters (as opposed to eight). 
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5.1.2. The transversity angle In the rest frame of 0, the decay B J/ipK^K' is 
planar. Let us define the transversity axis as the one perpendicular to this decay plane. 
Let the decay plane be the x — y plane, so that the z axis is the transversity axis. The 
transversity angle is the angle made by the spin of J/ip with this axis. 

The final state {J/iIjK~^K~} in the rest frame of may be represented in the 
basis of \pk+'Pk-V^IiJz^),p^ where J^^ is the z component of the spin of J/'?/'. Clearly, 
'Pk+ = —PK-- Consider the operator TZcp, which combines charge conjugation with the 
reflection in the x — y plane: 

7^cp = Clip = CVe'^^^ (84) 

where C and V are the charge conjugation and parity transformation operators 
respectively. Since the total angular momentum of the final state is zero, from (|8^) , 
we have 

T^Cp\PK+PK-P4>Jzi,)<P= iCP)f\pK+PK-P4>Jz'4,)<t> , (85) 

where {CP)f is the CP parity of the final state. 

On the other hand, both \pK+PK-)(t> and \p^i)Jzip) (f, are separately eigenstates of TZcp, 
and TZcp commutes with any boost in the decay plane. Then, if we denote the boost 
from rest frame to J/ip rest frame by B, so that B\p^Jz^)^ = \p^Jz'i(>)j/ip, then 

T^cp\PK+PK-Pi>Jz^p)<t> = 'R'Cp\pk+PK')4> B'^TZcpBIp^J^^)^ (86) 

= {CP)^\pk+Pk-)<p B~\CP)j/^e'^^\p^J,^),ji^ (87) 

where r is the transversity of J/V", i-e. the projection of its spin along the transversity 
axis. Using (CP)^ = {CP)j/^ = +1 and the fact that B commutes with the other 
operators, we can write 

T^cp\PK+PK-P'4>Jzi)<t> = \PK+PK-)<j> e^'^'^lp^J,^)^ (88) 

= (-1)^ \pK+PK-Pi}Jzip)^ ■ (89) 

From (|85|) and (|9]), we get 

iCP)f\pK+PK-Pi,Jzij)<t, = {-'i-Y \PK+PK-PijJzi,)4> , (90) 

which shows that the final states with even and odd CP parities correspond to the 
states with the transversity of J/ip equal to and ±1 respectively. There two states can 
easily be separated by using the parity conserving decay J/ip —>■ as an analyzer, 
and studying the angular distribution in the transversity angle, the angle made by the 
decay products with the transversity axis in the rest frame of J/ip. 

The above argument may be generalized and shown to be applicable in the case 
of all decays of the form B — > AC{^ C1C2), where (a) A and C are self-conjugate 
particles with Ci and C2 spinless and charge conjugates of each other (as in the case 
under discussion), (b) A,Ci,C2 all are self-conjugate particles, or (c) A and C are 
charge conjugates of each other, with Ci and C2 spinless particles. In all such cases, a 
transversity angle can be defined such that the angular distribution separates CP odd 
and even final states P]. 
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Figure 24. The definitions of angles 9, ip, ip. Here is the transversity angle. 



5.1.3. Transversity angle distribution in J/ipcf) Let us define the angles as shown 
in Fig. |2^. Let x axis be the direction of in the J/ip rest frame, and let z axis be 
perpendicular to the decay plane of — > K~^K~ , with positive y direction chosen such 
that py{K^) > 0. Then we define {0,ip) as the decay direction of in J/tp rest frame 
and ip as the angle made by with x axis in rest frame. 
The angular distribution is then given by 



dT[Bs ^ J/ij4>] 



l\A+it)\\l + cos' 6) + ^\A4t)\' sin' 9 
o 4 



where 



d cos 9dt 
|Ao(t)P + |A||(t)|2and |^_(t)p 



(91) 



\A^{t)\' are the CP even and CP 



odd components respectively. The time evolution of these components is given by 



cos^(57)e"'^* + sin^(57)e-'«* + e-'*sin(Am,t)sin(2(57) , (92) 



\A^it)\' = \A^ 



sm'{S'j)e + cos^(57)e' 



e"^* sin(Am,t) sin(257) . (93) 



The time evolution of either one (or both of) the components lets us determine the 
values of Th, Tl and ^7. 

The value of ^7 is small (~ 0.015) in the standard model. To first order in ^7, eqs. 
(p3) and (p3D may be written as 



\A-{t)\' 



iA(o)r 



3"^^* + 2e"^* sin(Am^t)57 



'^Ht _ 2e-^* sin(Am,t)57 



(94) 
(95) 
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The time evolution (|9^ ) is exactly the one used for the determination of 57 through 
J/ipT] in (p^). Since only two angular terms are present in (0), the separation of their 
coefficients is not hard. Therefore, the accuracy obtained through these two modes 
would be similar if the number of events were the same. However, the mode J/ipip is 
expected to have almost an order of magnitude more number of events than J/iprj. 

An interesting situation arises when the measured value of ^7 is large. This is 
clearly a signal of physics beyond the SM, in particular of an extra phase appearing in 
the Bs — Bs mixing. Using the exact expressions for time evolution (eqs. ^ and ^3|) will 
give a measurement of this physically interesting phase, whereas The approximate time 



evolution (eqs. |9J and ^Sj) will fail to do so. In this case, any advantage J / i/jt] would 
have had due to its being close to a CP eigenstate is lost. The non-oscillatory part of 
(^) can no longer be used to give a clean measurement of Tl, since the time evolution 
now consists of two exponential decays with similar lifetimes. In [Q, the information 
content per event in such decays has been quantified, and it has been shown that the 
mode with angular information may have orders of magnitude more information per 
event than the mode with time information alone. Adding this to the larger number of 
expected events in J/ipcf)., the case for this mode becomes even stronger. 
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5.2. Width Difference in the Bd - Bd system 

A. Dighe, CERN, T. Hurth, CERN, C.S. Kim, Yonsei University 

Within the standard model, the difference in lifetimes of Bd mesons is CKM 
suppressed compared with the ones of Bg mesons. A rough estimate leads to: 



^0.5% , (96) 

Vcd 

where Td is the average lifetime of the light and heavy Bd mesons {B^ and B^ 
respectively). We denote these lifetimes by rf,r^ respectively, and define AF^ = 

More precisely, we get for at leading order using the standard notations |[ 

ri2 = _g|g!^[5 ml _ K,)flBs{V,,V:df (97) 

247r 3 [nil, + rudY 



+\{K, + ^)fiBsiy^,v:,f + mi + K,)fiBB'^v,,v:dVt,v:d]. m 

Ki = —0.39 and K2 = 1.25 are combinations of Wilson coefficients. Bs and 
Bb are the bag factors corresponding to the matrix elements of the operators Qs = 
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{hd)s-p{hd)s-p and Q = {bd)Y_A{bd)v-A- Using the known expression for M12 (see i.e. 
|1|]) and = 4.5 GeV, we have 



- -5.0x10- ^1.4^. 0.24 + 2.5^^J. (99) 

In the vacuum saturation approximation one has Bs/Bb = 1 at some typical hadronic 
scale. 

With Re(ri2/Mi2) — AF/Am, and the measured value Xd = Amrf/F^ = 
0.73 ± 0.05 we confirm the rough estimate (^) taking into account the large hadronic 
uncertaintities. 

No experimental measurement of AF^^ is currently available. Indeed, any new 
physics contribution can only decrease AF^ 0, thus taking it still out of the range 
of sensitivity of the B factories or at the present hadronic machines that concentrate 
on B mesons. Moreover, no motivation for its measurement (other than just measuring 
another number to check with the SM prediction) has been discussed, and hence the 
study of the lifetime difference between B^ mesons has hitherto been neglected as 
compared to the Bs mesons. The corresponding lifetime difference in the case of Bg 
is expected to be AF^/F^ ^ 15%. 

The time resolution and the high statistics expected at the LHC Q might make it 
possible to measure AF^. Here we shall discuss the possible measurement of AF^ at the 
LHC, the effect of AF^^ on the measurement of sin(2/5) through B^, J/ipKs and the 
possible resolution of the discrete ambiguity in (3 through the decay Bd J/ipK* . 

The lifetime of is F^ 1.5 ps It is conventionally measured through the 
self-tagging semileptonic decays, and in terms of the lifetimes B^ and B^ , is given by 



(Fi + n)TiTi, ^ 1 - i (f" 



2 

d 

d 



^ V J 

Since the difference between Td{sl) and = (F^ + F^)/2 is quadratic in the small 
quantity AF^/r^, we shall neglect it and shall take F^ = Td{sl) = Td in our analysis. 

The lifetime difference AF^ affects the measurement of (3 through Bd —>■ J/ipKs. 
The CP asymmetry measured in order to determine (3 through this mode is 
, _ T[Bd{t) J/i^K,] - T[Bd{t) ^ J/i,K,] 

T[Bd{t)^J/,ljKs]+T[Bd{t)^J/^ljK,] ^ ^ 

e-^''*sin(AMrft)sin(2/j) 
cos^ p e + sm pen 
where (3 is the experimentally measured value. In the limit ATd —>■ 0, the expression 
reduces to Acp = sin(AMdt) sin(2/3), which is the approximation used in the absence 
of accurate enough time measurements. The error introduced in the measurement of /3 
due to this aproximation is 

^j^^^ = cos^ f3 e-^^^'/' + sin' f3 e^'^'^'/' - 1 (103) 
sin(2/?) MP y , 

^ - cos(2/5)AFrft , (104) 
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so that for t ~ 1 /F^^, the error introduced due to neglecting the hfetime difference is of 
the order of ATd/Ta ~ 0.5%. 

At LHCb, the proper time resolution is expected to be Ar ^ 0.03 ps [Q, at least 
when the decay vertex lies within the silicon strips and its position can be accurately 
determined. This would imply that, if the number of relevant events with the proper 
time of decay measured with the precision Ar above is N , then the value of Ar^/r^ is 
measured with an accuracy of 2 ■ 10^^/-\/iV. With a sufficiently large number of events 
A^, it would be possible to reach the accuracy of 0.5 • 10^^ or better. 

The above naive argument works if there is a final state that decays with a given 
lifetime or F^. However, such a state does not exist (e.g. decays into CP eigenstates 
does not help since the — Bd mixing phase (2/5) is large and the CP eigenstates are 
far away from the lifetime eigenstates). So we have to find a way to separate a mixture 
of two exponential decays with very similar lifetimes. 

The time evolution of the decay rate of B^ into CP eigenstates is given by 

\A+{t)\^ = 1^(0) |2 \cos^(3 e-^'* + sm^(3 e"^^* + e-^-** sin( AM^t) sin(2/?)j , (105) 

\A-{t)\^ = \A-{0)\^ [sin^p e-^i' + cos^(3 e"^^* - e"^^* sin(AMdt) sin(2/3)] (106) 

for CP even and odd final states respectively. Separating the two lifetimes from the non- 
oscillating part of the time evolutions above is a formidable task. However, if instead 
we use a decay to two vector mesons {B VV), the angular distribution of the events 
can help us in separating the two lifetimes. It has been shown in that the additional 
information due to the measurement of the angles can be an order of magnitude more 
than the information in the decay time alone. 

Let us have a close look at the decay B^ — >■ J/ipK*{—^ KsTt^), which will be useful 
in the determination of AF^^. The most general amplitude for this decay is given in 
terms of the polarizations ej/^, ex* of the two vector mesons: 

A{B, J/ijK*) = A, efj^e*^. - ^ e}J^ ■ e^^, - e}/^ x e^,. ■ p (107) 

where E^* is the energy of the K* and p the unit vector in the direction of K* in the 

J/ip rest frame. The superscripts L and T represent the longitudinal and transverse 

components respectively. Since the direct CP violation in this mode is negligible, the 

amplitudes Aq and A\\ are CP even whereas A± is CP odd. Let us define the angles 

as follows. Let x axis be the direction of K* in the J/ip rest frame, and z axis be 

perpendicular to the decay plane of K* K^tt^, with positive y direction chosen such 

that Py{Ks) > 0. Then we define {6,ip) as the decay direction of in J/ijj rest frame 

and as the angle made by Ks with x axis in K* rest frame. 

The transversity angle distribution, which is sufficient to separate the CP odd and 

even final states, is given by 

dT\B ^ J/M+i-)K*(KsTT°)] 3,, . , 3,, . , , 

J " «-\A4mi + co.^e) + -lA4t)?.inH , (108) 

where |^+(t)P = |^o(^)P + \A\\{t)\^ and \A-{t)\^ = |Ax(t)P are the CP even and odd 
components respectively. 
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The complete solution to the time evolutions involves a maximal likelihood fit to the 
five parameters F^, F^, /3, AM^, |^_(0)/^_|_(0)|. The value of AM^ can also be taken 
from other experiments. This not only would give the values of F^ and F^ separately, 
but it will also determine the value of sin^ /3, so that /? is measured without the discrete 
ambiguity (3 n /2 — (3 present in its determination through the gold-plated decay 

More ways for the determination of AF^ as well as (3 are available through the 
angular distribution in the angles 9 and yj: 



+ 1^11 



dcosO dip dt 
— sin^ 9 sin^ ip) 



2(1 -sin^^cosV) 
2 sin^ 6^ - Im (A| sin 29 sin p] 



(109) 



The time evolution of the coefficients of the four angular terms is 



\M 



\Mt)\' 



lm{A;{t)A^{t)} 



X 



= l^o(0)p 

cos^ /3 e"^'* + sin^ p e'^^* + e"^^* sin(AMdt) sin(2/?) 
= |A||(0)px 

cos^ (3 e'^^t + sin2 f] g""^' 
= l^±(0)rx 
sin^ P e'^L^ jf- cos^ (3 e 
= \Am\\A^{0)\ X 



Ht + g-rdt sin(AMdt) sin(2/3) 



e"^-^* sin(AMrft) sin(2/3) 



e"^ {sin 6i cos(AMrft) - cos 6i sin(AMdt) cos(2/3)} 



e ^ 



cos(5i sin (2/9) 



(110) 



where (5i = Arg(Ajj(0)A_L(0)). Note that even before reaching the precision to able 
to separate F^ and F^, the above can already measure the value of sin(2/9) through 
lm{Al{t)A^{t)}. 

The discrete ambiguity (3 7i/2 — f3 would remain unresolved in the absence 
of the lifetime separation, since the sign of cos^i, and hence the sign of cos(2/?) is 
undetermined. One way to resolve this ambiguity by using Bg J/'4^'P angular 
distribution and ?7-spin symmetry is suggested in [0]. But the time evolution ( |110| ) above 
offers a way without having to take recourse to any other decay or flavour symmetry. 

The non-oscillating part of ( |110| ) is 



C=[\m{Al{t)A^{t)}] 



NO 



— pd f —rd f 

e H - e 



cos(5isin(2/?) , (111) 



which is also the non-oscillating part of the corresponding charge conjugate decay 
Bd J/iIjK*{—>- KgiT^). The sign of this quantity is the same as the sign of cos^i, 
since F^ > F^. This in turn establishes the sign of cos(2/5) through ( |110| ). Note 
that, in the absence of any Ba — Bd production asymmetry, the non-oscillating part of 
( |110|) is exactly the quantity measured if the initial B meson was not tagged. Then 
determination of the sign of C would need neither tagging nor time measurements. 
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The Bj, — Bfi production asymmetry may be measured through the asymmetries 

in the pairs of self-tagging decays (a) DfD^ vs. B^i D~ or (b) 

Bd — ^ J/ipK*{K~^7T~) vs. Bd J / ip K* {K^ -K^) . The asymmetry in both the pairs 

of decays is a combination of the production asymmetry and a very small amount of 

direct CP violation. However, though the production asymmetry is the same in both 

the pairs, the direct CP violation is expected to be smaller for the pair (a). This can 

give us a handle on the B^ — Bd production asymmetry at the hadronic machines. In 

the case of B factories that run at the v{4:S) resonance, the production asymmetry is 

absent and the sign of C can be cleanly measured. For further issues we refer the reader 

to a forthcoming paper [§] . 
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6. Trigger and Event Selection Studies 

6.1. Overview of the CDF B-physics trigger strategy 

M.S.Martin, Oxford University 



6.1.1. Trigger hardware The CDF Run2 Trigger is a significant improvement over 
its Runl counterpart. The luminosity expected in Run2 is expected to go up to 
L = 2 X 10^^cm~^s~^|^. This will require the trigger to be able to take the interaction 
rate (about the 7.6MHz bunch crossing rate) and convert this into the rate to tape 
(about 50 ev/s). 

The physics capabilities of the Run2 trigger are greatly increased over Runl. We 
will be able to examine high pt tracks at Level 1 (hardware). We will be able to 
reconstruct (in 2D) high impact parameter tracks at Level 2 (hardware), and we hope 
to be able to perform offline quality reconstruction at Level 3 (software). 

Level 1 At Level 1, we will seek to find physics objects based on a subset of detector 
information. These include high energy muons and electrons, but also high pt tracks 
using the Extremely Fast Tracker (XFT). The XFT takes a subsection of the main 
tracking chamber P| (Central Outer Tracker) data and applies a two stage algorithm 
to locate high pt tracks. The COT is a large open cell drift chamber consisting of 96 
layers which are subdivided into 8 superlayers. For each one of four layers of the COT, 
the first stage of the XFT algorithm applies a mask of possible hit configurations for 
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pt > 1.5GeV. Then the second stage of the algorithm then examines 1.25° angular bins 
to put the high pt hits in these layers together to form high pt track candidates. This 
trigger forms the basis for all of the hadronic B-Physics triggers at CDF Run2. 

Level 2 Level 2 uses programmable hardware processors!^ to perform more 
sophisticated cuts on the physics objects obtained from Level 1. It will have an input 
rate of about 40 kHz. A significant Level 2 improvement over Runl for B-Physics is 
the Secondary Vertex Trigger[^ (SVT) which uses the main silicon detector (SVX) to 
find high impact parameter tracks. 

The SVT operates on a subset of the data from the SVX and the COT. It takes 
the high pt track candidates from Level 1 and forms superstrips which are used in 
conjunction with SVX data in the pattern recognition stage. The digitised data from 
the SVX is fed into the Hit Finder which clusters the strips into hits. The combination 
of hits and superstrips are then fed in parallel to a number of Associative Memory 
chips which compare the data to a list of previously computed legitimate combinations. 
Finally, these legitimate combinations (roads) are sent to the Track Fitter Farm where 
they are fed into a linear approximation fit to yield 2D tracks with impact parameter 
information. 

Level 3 Level 3 consists of a fast switch serving a farm of Linux PC's. The output 
rate of Level 2 is going to be about SOOifz, and the rate to tape is going to be about 
50ef/s. The necessary rejection will be obtained with a combination of confirmation of 
Level 2 measurements and more sophisticated cuts. We will only track in the parts of 
the detector which are necessary since there are too many channels (eg in the silicon) 
to do global tracking. 

6.1.2. B-phy sics triggers There are two types of B-Physics trigger to be implemented 
in Run2, Leptonic and Hadronic. The following sections describe the specific trigger 
strategies we intend to pursue for B-Physics. The modes mentioned are by no means 
exhaustive, but are rather the modes more relevant for CKM element measurements 
and CP- Violation. 

Leptonic triggers Firstly, building on extensive experience in Runl we will be 
triggering on B-decays with leptons in the final state. The J/ip ^ trigger [§] will be 
the most powerful version of this trigger and will make use of Level 2 muon objects to 
form J/%1) candidates. We will also implement a, J/if) ee^ trigger. 

Secondly, we will implement a lepton plus displaced track trigger to select more 
generic leptonic B-decays. The lepton will come from the Level 2 muon and electron 
objects, and the displaced track will come from the SVT. 

Finally we are going to try to implement a Radiative B-Decay Conversion trigger [p!0| 
pn| , where we attempt to trigger on an electron from the conversion, and a nearby 
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displaced track from the B to bring the trigger cross-section down. This has a large 
overlap with the Lepton plus displaced track trigger. 

J/ip ^ /i/i triggers There are several interesting physics modes which are 
accessible via a. J/ip trigger. The two I will highlight here are the B — >■ J /if} Kg and 
B, ^ J/V^0. 

The baseline requirements at Level 2 (Level 1 mirrors Level 2) for a, J/ip candidate 
are two central muons with pt > 2.0GeV. A central muon in the trigger is defined 
as a stub in the CMU0 (Central Muon) or CMX (Central Muon Extension) which 
matches a track in the COT. The muons must also pass a transverse mass requirement 
L5 < rrirp < 3.25GeV. We are also investigating lowering the CMU pt requirement to 
1.5GeV to improve this trigger. 

At Level 3 we will require an invariant mass cut 2.85 < m < 3.25GeV. It is 
expected that the cross-section out of Level 3 for this trigger will be 20nb leading to 
40 million J/ip signal plus sideband events in the first 2fb^^ of Run2. The number of 
signal events we expect for 2fb^^ are given in table These estimates are obtained 
from scaling observed Runl samples of these two modes. 

Table 12. Flagship J/'0 mode sample estimates 



Mode 


Expected sample size 


B ^ J/^K, 


20,000 


Bs ^ J/i^(p 


4,000 



J/ip ^ ee triggers]]^ Although the J/ip —>■ fifi trigger is by far the most powerful, 
we will seek to use electrons to trigger on J/ ip candidates. 

At Level 1 we require 2 electron candidates with Et > 2GeV and pt > 2GeV. 
We also require Had/ Em < 0.125 where Had is the deposit of energy in the hadronic 
calorimeter, and Em is the same for the electromagnetic calorimeter. 

At Level 2 we make positional requirements on the shower using the Shower- 
max detector [|12| (a proportional wire chamber), and require opposite sign tracks with 
A0 < 90° (transverse angle). There is also a full invariant mass cut (in contrast to 
the transverse mass for muons) imposed 2.5 < m < 3.5GeV. This strategy could be 
changed to match the muons though. 

At Level 3, we will impose the same mass requirement as for the muon version. In 
addition we will require further cuts based upon soft lepton tagging [§]. 

After Level 1, the cross-section is 7 to 18 fib depending on the luminosity. This 
could be reduced by additionally imposing opposite charge and an angle requirement 
at Level 1. After Level 2, the cross-section (without the extra Level 1 cuts) is 100 nb, 
with negligible luminosity dependence. The Level 3 cross-section is estimated to be 6 
nb leading to an estimate 12 million of these events in the first 2fb^^ of Run2. 
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Lepton plus displaced track triggers |j^ Triggers of this type gain access 
primarily to semi-Ieptonic B-Decays. Tlierefore tliis trigger is relevant for collecting 
samples to constrain the CKM elements. This trigger subdivides into Muon, and 
Electron plus displaced track. 

For the muon path at level one, we require a central muon matched to an XFT 
track (the Level 1 track trigger) with pt > AGev. For the electron path, we require a 
AGeV Central Electromagnetic calorimeter tower with Had/ Em < 0.125. Also, as for 
the muon we require a 4GeV XFT track matched to the calorimeter tower. These Level 
1 triggers fully overlap with other triggers. 

At Level 2 for the muon path we confirm Level 1 and require one SVT track with 
120/im < \dO\ < 1mm which is correlated to the muon by opening angle and transverse 
mass: < 90° and < 5GeV. For the electron path at Level 2, we require the 
same correlation with the associated track, and some positional information about the 
shower from the shower-max wire chamber. 

For the muon path the cross-section after Level 2 is 39 nb while for the electron 
path it is 37 nb. The cross-section after Level 3 needs more study, but the estimated 
data set size is 60-134 million for 2fb^^. 



Conversion trigger ||1C)| | | |11|| The physics goal of this trigger is to access the 
B ^ s'y type of decay (for example B K*'y, or B 07). 

This trigger is very similar to the Lepton plus Displaced Track trigger, except that 
it doesn't have the correlation requirements with the associated track. This means the 
cross-section goes up, and so we may require a second SVT track at Level 2. 

Without the second SVT track required at Level 2, the cross-section out of Level 
2 is ~ 30nb. The Level 3 cross section is estimated to be 6nb, leading to a 12 million 
event data set for 2fb~^. 



Hadronic triggers ||13|| ||14|] ||15|] [|16|| In contrast to leptonic B-decays, this class of 
decay will be accessible solely via the SVT. The strategy is to perform a range of cuts 
on combinations of two tracks. The generic term for this is the Two Track Trigger. 
There are two flavours of this trigger: B vrvr, and Bg hadronic. 



B ^ mv The main physics goals of this trigger are to measure the CKM angles a 
and 7. The measurement of 7 will be attempted in conjunctions with another mode 0] 
which we hope to collect with this trigger, Bg K^K^ . We also require B J / ipKg 
for this measurement. We also hope to collect samples of -Bs/d nK with this trigger 
with a view to a direct CP- Violation Limit /Measurement. 

The cuts for this trigger are summarised in the table below. 

The A0 cuts are made in 2D on the transverse angles. The difference between this 
trigger and the next one is that the opening angle of the two tracks is constrained to be 
quite large whereas the opening angle in the Bs hadronic trigger can be more restricted. 
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Table 13. B ^ tttt Two Track Trigger (dO is impact parameter) 



Trigger level 


Cuts on each track 


Cut on the pair 


LI 


pt > 2GeV 


Spt > 5.5Gev, A(p < 135° 


L2 


100/i < do < 1mm 


20° < A0 < 135° 


L3 


Under discussion 





This effectively makes this trigger into a high Q- Value trigger, and correspondingly the 
subsequent one into a low Q-Value trigger. 

At Level 3 the base-line is verification of Level 2, which will require tracking 
regionally in the silicon as mentioned in section |6.1.1| . The cross-section of this trigger 
is considered with the subsequent one below. 

Bg hadronic The main purpose for this trigger is to collect events for a 
measurement of Bg mixing. However, many other fully reconstructed B decays will come 
in on this trigger. One example is Bg oi*^ oi*^ in which CP- Violation is predicted to 
be zero. 

The cuts for this trigger are summarised in the table below. 



Table 14. Bg hadronic Two Track Trigger (dO is impact parameter) 



Trigger level 


Cuts on each track 


Cut on the pair 


LI 


pt > 2GeV 


Spt > 5.5Gev, A(p < 135° 


L2 


100/i < do < 1mm 


2° < A0 < 90°, {pt ■X^)>0 


L3 


Under discussion 





The {pt ■ X^) > cut is defined as follows. Vj, is the vector pointing from the 
primary vertex to the secondary vertex, and so, the dot product pt ■ Xy will be positive 
for tracks originating in front of the primary, and negative for tracks originating behind 



it. See figure 3.1.2 




Dot product 
Negative 

Primary Vertex 



Figure 25. Diagram explaining the effect of {pt • X^,) > cut 



At Level 3, as for B ^ tttt we will at the very least verify Level 2, but we will 
probably need to do something more sophisticated than this to beat the cross-section 
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down to an acceptable level. 

The trigger rates are estimated to be (for the initial configuration of the accelerator) 
250n6 at levell, 560 ± 125n6 at Level 2, and 200nfe at Level 3. This would lead to a 200 
million event data set size. 
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6.2. A brief outline of the LHCb trigger 

C. Shepherd- Themistocleous, Cambridge University 



6.2.1. Introduction B meson decay processes of interest to CP violation physics have 
small branching ratios, typically of the order of 10~^. The production of large quantities 
of B mesons is therefore required for a study of these processes. Such large rates, of the 
order of 10^ bb events/s will be available at the LHC collider. The price that has to be 
paid for these large rates is a very high background rate. The rate of bunch crossings 
with at least 1 interaction at the LHCb interaction point will be of the order of 12MHz. 
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One of the greatest challenges of the experiment is therefore to be able to trigger on 
interesting events. The event rate is so high that it is not even possible to record all B 
events. This therefore means that a highly selective and efficient trigger is required. It 
also has the important consequence that the designers of the trigger should be aware of 
all potentially interesting final states while the trigger is being designed. If a particular 
channel is only found to be of interest at a later date then it is likely that the trigger 
will be inefficient for this particular channel. While the parts of the trigger that are 
algorithms running on processors may easily be changed, changing the hardware is a 
much more difficult task. The strategy that will be followed at the LHCb detector is 
briefly outlined below. 

6.2.2. Trigger strategy at the LHCb experiment |l| 

The beam crossing rate at the LHC is 40MHz and the LHCb trigger is designed to 
be ready to receive data at every beam crossing. This means that the detectors must be 
readout at a rate of 40MHz and a trigger decision must also be made a this rate. It is 
clear that a sophisticated trigger algorithm cannot be run in the 25ns available between 
beam crossings. The strategy followed therefore is to use a pipeline trigger and to use 
several different levels of trigger. In a pipeline trigger data are stored for a fixed time 
while processing is performed to make a trigger decision. The processing takes place 
in parallel for different events. Therefore, in the first level of the trigger for example, 
while a decision is made at a frequency of 40MHz, 4/is are available for each event to 
be processed. 

The Level trigger utilizes the fact that decay products of B decays have on average 
higher momenta transverse to the beam (pt) direction than background particles. This 
trigger looks for high pj muons, electrons and hadrons as well as photons and possibly 
neutral pions. These triggers take place in parallel and are then combined into a final 
Level decision. At this stage a veto trigger on beam crossings that contain more than 
1 interaction may also be used. This uses part of the vertex detector to identify multiple 
primary vertices. This level of the trigger is designed to produced an output of IMHz. 

The next trigger level utilizes the other major characteristic of B mesons, namely 
their long lifetimes. This leads to a mean decay length of 7mm at LHCb and hence 
to decay vertices clearly separated from the primary vertex. A silicon vertex detector 
placed close to the beam pipe is used to look for these secondary vertices. 

These first two levels of the trigger are "hardware" triggers in the sense that it 
will be hard to change these triggers once the detector is built. Any changes would 
have implications for aspects such as the physical configuration in which the detectors 
are read out and the speed of operation of specialized chips on electronics boards. The 
first two levels of the trigger have to operate a high rates and hence essentially use one 
subdetector in each trigger for speed. 

At this point all the data is made available to a processing farm and all subsequent 
triggers are software algorithms. Data from several subdetectors can now be combined 
to refine the earlier triggers. An example is the refinement of the vertex trigger. The 
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vertex detector is not in a magnetic field and a source of fake secondary vertices is 
from multiply scattered low momentum tracks. This background can be removed by 
using information from tracking chambers. Complete tracking information will also be 
available to selection algorithms in the trigger. After all generic B selection algorithms 
the event is rate still too high to allow the recording of all events selected at this point. 
Once this stage has been reached algorithms to select specific B decays or classes of B 
decay are performed. Data will finally be written to tape at rate of about 200Hz.|J| An 
important difference between this trigger and the trigger used in CDF is that analysis of 
data in the later stages of the trigger is not restricted to regions of interest in the detector 
that have been defined by earlier stages of trigger. The processing power available here 
is such that this time saving device is not required. 

The discussion about triggers that took place in the CP violation working group 
proved to be a very useful one. In particular many people where made aware of the 
difficulty of triggering in the hadron collider environment and of the importance of 
knowing early in the design phase what final states must be triggered on. Many useful 
discussions on what ought to be included in the LHCb trigger took place. 
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6.3. Separating T(4S) decays from continuum events using a neural network 

Fabrizio Salvatore and Glen Cowan, Royal Holloway 

6.3.1. Introduction Artificial Neural Networks (NN) have been applied to a variety 
of problems in high energy physics in order to discriminate between different classes of 
events. Here the technique is used to separate hb events from T(4S) decays ('signal') from 
the continuum background for e"'"e~ interactions at a centre-of-mass energy £^cm = 10.58 
GeV, i.e., at the T(4S) resonance. The analysis relies on the different topologies for 
signal and background at a B Factory: while hh events are more 'spherical', uds and 
cc events are more 'jet-like'. This is due to the fact that in the T(4S) rest frame, the 
momenta of the produced B mesons are small and they decay isotropically. On the 
other hand, events where light quarks (mm, dd or ss) are produced are characterized 
by a preferred direction with hadron jets following roughly the quark and antiquark. A 
similar jet-like structure is also present in cc events, although there it is less pronounced. 

For this study, two samples of 600k events of the type e"'"e~ qq have been 
generated using PYTHIA all at a centre-of-mass energy equal to the mass of the 
T(4S), where the quark flavour g is a mixture of u, d and s in the first set and c in the 
second. A third data sample of the same size consisted of T(4S) decays. The analysis 
has been performed at the level of the generated particles and no attempt has been 
made at this stage to simulate detector effects. 

ff This number is under revision and may well change in the near future. 
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6.3.2. The shape variables The variables used to discriminate jet-like from isotropic 
events are described below. Further information on the variables can be found in |^. 

(a) The ratio R2 = H2/H0 of the 2nd to the 0th Fox- Wolfram moments. Neglecting 
particle masses, 4-momentum conservation requires that Hq = 1. For a two-jet 
event, Hi = and ~ 1 for / even and if/ ~ for I odd. The R2 distribution is 



shown in Fig. where the different behaviour for signal and background events 
can be observed. 

(b) The thrust T = J2i 1^ ' Pi|/ J2i I Pi I; where T is the thrust axis of the event 0. For 
an isotropic event one has T = 0.5 while T ~ 1 for a highly directional one, as 
shown in Figs. (a) and (c). 





Figure 26. R2 distribution for (a) T(4S) signal and (b) continuum events. 



The variables defined above are quite powerful for identifying events where the two- 
jet topology is well defined and clearly distinguishable from the isotropic case. The signal 
to background separation decreases, however, when one or more gluons are emitted, 
which is more likely to happen when light quarks are produced. In this case there are 
three or more jets, so that the event shape resembles more that of the signal. 

To discriminate further between signal and background, the Durham clustering 
algorithm is applied to group particles into jets. For every pair of particles i and j, 
the 'distance' i/ij is computed as 

2 mm{Ef,E]) (l-cos%) 



(112) 



where Ei and Ej are the particles' energies, 9ij is their opening angle and i?vis is the total 
(visible) energy in the event. The pair with the smallest yij is replaced by a pseudo- 
particle with four-momentum = +Pj- The procedure is repeated until jets have 
been found. In our case iV = 3 or 4 and the smallest yij values are called and 7/4, 
respectively. Once the clustering has been done, the following discriminating variables 
can be added to the list: 
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(c) Z/3, Va] 

(d) 6'bz, the angle between the highest and second highest energy jet; 

(e) 6'ks5 the angle between the lowest and second lowest energy jet; 

(f) ^NR, the angle between the planes defined by the highest and lowest energy jets 
and by the other two jets; 

(g) ^34) the angle between the planes defined by the highest and second highest jets 
and by the other two jets; 

(h) the QCD four-jet matrix element squared, \Mqcd\'^i averaged over all possible 
assignments of the jets to the final-state partons |^. 

As an example, the distribution of In \Mqcd\'^ is shown in Fig. |2^(b) and (d) for 
signal and background events. 




0,5 0,6 0,7 0,8 0,9 1 3 4 5 6 7 



T InlMI^ 

Figure 27. Distribution of (a) thrust and (b) In \ Mqcd\'^ (right) for signal events; (c) 
and (d) show the same for continuum events. 

6.3.3. The analysis results All the variables defined in the previous section have 
been combined in a Neural Network using the package JETNET 3.51 0. The network 
is a 3 layered feed- forward one, with 9 input nodes, 12 hidden nodes and one output. 
The network has been trained using 100k 66, 100k uds and 100k cc events. Each input 
variable has been scaled linearly so that it lies in the range [-1,1]; in this way at the 
start of the training no single variable dominates the inputs to any neuron. The training 
algorithm is 'error back-propagation'. 

After the training, the NN has been used with an independent sample of events, 
500k for each type, and the output distributions are shown in Fig. The signal 
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0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 

NN output NN output 




NN output 

Figure 28. Distributions of the NN output for (a) signal and (b) background events. 
The distributions are superimposed in (c) after being normalized to unity. 



and background curves are very well separated, with continuum background events 
accumulating around and signal events around 1. 

To evaluate the rejection power of the NN we can set a cut on the output such that 
the efficiency for the signal is 90%; this corresponds to requiring an output greater than 
0.35. The fraction of background events surviving the cut is 14% for uds events and 
18% for cc events. 

6.3.4- Conclusions A feed-forward neural network has been applied to separate qq 
continuum events from T(4S) decays produced in simulated e^e^ interactions at the 
T{4S) resonance. Preliminary results are quite encouraging. For a signal efficiency of 
90%, only ~ 16% of background events survive the NN cut. More work needs to be done 
to verify that the results are still valid after a full detector simulation; nevertheless, the 
NN technique seems a promising option for this kind of analysis. 
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7. QCD Factorization 

7.1. Factorization in charmless B decays 
W.N. Cottingham 

Apart from electro-weak terms and other small corrections and taking account 
of QCD corrections the effective Hamiltonian for b quark decay is a sum of current 
current terms j ■ j. Factorization approximates the quasi two body sum of matrix 
elements (/ii/i2 | j ■ j \ B) by a sum over products of single body matrix elements 
{hi I I B){h2 I I 0). So that the transition amplitude ||, |^ becomes 

{hih2 I H,ff I B) = GF/V2J:fP,[Qi{hi, h2) + Qi{h2, h^)] 



Qiih 


hi) 


= {hi wy'^ 


[1 
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s 







The coefficients Pi are constructed from Wilson coefficients 0, ^] and CKM 
parameters. The Wilson coefficients have been quite reliably estimated from QCD. 
The single particle matrix elements, for example (vr"^ | 'U7'^(l — 75)^ | 0) = fnpf^, involve 
soft QCD parameters, in this case /tt and many are quite precisely known for example 
from measured leptonic decay rates. The matrix elements {hi \ j'^ \ B) are in principle 
measurable from B semi leptonic decays, they involve only a few parameters. For the 
moment they are only available as theoretical estimates. 

The physical justification for the factorisation approximation is the idea of colour 
transparency. It is argued that the b quark by virtue of its large mass imparts a large 
velocity to the light quark decay products, the pair which form the light particle hi 
move away but stay together as a colour singlet so escape from the gluonic environment 
of the B meson without interaction as does a lepton pair in the semi leptonic decay 
B ^ hi + I + V which is described by a matrix element {hi \ \ B). This rationale 
has recently been given theoretical justification [Q, it has been claimed that for decays 
into a pair of light mesons and in the heavy quark limit that the corrections to the 
factorisation approximation can be calculated from first principles of perturbative QCD. 
The corrections involve, as you might expect, hard gluon interactions with the spectator 
quark but the only long distance physics, other than the current matrix elements, 
involves the meson light cone distribution amplitudes. The beauty of the results 
presented in for tttt decay is that the coefficients Pi acquire additional perturbative 
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corrections. The difference between the decay rate formulae of the factorisation 
approximation and the formulae with the QCD corrections lies only in the coefficients 
Pj, the factorisation matrix elements are common to both. Generalisations to other 
decay channels can be found in |Q ^, |^. 

The corrections to the coefficients Pi so far presented are to first order in a^. An 
encouraging feature is that they are not large, which leads one to hope that the precision 
with which the standard model can be tested will be determined by the proximity of the b 
quark mass to the heavy quark limit and the precision of our knowledge of the soft QCD 
parameters, B meson semileptonic transition form factors, meson light cone distribution 
amplitudes etc. Because of the involvement of soft QCD parameters, some of which are 
only poorly known, and to be confident of the conclusions about the standard model, it 
will be important to have a consistent picture of as many channels of charmless B decay 
as possible. 

Limiting the discussion to the light hadrons hi and /i2 being the lowest mass pseudo 
scaler and vector mesons there are only a few poorly known soft QCD parameters. 
Because the B has zero spin a spin one meson with a spin zero partner is confined to 
the helicity zero state. Also in the heavy quark limit and to lowest order in (m/j/ms)^ 
only helicity zero states should be involved in transitions to two spin one mesons . 
The poorly known B transition form factors should then involve only five parameters 
that describe the transitions B tc, B p, B ^ uj, B ^ K, B K* Setting aside 
transitions involving the 7] and i]' there are about forty channels of B^ and B~^ decay 
with branching ratios anticipated to be greater than 10~^ [0]. It can be expected that 
the basic validity of the improved factorisation approximation will be tested within the 
next year or two, and the proximity of the b quark mass to the heavy quark limit 
will determine the sharpness of focus that this method provides on the standard model 
description of the weak interaction. 
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7.2. Phenomenological impact of the QCD improved factorization approach 

Content identical to JJ/ 

M. Beneke, RWTH Aachen 



7.2.1. Introduction The observation of B decays into nK and tttt final states has 
resulted in a large amount of theoretical and phenomenological work that attempts to 



CP Violation, UK Phenomenology Workshop 2000 



86 



interpret these observations in terms of the factorization approximation (FA), or in terms 
of general parameterizations of the decay amphtudes. A detailed understanding of these 
amplitudes would help us to pin down the value of the CKM angle 7 using only data 
on CP-averaged branching fractions. Theoretical work on the heavy-quark limit has 
justified the FA as a useful starting point 0, but predicts important and computable 
corrections. Here we discuss the most important consequences of this approach for the 
T{K and tttt final states. 

To leading order in an expansion in powers of Aqcd/'"^;), the B txK matrix 
elements obey the factorization formula 

+ /+^''(0)/.T^*<|.. (113) 

where Qi is an operator in the weak effective Hamiltonian, f^^^\Q) are semi-leptonic 
form factors of a vector current evaluated at = 0, are leading- twist light-cone 
distribution amplitudes, and the *-products imply an integration over the light-cone 
momentum fractions of the constituent quarks inside the mesons. When the hard- 
scattering functions T are evaluated to order Eq. ( |113D reduces to the conventional 
FA. The subsequent results are based on kernels including all corrections of order as- A 
detailed justification of ( |113| ) is given in Ref. p. Compared to our previous discussion 
of TTTT final states the present analysis incorporates three new ingredients: 

i) the matrix elements of electroweak (EW) penguin operators (for t:K modes); 

ii) hard-scattering kernels for general, asymmetric light-cone distributions; 

iii) the complete set of "chirally enhanced" l/ruh corrections 0. 

The second and third items have not been considered in other 0] generalizations of 
Ref. to the nK final states. The third one, in particular, is essential for estimating 
some of the theoretical uncertainties of the approach. 

We now briefly present the input to our calculations. Following Ref. |Q , we obtained 
the coefficients ai{7[K) of the effective factorized transition operator defined analogously 
to the case of tttt final states, but augmented by coefficients a-/_iQ{7rK) related to 
EW penguin operators and electro-magnetic penguin contractions of current-current 
and QCD penguin operators. A sensible implementation of QCD corrections to EW 
penguin matrix elements implies that one departs from the usual renormalization-group 
counting, in which the initial condition for EW penguin coefficients is treated as a 
next-to-leading order (NLO) effect. Our NLO initial condition hence includes the 
corrections computed in Ref. 0. 

Chirally enhanced corrections arise from twist-3 two-particle light-cone distribution 
amplitudes, whose normalization involves the quark condensate. The relevant 
parameter, 2fiT,/mh = —4:{qq)/{f^mb), is formally of order Aqcri/^b, but large 
numerically. The coefficients ag and as are multiplied by this parameter. There are 
also additional chirally enhanced corrections to the spectator-interaction term in (|113|) , 
which turn out to be the more important effect. In both cases, these corrections involve 
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logarithmically divergent integrals, which violate factorization. For instance, for matrix 
elements oi V — A operators the hard spectator interaction is now proportional to 
{u = 1 — u) 

^ dudv ^ , . f ^ , . 2u7r u 



Iq u V \ rrih 

when the spectator quark goes to the pion. (Here we used that the twist-3 distribution 
amplitudes can be taken to be the asymptotic ones when one neglects twist-3 corrections 
without the chiral enhancement.) The divergence of the u-integral in the second term 
as V —>■ implies that it is dominated by soft gluon exchange between the spectator 
quark and the quarks that form the kaon. We therefore treat the divergent integral 
X = ^^[dv/v) as an unknown parameter (different for the penguin and hard scattering 
contributions), which may in principle be complex owing to soft rescattering in higher 
orders. In our numerical analysis we set X = ln(mB/0.35 GeV) + r, where r is 
chosen randomly inside a circle in the complex plane of radius 3 ("realistic") or 6 
("conservative"). Our results also depend on the 5-meson parameter]^ Ab, which we 
vary between 0.2 and 0.5 GeV. Finally, there is in some cases a non-negligible dependence 
of the coefficients ai{TrK) on the renormalization scale, which we vary between mb/2 
and 2mb- 

7.2.2. Results We take = 0.085 and m,(2GeV) = 110 MeV as fixed 

input, noting that ultimately the ratio iKifc/Kfe], along with the CP- violating phase 
7 = arg(V^*^), might be extracted from a simultaneous fit to the B ttK and B tctt 
decay rates. 

SU{3) breaking Bounds[^ |^ on the CKM angle 7 derived from ratios of nK 
branching fractions, as well as the determination of 7 using the method of Ref . |^ , rely 
on an estimate of SU{3) fiavour-symmetry violations. We find that "non-factorizable" 
5'f/(3)-breaking effects (i.e., effects not accounted for by the different decay constants 
and form factors of pions and kaons in the conventional FA) do not exceed a few percent 
at leading power. 

Amplitude parameters The approach discussed here allows us to obtain the decay 
amplitudes for the vrvr and ttK final states in terms of the form factors and the light-cone 
distribution amplitudes. The vr^vr^ final state is very poorly predicted and will not be 
discussed here. We write 

AiB'^n+n-)=T[e'^ + iP/TU] 

and parametrize the ttK amplitudes by0] 

A{B^ -> 7r+i^°) =P{l-ea e'^^e'^) , 

-V2A{B+ 7r°ir+) = p\l-ea eV^ - £3/2 e'^ie'^ - q e'")] , (114) 



-A{B^ n-K+) = P 



1 - Ea e'^^e'^ - St e'^^{e'^ - qc e*"^) 
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Range, NLO 


LO 




(0.017-0.020) e^[i=^'2i]° 


0.02 


£3/2 e'^ 


(0.20-0.38) e*[-^°'^l 


0.36 




(0.53-0.63) e*[-7'3r 


0.64 


St e'^T 


(0.20-0.29) e^[-i9'3l° 


0.33 


qc e^'^c 


(0.00-0.22) [-180,180]° 


0.06 


{P/TU 


(0.19-0.29) [-1-231° 


0.16 



Table 15. Parameters for the B ttK amplitudes as defined in (lf4), for conservative 
variation of all input parameters (see text). 



and y/2A{B^ n'^K'^) = A{B+ tt+J^O) + ^/2A{B+ n'^K+) - A{B'^ tx-K+). 
Table ^ summarizes the numerical values for the amplitude parameters for the 
conservative variation of X, and variation of the other parameters as explained above. 
The LO results correspond to the conventional FA at the fixed scale = mi,. They 
are strongly scale dependent. In comparison, the scale-dependence of the NLO result 
is small, with the exception of qc e^^'^ ■ One must keep in mind that the ranges may 
overestimate the true uncertainty, since the parameter X may ultimately be constrained 
from a subset of branching fractions. This is true in particular for the quantity £3/2 in 



Table 15, which can be extracted from data. 



Ratios of CP-averaged rates Since the form factor /+(0) is not well known, we 
consider here only ratios of CP-averaged branching ratios, discarding the tt'^tt^ final 



state. We display these as functions of the CKM angle 7 in Fig. [29 . 

Table |15| shows that the corrections with respect to the conventional FA are 
significant (and important to reduce the renormalization-scale dependence). Despite 
this fact, the qualitative pattern that emerges for the set of tiK and vrvr decay modes is 
similar to that in conventional factorization. In particular, the penguin-tree interference 
is constructive (destructive) m B ^ ti^ti~ [B ti^K^) decays if 7 < 90°. Taking the 
currently favoured range 7 = (60 ± 20)°, we find the following robust predictions: 



Br(7r+7r-) 
Br(7r=Fii'±) 

2Br(7rOis:o) 

Br(7r±iro) 
TB+ Biin^K^] 



0.5-1.9 [0.25 ±0.10] 

0.9-1.4 [0.59 ± 0.27] 

0.9-1.3 [1.27 ±0.47] 

0.6-1.0 [1.00 ±0.30] 



tbo Br(7r±K0) 

The first ratio is in striking disagreement with current CLEO data[[T^ (square brackets). 



The near equality of the second and third ratios is a result of isospin symmetry. []^] We find 
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Figure 29. Ratios of CP-averaged B ttK and tttt decay rates. The scattered 
points cover a realistic (dark) and conservative (light) variation of input parameters. 
The dashed curve is the LO result, corresponding to conventional factorization. 



Br(5 ^ n^K°) = (4.5 ± 2.5) x lO^^ (Kf,/0.039)2(/^^^(0)/0.3)2 almost independently 
of 7. This is three time smaller than the central value reported by CLEO. 

CP asymmetry in B tt+tt" decay The stability of the prediction for the 
TT+TT^ amplitude suggests that the CKM angle a can be extracted from the time- 
dependent mixing-induced CP asymmetry in this decay mode, without using isospin 
analysis. Fig. |3y displays the coefficient S of — sin( AM^^t) as a function of sin(2Q;) for 
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Figure 30. Mixing-induced CP asymmetry in B ^ tt+tt^ decays. The lower band 
refers to values 45° < a < 135°, the upper one to a < 45° (right) or a > 135° (left). 
We assume a, /3, 7 e [0, tt] . 



sin(2/5) = 0.75, which may be compared with the resuh in Ref. ||^. For some values 
of S there is a two-fold ambiguity (assuming all angles are between 0° and 180°). A 
consistency check of the approach could be obtained, in principle, from the coefficient 
of the cos(AmB^t) term. 



7.2.3. Conclusions We have examined some of the consequences of the QCD 
factorization approach to B decays into irK and tttt final states, leaving a detailed 
discussion to a subsequent publication. Here we have focused on robust predictions 
for ratios of CP-averaged decay rates. Our result for the ratio of the B —>■ tt+tt" and 
B 71^ decay rates is in disagreement with the current experimental value, unless 
the weak phase 7 were significantly larger than 90°. 
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8. Statistical Errors 



8.1. Feldman and Cousins for Pedestrians 

R. Barlow, University of Manchester 

"What confidence is this wherein thou trustest?" 

Isaiah Chap. 36 v. 4 



8.1.1. What is a confidence level? Confidence levels are the subject of much 
confusion and argument. To bring out what they are, and what they are not, consider 
two statements which might validly be made with 90% confidence. (For illustration this 
note uses 90% confidence levels throughout, but of course any value can be chosen.) 

A grocer: "Our potatoes weigh between 100 and 400(7." 

A physicist: "The Higgs Boson has a mass between 100 and iOOGeV / c'^ . 

The superficial similarity conceals a difference. The grocer has weighed many 
potatoes and found that 90% of them have weights in the limits given. The statement has 
a 90% probability of being true, with probability defined in the conventional (frequentist) 
way as 

j\Tobjects 

P= lim '"^""f (1). 

Total 

The physicist has not weighed any Higgs Bosons. If they had, they would discover 
either that all the bosons are within the range, or that all the bosons are outside it. 
Equation (1) will give a probability which is either or 1. What they mean by a 90% 
confidence level statement is not 

"A Higgs boson mass has a 90% probability of lying between 100 and AOOGeV/c"^" 

but 

"The statement 'The Higgs Boson Mass lies between 100 and AOOGeV/c"^^ has a 
90% probability of being true." 

If this physicist goes through life making statements of this type, they will be correct 
in 9 cases out of 10. The 'probability' of 90% is still given as the limit of a frequency, 
but a rather different one. 

T\T statements 

P = lim (2) 

M ^ AT statements v 

As a technical detail, a confidence level of 90% actually means that the probability is 
at least 90% (this covers cases in which it is not calculable exactly.) So they are right 
at least 9 times out of 10. 



8.1.2. Frequentist confidence levels The basic tool for constructing confidence 
intervals is the Confidence Diagrarn^. Suppose that a parameter of interest a and an 
observed quantity x are related by a pdf P{x] a). (For example: x could be a number 
of events, and a a branching ratio.) 

Choose a value of the probability (e.g. 68%, 90%, 95%) and a strategy (e.g. one 
sided upper, one sided lower, two-sided symmetrical). For any value of a, we can make 
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statements about the probability of x lying in certain regions, and determine a range 
] within which we say (at the desired level of confidence) that x will lie. 




Figure 31. A two-sided confidence band. 



We then plot x^ and x+ on the diagram as a function of a. Suppose we chose to 
quote 90% central confidence limits. Any horizontal line across the plot represents a 
particular value of a and thus a pdf with its limits, as shown in two particular cases. 
There is a 5% chance of a measurement lying below the lower limit, and a 5% chance of 
it lying above the upper limit. The unshaded region between the limits is the confidence 
hand (or heW) . 

The next step is a neat bit of logic. For any value of a, the probability that x lies 
in the belt is, by construction, 90%. So if you make an observation x the probability of 
the pair (x, a) lying within the belt is 90%. A vertical line at the observed x gives the 
limits for a where it intersects the edges of the confidence belt. 

A similar technique gives one sided (upper) limits: for each value of a the value 
of X is found for which the probability of a result this small or smaller is only 10% 
(or whatever), and for a measured x you read off the upper limit a+ for which the 
probability of such a low result is only 10% - for a higher value of a the probability is 
even smaller. 




Figure 32. A one-sided confidence band. 
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8.1.3. Poisson statistics This theory of confidence hmits apphes to any distribution 
function, but in B physics its main use is in deahng with rare decays, where observed 
numbers of events are governed by Poisson statistics. This means that the variable x is 
discrete (integer) - denote it by A^. Confidence diagrams have a continuous ordinate but 
a discrete abscissa and the smooth curves become staircases. Integrations are replaced 
by summations. For an observed number N the limits are given by 

N N-1 

^P(r;a+) = 0.05 ^ P(r; a_) = 0.95 (3) 



This means that in general it won't be possible to find a range which satisfies ^ P{r; a) = 
0.05. One has to use the conservative requirement ^ P(r; a) < 0.05. 

a is a branching ratio, lifetime, or similar quantity describing a source of events. It 
trivially translates into a number of events 5*. The actual expected number is then often 
given by S + B, where B is some background which we will take as being reliably and 
accurately estimated. The probability of observing r events is the Poisson expression 

e~(s+^){S + By/r\ (4) 



8.1.4. Problem: non-physical regions Often this works perfectly well, but there 
are problems if the parameter estimates lie near some physical boundary. (For example 
Branching Ratios cannot be negative. Another relevant example is the interpretation 
of the limit on the Higgs mass from precision electroweak measurements when much of 
the range is already ruled out by direct searches.) 

For example, if B is 5.1, and you observe 1 event, what can you conclude about 
S"? Following the above procedure, you calculate 0.05 = P{r; 4.74) = P{r; 0.05) 
so S -\- B lies between 0.05 and 4.74, so you can say with 90% confidence that S lies 
between —5.05 and —0.36. One would have to be an extreme purist to publish such a 
statement. It is technically correct: you expect one in ten of your 90% CL statements 
to be wrong. In this case you just have independent confirmation that this is one of 
them; you happen to have had a large downward fluctuation in the background. There 
is no way of including the knowledge that S must be non-negative. 

8.1.5. Bayesian confidence levels Bayesian statistics^ handles these limits without 
problems using Bayes' theorem. 

P{a-x) = ^^P{a) 

here P{x\ a) is the usual pdf. P{x) is a prior distribution for x, usually ignored as being 
taken care of by normalisation. P{a) is the prior probability for a. A full discussion of 
Bayesian versus Frequentist methodologies can be found in 

The probabilities P{a) and P{a]x) do not have a frequentist definition (except as 
a delta function) but can be defined as subjective probabilities. The form of P{a) is 
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often taken as uniform, under the 'principle of ignorance'; this title is disingenuous as 
the assumption of a prior uniform in, say a lifetime r will lead to different limits than 
those obtained from a prior uniform in, say, the width h/r. According to Jeffreys' 
prescription one should take a uniform distribution for a location parameter, whereas 
for a scale parameter A the distribution should be proportional to ^ (i.e. uniform in 
In A) 0. Jeffreys' criteria [0 suggests that for a Poisson source S the appropriate prior 
should be taken as proportional to 

With a step-function prior (zero for negative values and uniformly constant for 
positive ones), the requirement for a one-sided confidence limit S on the signal at the 
level a from observed events with expected background B 

a= f P{S'\N)dS' 
Jo 

produces (after repeated integration by parts) the requirement 



^N 

Ei' P{r; B) 



l_„ = nZ(^L£±^ (5) 



which is found in early versions of the PDG For large iV, the denominator approaches 
the standard value of 1. For small N it expresses the fact that we have manifestly got 
a downward fluctuation of the background, and one should temper the requirement on 
the signal accordingly. 

8.1.6. What Feldman and Cousins do Feldman and Cousins [0 get us out of this 
dilemma by pointing out another one. 

In real life, physicists (being human) do not decide beforehand whether they will 
quote a one-sided or a two-sided limit. If the observations are high, a two-tailed limit 
will be appropriate; if they are low then a single upper limit is the best information that 
can be given. 

For example: given a background of 1.1, if you observed 10 events you would very 
possibly quote as a result the range of 4.3 to 15.1 (at 90% confidence). If you observed 
2 events you would probably say (with 90% confidence) that any signal was less than 
4.2. There is some intermediate value at which you would 'fiip-fiop' between the two 
states. 

This may be sensible, but it destroys the confidence belt logic. For small x values 



we have the one sided belt of Figure ^; for large values the two-sided belt of Figure ^ . 
The resulting hybrid (Figure ^3[) violates the basis on which the curves were constructed, 
as for the intermediate values of a the probability content of the belt is not the specified 
90%. 

Feldman and Cousins thus argue that the strategy for deciding whether to quote 
a one-sided or two-sided result must be decided beforehand, and decided in such a way 
that for all values of a the probability content of the belt is 90%. They propose a 'unified 
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Figure 33. The effect of 'flip-flop'. 



approach' in which a single algorithm decides whether to quote one-sided or two-sided 
limits, and gives the value (s) in question. 

There are an infinite number of solutions for limits xi and X2 satisfying 



f 



P{x; a) dx = 0.90 



(For discrete processes Ylll Pi^'^ ^ 0.90) 

One initially plausible algorithm is to include the region where the pdf itself is 
greatest, adjusting C to satisfy 



P{x; a) dx = 0.90 



P{x;a)>C 



(6) 



This in fact gives the most compact region. 

But consider an outcome with 1 events seen, with an expected signal and 
background of and 15.1 respectively. This is very unhkely (4 10~® probabihty) and 
this algorithm is bound to exclude the possibility. But given one observed event and a 
background of 15.1, any probabilities are going to be pretty small, and we would want 
to include signal in the band, on the grounds that for any other values the probability 
is even lower! It would be fairer to include points by comparing their probability to, 
say, the best that can be achieved. 



P{x;a) = 0.90 



P{x;a)/P{x;abe3t)>C 



(7) 



For Poisson processes with a known background B and unknown signal S, Sf,est is 
X — B, a the observed number is bigger than the expected background, and if it is 
smaller. 

A Feldman- Cousins confidence interval for rare decays is thus constructed as 
previously, with a 90% (or whatever) confidence band being constructed for each value 
of the source strength S. For each S the Poisson probability is calculated for each 
number of observed events r = 0, 1,2. . ., and its ratio taken to the 'best' probability 
(with S = r — BoTS = Oa.s appropriate), r values are included in the confidence belt 
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N 


B 


Conventional 






Uniform Bayesian 




Jeffreys' Bayesian 


Feldman- 






1 Tailed 


2 tailed 


1 Tailed 


2 Tailed 


1 Tailed 


2 Tailed 


Cousins 


10 


0.0 


< 15.4 


5.43 . . 


. 17.0 


< 


15.4 


6.17 . . 


. 17.0 


< 


14.8 


5.80 . . 


. 16.3 


5.50 


. . 16.5 




1.0 


< 14.4 


4.43 . . 


. 16.0 


< 


14.4 


5.1/ . . 


. 16.0 


< 


13.7 


4.73 . . 


. 15.3 


4.50 


. . 15.5 




5.0 


< 10.4 


0.425 . . 


. 12.0 


< 


10.4 


1.43 . . 


. 12.0 


< 


9.24 


0.345 . . 


. 10.8 


1.19 


. . 11.5 




10. 


< 5.41 


4.0 / . 


. . o.yo 


< 


O.Do 


U.zoo . 


. o.Uo 


< 


4. i o 


U.UIUO . 


. . D.iO 


< 


D.OU 




15. 


< 0.407 


y.o / . 


. . i.yo 


< 


A QQ 

4.cSo 


n 1 o/i 
U.lz4 . 




< 


o.U ( 


U.UU4ycS . 


. . 4.zi 


< 


Z.DO 


5 


0.0 


< 9.27 


1.97 . . 


. 10.5 


< 


9.27 


2.61 . . 


. 10.5 


< 


8.64 


2.29 . . 


. 9.84 


1.84 


. . 9.99 




0.50 


< 8.77 


1.47 . . 


. 10.0 


< 


0.77 


2.11 . . 


. 10.0 


< 


8.07 


1 '71 

1.71 . . 


. 9.27 


1.53 


. . 9.49 




2.5 


< 6.77 


— 0.530 


. . 8.01 


< 


6.85 


0.558 . 


. 8.09 


< 


5.66 


0.0580 . 


. . 6.87 


< 


7.49 




4.0 


< 5.27 


— 2.03 . 


. . 6.51 


< 


5.72 


0.244 . 


. 6.92 


< 


4.25 


0.0127 . 


. . 5.42 


< 


5.99 




10. 


< — 0.725 


—8.03 . 


. 0.513 


< 


O. 1 I 


u.uyu ( . 




< 


Z.OV 


U.UUoOO . 


Q 1 O 


< 


1.83 


1 


0.0 


< 3.89 


0.0513 . 


. . 4.74 


< 


o.oy 




A 1 A 
.4.(4 


< 


Q 1 Q 

o. lo 


U. i ( D . . 


. o.yi 


0.104 


. . . 4.36 




0.50 


< 3.39 


—0.449 


. . 4.24 


< 


3.51 


0.141 . 


. 4.36 


< 


2.48 


U.UU/OZ . 


. . 3.25 


< 


3.86 




1.0 


< 2.89 


—0.949 


. . 3.74 


< 


Q 0'7 
O.Z ( 


U.IUU . 


y1 1 1 


< 


O 1 Q 

z.lcS 


n r\r\A ao 
U.UU44Z . 


O OQ 

. . z.yo 


< 


3.36 




5.0 


< — 1.11 


—4.95 . . . 


— 0.256 


< 




u.Uoio . 


Q /I K 
. . 0.40 


< 


1 an 


n nnoQ'7 


O OK 

. . Z.Zo 


< 


1.20 




10. 


< — 6.11 


—9.95 . . 


— 5.26 




2.51 


0.0564 


3.25 




1.48 


n 0091 7 


. 2.10 


< 0.672 





0.0 


<^ z.ou 


None 


< 


2.30 


0.0513 . 


. . 3.00 


< 


1.35 


0.00196 . 


. . 1.92 


< 2.43 




0.10 


< 2.20 


None 


< 


2.30 


0.0513 . 


. . 3.00 


< 


1.35 


0.00196 . 


. . 1.92 


< 


2.33 




1.0 


< 1.30 


None 


< 


2.30 


0.0513 . 


. . 3.00 


< 


1.35 


0.00196 . 


. . 1.92 


< 


1.61 




5.0 


< —2.70 


None 


< 


2.30 


0.0513 . 


. . 3.00 


< 


1.35 


0.00196 . 


. . 1.92 


< 


0.770 


9 


3.7 


< 10.5 


0.995 . . 


. 12.0 


< 


10.5 


1.82 . . 


. 12.0 


< 


9.49 


0.775 . . 


. 11.0 


1.49 


. . 11.6 


9 


3.7 


< 12.0 


0.415 . . 


. 13.4 


< 


12.0 


1.24 . . 


. 13.4 


< 


11.0 


0.306 . . 


. 12.3 


0.994 


. . . 13.1 


9 


3.7 


< 7.49 


2.37. . 


. 9.39 


< 


7.51 


3.25 . . 


.9.40 


< 


6.51 


2.22 . . 


. 8.39 


2.64 


. . 9.09 



Table 16. Some examples of 90% confidence limits for various values 



in order of their ranking in this ratio until the total probability condition is fulfilled. 
The result is a diagram where low values of will just give one limit whereas higher 
ones give a range, and the probability content of the belt is valid. 

There is a loss of universahty as such confidence diagrams, unlike those of the 
previous scheme, are non-trivially different for different background values. 

8.1.7. Some examples 90% confidence intervals are shown for various values for an 
observed signal and expected background. Both 1-tailed and 2-tailed limits are given for 
conventional frequentist confidence intervals, and for Bayesian intervals, and Bayesian 
intervals are shown both for a step- function prior (uniform for positive values) and a prior 
according to Jeffreys' prescription, proportional to (Taking a prior proportional 

to 1/s gives divergent integrals.) 

It can be seen that where the signal is safely larger than the background, all 
methods are in broad agreement. Differences are a salutary reminder of the inevitable 
arbitrariness of any analysis. 
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As signal sinks below background, the conventional limits become manifestly 
erroneous. The Bayesian limits appear more sensible, though the differences between 
the two versions become larger. The Feldman-Cousins value (s) flip nicely. 

8.1.8. Some objections The new technique has been warmly welcomed by some, but 
has been criticised by others. 

In the 1-sided regime the limit decreases with the value of the expected background, 
as can be seen from the table. This is apparently counterintuitive: for 2 experiments 
with the same number of events seen the one with worse (i.e. higher) background has the 
better (i.e. lower) limit! The counterargument is that the inferior experiment has to be 
very lucky to see so few events. On average, an experiment with a higher background 
will see more events and thus give a higher limit: some experiments will always be 
more fortunate than others, and by increasing B for fixed you are increasing the 
experimental luck. 

The method does not cope cleanly with the case of zero observed events. This is 
always going to be a special case because, uniquely, if you observe events you know 
the number of signal events (zero), and all the background information is irrelevant. 
Modifications of the method have been suggested to get over this problem |Q. 

The method can flip to the two sided form in cases the physicist knows are 
inappropriate. When counting decays forbidden by the standard model, or even by 
basic conservation laws, the algorithm will give a two sided limit in cases (like 10 events 
on a background of 5) which are obviously just upward fluctuations of background. Here 
the experimenter is probably overstretching the applicability of the method. For more 
plausible decay channels the same experimenter would probably be happy to report a 
limit from the same data (10 on 5). The method, by its nature, knows only about the 
numbers and not about the physics. 

Broadly speaking (and with some exceptions), frequentists accept the Feldman 
Cousins method as an improvement on the original form, though it does still have some 
drawbacks. Bayesians dislike it, perhaps because it expands the scope of frequentist 
analysis in an area where previously Bayesian technique had a monopoly. 

8.1.9. Some questions Can this method be adapted to the setting of limits on the 
ratio of two numbers, as was done by James and Roos for the conventional frequentist 
confidence level? 

When one is tuning cuts in an analysis, one aims to optimise signal/background or 
signal/ (signal+background), depending on whether you are questioning the presence of 
a signal, or measuring the size of a signal you know to be there. Could a similar unified 
approach be found here? 

The unified approach prevents the experimenter flip-flopping between limit 
strategies, but still leaves the choice of confidence level free. An experimenter who 
normally used 90% limits may choose to present 68% limits if they obtain a low result 
where they believe there should be a signal, or 95% or 99% limits if they obtain a high 
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result where they beheve no signal should be. Is there some natural way of incorporating 
this (undesirable?) freedom in an algorithm? 
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8.2. Statistical issues in heavy flavour physics 

Glen Cowan, Royal Holloway, University of London 



8.2.1. Introduction In this paper two statistical issues are addressed that come 
up in many physics analyses but which are particularly relevant for heavy flavour 
physics. The first, presented in Section |8.2.2| , is the problem of combining measurements 
of the branching ratios of rare decays where some of the measurements are based 
on only a few or perhaps no candidate events. The solution recommended is a 
straightforward application of the method of maximum likelihood. The next question, 
discussed in Section ^.2.3| , is much more open-ended, namely, how to quantify theoretical 
uncertainties. Both issues are of course not unique to heavy flavour physics, and the 
approaches recommended can be applied in many other contexts. 



8.2.2. Combining branching ratios for rare decays Consider first a single 
experiment in which one searches for rare decay and finds n candidate events out of 
total decays. Suppose the efficiency for the signal is e and there are h expected 
background events. In general, n should be treated as a binomially distributed variable 
where represents the total number of trials and the binomial probability B is the 
branching ratio. If B is very small and N very large, however, we can treat n as 
following a Poisson distribution with expectation value i.e., 

P(n; z/) = ^e"'' . (115) 
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Here the number of events observed n is the sum of Us signal and n?, background events 
and the expectation value of n, here called is also the sum of the signal and background 
expectation values s and h: 



n = Us + Ub , 
u = s + b . 



(116) 
(117) 



The branching ratio B, which is the parameter we want to estimate, is related to the 
expected number of signal events s, the efficiency e and the total number of events N 
by 

e = (118) 

For a single experiment, the likelihood function L{B) is simply the Poisson 
probability ( |115| ) evaluated with the observed number of events n and expressed as 
a function of the branching ratio B. Maximizing L gives the maximum likelihood (ML) 
estimator B, 

b 



B 



n 



eN 



(119) 



where here estimators will be denoted by hats. The statistical uncertainty in B can be 
quantified by giving an estimate of its standard deviation. 



(Jr. 



In 
eN 



(120) 



or by reporting a confidence interval (see e.g. |I[]). 

Now consider the case where we have m independent measurements of the 
branching ratio B. The m experiments report ni,...,nm candidate events out of 
Ni, . . . , Nm total decays; they have signal efficiencies ei, . . . , and expected numbers 
of signal and background events si, . . . , Sm and bi, . . . ,bm- Assuming the measurements 
are statistically independent, the likelihood function is the product of the Poisson 
probabilities for each. 



u 



i=l 
m 

n 

i=l 



[EiNiB + bi 



-(e^N,B+b,) 



:i2ii 



Hi 



The value of B that maximizes (|121|) is the ML estimator i3, which in practice will need 
to be determined numerically, e.g., by solving 



(91nL 
dB 



i=l L 



. 



(122) 



InL 
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As in the case of a single measurement, the statistical uncertainty in B can be 
quantified by estimating its standard deviation using, for example, the curvature of the 
log-likelihood function at its maximum, 

-1 

(123) 

B=B/ 

or by constructing a confidence interval. 

As the final result must be determined numerically it is difficult to see what the 
properties of the ML estimator are. If, however, we consider the case where there is no 
background, i.e., hi — Q for all i, then the estimator B takes on the simple form 

for which the standard deviation can be estimated by 



= Y^m ^ ■ (125) 

For this case, it is easy to see that there is no problem if some of the n,i are small 
or even zero. Such problems would arise, for example, in least squares averaging of 
individual estimates of B. Furthermore, in the real case with nonzero backgrounds, the 
ML estimate is not spoiled by including experiments that carry less information than 
others because, for example, they have larger backgrounds or are based on smaller data 
samples. 

It should not be controversial that the ML estimator has many advantages and 
should be used unless there is a good reason not to. The important point to note 
here is that it is not sufficient for experimenters to report only an estimate of the 
branching ratio and a standard deviation. In order to combine the measurements using 
the method of maximum likelihood, one requires the number of observed events, the 
efficiency, the expected number of background events and the total number of events 
for all experiments. 



8.2.3. Quantifying theoretical uncertainties In many experimental analyses, 
the dominant uncertainty in the determination of a parameter does not stem from 
measurement errors but from the approximations used in deriving the theoretical 
prediction for the quantities observed. An important example from heavy flavour 
physics is the theoretical prediction for the differential decay rates of B mesons into 
various final states. Although the Standard Model makes a well-defined prediction for 
the functional form of the probability distribution for measured decay times, the exact 
form of this prediction cannot be derived because of calculational difficulties, especially 
for multi-body final states. Instead, approximations are used, based, for example, on 
a perturbation series calculated only to some finite order. As a result, the estimated 
parameter values are in general uncertain even in the absence of measurement errors. 
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statistical or otherwise. This is what is meant here by the 'theoretical uncertainty' in 
an estimated parameter. 

Here we will argue that the Bayesian approach to statistics provides a potentially 
useful framework for quantifying theoretical uncertainties. To keep the discussion 
as general as possible, suppose we measure a vector of data x, whose probability 
distribution function (pdf) is predicted as a function of a vector of parameters d. For 
example, the vector x could represent a set of measured decay times for B meson decays 
into J/V'Kg and the vector d would then be related to CKM matrix elements as well as 
other Standard Model parameters that may enter. 

Now suppose that we do not know the exact form of the pdf /(x; 6) but the function 
that we have to work with is only an approximation. This approximation can be viewed 
as belonging to a family of functions, one of which is the true prediction of the theory, or 
at least close to it. Suppose wc can parametrize this family by some vector of nuisance 
parameters /j,. The elements of (j, could represent, for example, higher order coefficients 
in a perturbation series that are not included in the prediction actually available, but 
which could be calculated in principle. It is necessary in the procedure presented here for 
one of the members of the family of functions to be sufficiently close (in principle, equal) 
to the true prediction. The family can always be made arbitrarily large by including 
more nuisance parameters, but this will have the end result of inflating the uncertainties 
in the parameters of interest. 

In Bayesian statistics, all knowledge about the parameters of a theory is summarized 
by giving the posterior pdf, p(0,/x|x). This is the conditional pdf for the parameters 
given the data x. In our case, the parameters consist of those in which we are interested, 
0, as well as the nuisance parameters fi. According to Bayes' theorem, they are related 
to the likehhood function l/(x|0, //) and the prior pdf 7r(0, fx) by 



The likelihood function L(x|0,/Lt) is the conditional probability for the data given the 
parameter values, and the prior pdf 7r{6, /*) summarizes our knowledge of the parameters 
before carrying out the experiment. 

It should be emphasized here that the probability associated with the value of 
a parameter cannot be interpreted meaningfully as a frequency of an outcome of a 
repeatable experiment. Instead it is understood to reflect the degree of belief that 
the parameters have particular values, an interpretation called subjective probability. 
In order to determine this degree of belief given the data, we need to say what it 
was before carrying out the experiment, i.e., we need the prior pdf 7r(0, /Li). Bayesian 
statistics provides no fundamental prescription for determining the prior pdf. Given 
n{d,fj,), however, Bayes' theorem specifles how our knowledge about the parameters 
should change in the light of the data. 

The usefulness of the Bayesian approach in the present problem will depend on 
our being able to parametrize the pdf /(x;0,/Lt) in terms of the nuisance parameters 
/J,, and to make a reasonable statement about their probabilities. Often this means 




(126) 
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interrogating the theorist to the point where he or she tells you how large, say, the 
next order coefficient in a perturbation series could be. At ffist, the theorist may say 
something like 'Until I've calculated it, all values are equally likely'. Under pressure, 
however, it will be admitted that if it were to differ by many orders of magnitude 
from the previous term, or from one's naive expectation of the value based on past 
experience with perturbation theory or Pade approximants or some other procedure, 
then this would be surprising. 

The key now is to quantify this level of surprise, by calibrating it against something 
for which probability can be defined in terms of a frequency. The theorist should state 
the range for the nuisance parameter such that the degree of surprise at finding out that 
it is outside this range is equal to, say, the level of surprise that one would experience at 
drawing a white ball out of an urn of balls, 32% of which are white. In most cases this 
level of surprise will be a smoothly varying function of the parameter, which could be 
represented by some sort of a bell-shaped curve like a Gaussian. It would be rare, for 
example, that a uniform distribution between fixed limits would provide a reasonable 
model for quantifying the uncertainty in such a parameter. Assuming a Gaussian model, 
the central interval containing 68% of the probability corresponds to the one-sigma error. 
In Bayesian parlance these are often called 68% credible intervals. More information on 
subjective probability and Bayesian statistics can be found in 0; its application to 
problems of data analysis in physics is treated at various levels of detail in p, ||, ^ ^. 

It is clear that this calibration of a degree of belief will never be perfect, and we 
should be happy if it is reasonable to within a factor of two or even a factor of ten. 
In principle the level of calibration could be measured by forcing theorists to guess the 
68% credible intervals, waiting long enough for the calculations to be made and then 
determining the frequency with which they are right. An alternative proposed at this 
workshop would be to force them to bet (large) sums of money on the values of nuisance 
parameters that at some point will be calculated. The Bayesian framework will be useful 
if reasonable prior probabilities can be assigned; there are no doubt cases in which this 
is not possible. 

The importance of making reasonable statements about prior probabilities holds as 
well for the model parameters 6, not only the nuisance parameters ^i. Here the posterior 
will be insensitive to the prior pdf for d as long as it is a more slowly varying function 
of 6 than is the likelihood function. If this is not the case, then the Bayesian approach 
is less likely to yield useful results, since the posterior pdf will be determined more by 
the prior beliefs than by the data. 

Providing one can quantify the necessary prior probabilities, it is a straightforward 
exercise to apply Bayes' theorem to obtain the posterior pdf p{d, /x|x). We are then in a 
position to apply one of the most powerful tools of Bayesian statistics: we eliminate the 
nuisance parameters from the problem simply by computing the marginal pdf p{0\x.), 
i.e.. 




(127) 
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One can then summarize the posterior pdf by, for example, presenting contours 
of constant probabihty in two-dimensional subspaces or by characterizing it by various 
sorts of location parameters such as the mode. Regions in the parameter space can be 
determined which contain specified probabilities, e.g., 68%, 90%, etc. 

The sizes of the contours of constant probability reflect both the statistical and 
systematic errors. These two components of the uncertainty can be separated by 
constructing estimators 0, e.g., by taking the mode of the posterior pdf. These will 
be specific functions of the data, and as such they are random variables in the sense of 
frequentist probability. One can then ask, for example, how would be distributed upon 
repetition of the experiment, and in this way the statistical and systematic components 
can be separated, if desired. 

Although up to now we have only taken /x to quantify errors due to approximations 
made in deriving theoretical predictions, it can just as easily represent other types 
of systematic uncertainties. In either case, the assumptions made in quantifying the 
uncertainties should clearly be reported along with the result. It should be emphasized 
that the posterior pdf by itself does not in general provide a useful summary of the 
result of the experiment, as it combines both the data and the subjective beliefs needed 
to construct the prior pdf. If only the posterior pdf is reported, it is impossible for a 
consumer of the result to separate the infiuence of the experimenter's prior beliefs from 
the infiuence of the data. 

The experimenter should therefore separate as well as possible a summary of the 
observations from their interpretation. In a paper, for example, this may be done in 
sections titled 'results' and 'discussion' or 'interpretation'. In the section on results one 
can report the values of observables that to the greatest reasonable degree are free of 
systematic effects, theoretical or otherwise. For example, one should report the various 
numbers of events observed that satisfy specific criteria, and not only the estimates 
of parameters that are derived from these data. In the interpretation of the results, 
however, there is no reason for experimenters to shy away from making reasonable 
(stated) assumptions based on subjective probability and showing what they imply for 
the values of parameters of interest. 

8.2.4- Conclusions A straightforward technique based on the method of maximum 
likelihood has been presented for combining measurements of branching ratios. In order 
to employ this method, experimenters must report the number of candidate events seen, 
the efficiency of the selection procedure, the expected background and total number of 
decays observed. It cannot be used in general if only an estimate of the branching ratio 
and its standard deviation arc reported. 

The use of Bayesian statistics in quantifying systematic errors, in particular those 
due to uncertainties in the functional form of theoretical predictions, has been discussed. 
In this framework one must quantify the uncertainty of unknown constants using 
subjective probability, i.e., probability based on a degree of belief. The method is 
useful to the extent that the uncertainty can be parametrized and quantified, e.g., by 
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interrogation of the theorist responsible for the prediction. Owing to difficulties in 
calibrating one's degree of belief relative to frequentist probabilities, this is necessarily 
an approximate procedure and all assumptions made in quantifying the uncertainty 
should be reported along with the result. 
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